Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,782,913 A→C 22.0% V250G (GTG→GGG)  gpsA ← glycerol‑3‑phosphate dehydrogenase (NAD+)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,782,9130AC22.0% 96.7 / 10.0 73V250G (GTG→GGG) gpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base C (3/14);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 2.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGC  >  NC_000913/3782678‑3783149
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
tAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCa                                                                                                                                                                                                                                              <  2:1162383/239‑1 (MQ=255)
        tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCCCC‑ACCGGGTGTCGG‑GTCCGGGAAGCACCAGaaccccc                                                                                                                                                                                                                                      >  1:1037945/1‑233 (MQ=255)
        tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑CCCGGTTGTCGG‑TTCCGGGAAGCACCACATCG‑ccc                                                                                                                                                                                                                                     >  1:505674/1‑238 (MQ=255)
        tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCG‑CCCGGGTTTTGG‑TACCGGTAAACACCCGATTG‑ccc                                                                                                                                                                                                                                     >  1:616706/1‑238 (MQ=255)
         cAATGCTGCGTGGCGGTCTTTTTTCCCGAAATATAAACCTTGAAAATTTCCCTGGTTTATTGGCATTCCAACGCCG‑AGCCGTGCCCCCGTTTCGCGGCTTTCTTTCGTATTGCGGTACCCTCCCCCCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACACCCAGGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCCCCAGATCG‑cc                                                                                                                                                                                                                                      <  2:616706/237‑1 (MQ=255)
         cAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCa                                                                                                                                                                                                                                     <  1:812018/238‑1 (MQ=255)
            tgctgcCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCC‑ACTGGTTGTCGG‑GACAGGTAAGAACCCAAACG‑CCa‑aac                                                                                                                                                                                                                                 >  1:484379/1‑236 (MQ=255)
                   ccggccgtgTTTTTCCCGCAATATAATCCTTGAAAAATTTCCTGTGTTTTTGGCATTTCAGCGCCGAACCGATGCCCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCCCCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGGGCATCATGCCAAAACGGCGGTTCCGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCAt                                                                                                                                                                                                                         <  2:944065/231‑1 (MQ=255)
                        cgcgTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGGTTAGCGCCCGGGTTTTGGG‑GACGGGAAAGCCCCAAAAGC‑CCA‑AACCCCCCCATGccc                                                                                                                                                                                                                     >  2:147341/1‑239 (MQ=255)
                           gTTTTTCCCGAAATAAAATCCTTAAAAATTTTCCTCGGTTATTGGCATTTCAACTCCGAAGCGATGCCCCAG‑TTCGCGGCCTTCTTTCGTATGGCGGTAGCCTTCCCCCCCCGGACCAATCTTCCCCGGCGCGCTTTGTACATCCATGCCCTGCCCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATaa                                                                                                                                                                                                                 <  2:1037945/239‑1 (MQ=255)
                                       atatAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGGAAGAACCAGATAG‑CCC‑AGCCCCCCCCCGCCCCTACAAGTGTGAGGGt                                                                                                                                                                                                     >  2:155776/1‑239 (MQ=255)
                                        tataATACTTGAAAATTTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCAt                                                                                                                                                                                                                         <  1:1051276/218‑1 (MQ=255)
                                        tagaataCTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCAt                                                                                                                                                                                                                         <  1:536292/215‑1 (MQ=255)
                                         ataatcCTTGATAATTTCCCCCGTTTTTTGACATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTACCCTTCCCCCCCCTGCCCAATCTTCTCCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCCCCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑Cg                                                                                                                                                                                                  <  2:484379/239‑1 (MQ=255)
                                            acccTTGATAATTTCCCATGGTTTTTGGCTTTACAACGCCGAAGCGATGCCCCTG‑TTCGCGGACTTCTTTCGTTTTGCGGTAGCCTTCCCCCACCGGCCCAATCTTCCCCGGCGCGTTTTGTACCTCCATCCCCTGCCCGGGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAAACCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCa                                                                                                                                                                                               <  2:669049/236‑1 (MQ=255)
                                              tcTTGATAATTTCCCTCGGTTGTTGCCATTTCACCCCCGAAGCGATCCCCCAG‑TTCGCGGCCTTCTTTCGTATTGCGGTAGCCTTCCCCCCCCTGCCCAATCTTCCCCTGCGCGCTTTGTACACCCATGCCCTGCCCGACCATCATGCCAAACCGGCGGTTCCGCG‑CTTGTTTGTCGG‑TACAGGTAAGCCCCAGATCC‑CCA‑ACCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCAcc                                                                                                                                                                                             <  2:176200/238‑1 (MQ=255)
                                              aCTAGATAATTTGCCTCGGTTTATGGCGTTTCACCGCAGAACCGATGCGCCAG‑TTGGCGGCCTTCTTTCGTATTGCGGTACCCTCCCCCCCCCTGCCCAATCGTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCAGCCCAAAACGGCGGTTACGCA‑ACTGGTTGTCGG‑TACGGGTAAGCACCAGACCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCAcc                                                                                                                                                                                             <  2:176547/239‑1 (MQ=255)
                                                ttGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCC‑ACTGGGTGTCGG‑TACAGGGAAGCACCCGATCG‑CCA‑AGACCCCCCATTCCCCCAAAGGGGGCAGGGT‑CGGCAccc                                                                                                                                                                                            >  1:837181/1‑238 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑CCTGGTTGTCGG‑GTAAGGGTAGAACCCAAACC‑CCACACCCCCCCCATCCCCCAAAAGG‑GGCAGGGT‑CGGGAccccc                                                                                                                                                                                          >  1:1022929/1‑236 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑CCTGGTTGTCGG‑GACAGGGAAGCCCCCAAACG‑CCA‑AACCCCGCCATGCCCCTACAAGTGGCAGGGG‑CGGCACCCac                                                                                                                                                                                          >  1:425071/1‑237 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTTTCGG‑TACAGGGAAGCACCAGATCG‑CCC‑AACCCCGCCCTGCCCCTAAAGGGGGACGGGG‑CGGGAcccca                                                                                                                                                                                          >  1:176547/1‑236 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AACCCCGCCCATCCCCTAAAGGTGGCAGGGT‑CGGCACCCAgc                                                                                                                                                                                         >  1:858266/1‑239 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACCTGGTGTCGG‑GACAGGTAAGCACCAGAAAG‑CCA‑AACCCCCCCATGCCCCTAAAGGGGGCCAGGG‑CGGGACCCAgc                                                                                                                                                                                         >  1:632862/1‑239 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACCGGTTGTCGG‑TACAGGGAAGCACCAAATCG‑CCA‑AGCCCCGCCACGCCCCCAAAAGGGGCAGGGT‑CGGCACCCAgc                                                                                                                                                                                         >  1:318431/1‑239 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCG‑ACTGGTTGTCGG‑GTCAGGTAAGGACCAGAACT‑CCC‑AACCCCGCCATGCCCCCCAAGAAGGCAAGGT‑CGGccgcccc                                                                                                                                                                                          >  2:468665/1‑232 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCG‑ACTGGTTGTCGG‑GACAGGGAAACCACAAAACG‑CCA‑AACCCCCCCATGCCCATAAAAGTGGCAGGGT‑CGGCACCCCgc                                                                                                                                                                                         >  1:283221/1‑239 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTAAGCG‑ACTTGTTTTCGG‑GACAAGGAAGCACCCGAATG‑CCA‑AACCCCCCCCAGGCCACTAAAGTGGCAGGGG‑CGGCACCCac                                                                                                                                                                                          >  1:669049/1‑237 (MQ=255)
                                                  gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGGTACGCG‑ACTGGGTGTCGG‑TACAGGTAAGAACCAAAACG‑CCA‑AGACCCGCCATGGGCCTAAAAGTGGCAgggggcgggccccc                                                                                                                                                                                           >  2:365181/1‑227 (MQ=255)
                                                   aaaaaTTTCCGGGGTTATTGCCATTCAAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCAGGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAgcg                                                                                                                                                                                        <  2:425071/237‑1 (MQ=255)
                                                         ttCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCg                                                                                                                                                                                                                             <  1:802272/197‑1 (MQ=255)
                                                                  tATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAACCGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGAcgcg                                                                                                                                                                         <  1:365181/239‑1 (MQ=255)
                                                                           ttCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCG‑CCTGGTTGTCGG‑GGCAGGGAAGGACCAGAAGGCCCA‑ACCCCCCCCAGCCCCCTAAAGGGCGGAGGGTGGGGCCCCCAGCGCCCCCCAAAGCcgcacctact                                                                                                                                                                   >  2:1133558/1‑231 (MQ=255)
                                                                            tCACCGCCGAAGCGATCCGCAAT‑TCCCCGGCTTTCTGCGGTATTGCGGTAGCCTTCCCCCACCGGACCAATCTTCTCCTGCGCGTTTTGTACATCCATGCCCGGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TAAAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                               <  2:505674/239‑1 (MQ=255)
                                                                                       gCGATGCCCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACgg                                                                                                                                                     <  1:159351/238‑1 (MQ=255)
                                                                                                     tCGGGGACTTCTTTCGTATTGCGGTAGCCTTCCACCCCCTGACCAATCTTCCCCGGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAgcg                                                                                                                                           <  1:1241212/235‑1 (MQ=255)
                                                                                                     tCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAgcgc                                                                                                                                          <  1:434305/236‑1 (MQ=255)
                                                                                                         ggATTTCGTGCGTATTGGGGTAGCCTTCCACCCCCTGACTAATCTTCCCCTGCGCGTTTTGAACTTCCAGGCCCGGACCGAGCATCAGGCCAAAACGGCGGTTACGCG‑ACGGGTTGTCGG‑TACAGGTAAGCACCAGACCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTAcg                                                                                                                                    <  2:862122/238‑1 (MQ=255)
                                                                                                                    gTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCAAAAACGGCGGTTCCGCG‑ACTGGTTGTCGG‑TACAGGTAAGCCCCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCAcc                                                                                                                        <  1:468665/239‑1 (MQ=255)
                                                                                                                     tATTGCGGTAGCCTTCCCCCCCCTGACCAATCTTCCCCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCCCCAGATCG‑CCA‑AGCCCCGCCATGCCCATaaa                                                                                                                                                                                                                <  2:488724/151‑1 (MQ=255)
                                                                                                                     tATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTATCC‑CCCGGGTGTCGG‑TACAGGTAAGCACCAAAACG‑CCA‑AACCCCCCCCTGCCCCTCAAAGGGGCCGGGG‑CGGGACCCAGCGCCGCACCCACACGCGACCTTTTATCCCGCCCACCGGCGGTCAAGGCCGATCgcgccgtcgtcccc                                                                                                                       >  2:1010800/1‑229 (MQ=255)
                                                                                                                           ggTAGCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCAcc                                                                                                                                                                                             <  1:330607/163‑1 (MQ=255)
                                                                                                                               gCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTTTCGG‑GACAGGGAAGCACCAGAACG‑CCA‑AACCCCCCCCTGCCCCTAAAAGGGGCAggg                                                                                                                                                                                                      >  1:488724/1‑151 (MQ=255)
                                                                                                                                ccTTCCACCACCGGACCAATCTTCCCCGGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATg                                                                                                                                                                                                                        <  2:259762/132‑1 (MQ=255)
                                                                                                                                    ccacccCCTGCCCAATCTTCCCCGGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGGCGG‑TACAGGTAAGCACCAGACCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAgcg                                                                                                                                           <  1:438757/204‑1 (MQ=255)
                                                                                                                                         cccTGACCAATCTGCCCCTGCGCGCTTTGTACATCAAGGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACGGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGa                                                                                                      <  1:529693/235‑1 (MQ=255)
                                                                                                                                          ccTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGGTGTCGG‑GACAGGGAAGCACCAGATCG‑CCA‑AGCCCCCCCCTGCCCCTAAAgg                                                                                                                                                                                                              >  1:259762/1‑132 (MQ=255)
                                                                                                                                            tGACCAATCTTCCCCTGCGCCTTTTGTCCATCCATGCCCTGACCGAGCATCATCCCAAACCGGCGTTTCCGCG‑ACTTGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCTCCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATccc                                                                                                <  1:995961/239‑1 (MQ=255)
                                                                                                                                              cccACTCTCCACCTGCGCGCTTTGTCCATCCATGCCCTAACCGGGCTGCATGCCAAAGCGGCGGTCCCGCG‑ACTGTTTTCCGG‑TACGGGTAAGCCCCAGACCC‑CCA‑GGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCCCCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCg                                                                                              <  2:318431/238‑1 (MQ=255)
                                                                                                                                              aCCAATCTGCCCCGGGCCGTTTTGTACATCCACTCCCTACCCGACCATCATGCAAAAGCGGCGGTTCCGCG‑ACTGGTTGTCGG‑TACAGGTAAGCCCCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCg                                                                                              <  2:285254/239‑1 (MQ=255)
                                                                                                                                                 aaTCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCGCG‑ACCTGTTTTCGG‑GACAGGTAAACACCAAAACG‑CCA‑AGCCCCGCCCTGCCCCTAAAAGGGGCACGGG‑CGGGCCCCCGCGCCCCCCCCACACGCGGAACTTCAGCCAGACCACGCGCGCTCAACCCCCGACGCGCACTCGCACCAAAACACAAACCGTCGGGCagtcccaca                                                                                            >  2:687014/1‑230 (MQ=255)
                                                                                                                                                   ttttctCCTGCGCGTTTTGTCCATCCATGCCCTGCCCGACCATCATGCAAAACCGGCGGTTCCGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCAcc                                                                                                                        <  1:824167/206‑1 (MQ=255)
                                                                                                                                                   tcttctCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCAGGCAAAAGCGGGGGTTACGCG‑ACTGGTTGTGGG‑TACAGGTAAGCCCCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                               <  1:588800/169‑1 (MQ=255)
                                                                                                                                                                         aTCCATGCCCTGACCGAGCATCAGGCCAAAACGGCGGTTACGCG‑ACGGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcg                                                                   <  1:1133558/239‑1 (MQ=255)
                                                                                                                                                                           ccATGCCCTGACCGAGCATCATGCCAAACCGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACcgcg                                                                 <  1:888145/239‑1 (MQ=255)
                                                                                                                                                                           ccATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCCCCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGa                                                                                                                                                                             <  1:989866/131‑1 (MQ=255)
                                                                                                                                                                                 ccTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCCCCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCa                                                           <  1:147341/239‑1 (MQ=255)
                                                                                                                                                                                     aCCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑GACAGGGAAGCACCAAAACG‑CCC‑AGCCCCCCCATGCCCCTAAAGGGGGCCGGGG‑CGGCCCCCCGCCCCCCCCCCAAACGCGGCAttt                                                                                                                                                                   >  2:989866/1‑131 (MQ=255)
                                                                                                                                                                                     aCCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGGTGTCGG‑GACAGGTAAACACCAAAACG‑CCC‑AACCCCCCCCCGCCCCTAAAAGGGGGAGGGG‑CGGCACCCCGCGCCCCCCCAAGAACCCACCTTTTAGCCCCCCCCCCGGGGAGCAGCGCCGGACGCGCCTTCGCACCCAAACCGATACCCCCGGACCTCCCCGCCCCCCTGGCCACAAAGTATTTTATCGCgccgccccgc                                                        >  2:886736/1‑234 (MQ=255)
                                                                                                                                                                                      cggAGCATAACGCCAAACCGGCGGTTCCGCG‑CCT‑GTTGGCGGTTACAGGTAACCCCCAGCCCG‑CAA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg                                                      <  2:858266/237‑1 (MQ=255)
                                                                                                                                                                                       cGAGCATCATGCCAAAACGGGGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg                                                      <  1:1150249/238‑1 (MQ=255)
                                                                                                                                                                                         aGCATCATGCCAAAACGGCGGTTACGCG‑GCTGGTTGTCGG‑GACAGGTAAGCACCAGAACG‑CCA‑AGCCCCCCCCTGCCCCTAAAGGGGGCAGGGT‑CGGCACCCAACGCCCCCCCCAGACGCGACATTTCAGCCAGCCCCCGGGTGATCAGCGCCCTACCCGCCCTCGCACCCAAACAGCCAACCTCACAAATCCCCTCCACAATAGCCATAACGGATTTTACAGCGCCGCCaccagtca                                                    >  2:847472/1‑231 (MQ=255)
                                                                                                                                                                                          gCATCATGCAAAAACGGCGGTTACGGG‑ACTGTTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTg                                                      <  1:1176027/235‑1 (MQ=255)
                                                                                                                                                                                           catcatGCCAAAGCGGCGGTTCCGCG‑ACTGGTTGCCGG‑TACAGGAAAGCCCCAGACCC‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACg                                                  <  1:513776/238‑1 (MQ=255)
                                                                                                                                                                                               atGCCAAACGGGGGTTTCCGCG‑CCTGGTTGGCGG‑TACAGGTAAGCCCCAGATCC‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAt                                             <  1:155776/239‑1 (MQ=255)
                                                                                                                                                                                                         ggccttcccgcg‑TCTGGTTGTCGG‑AACAGGAAAGCACCAGACCG‑CCA‑ACCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCCCCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGa                                   <  2:283221/231‑1 (MQ=255)
                                                                                                                                                                                                                          gTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGa                   >  1:5143/1‑239 (MQ=255)
                                                                                                                                                                                                                                     cAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGcc         >  1:1141603/1‑239 (MQ=255)
                                                                                                                                                                                                                                     cAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGc          >  1:1228092/1‑238 (MQ=255)
                                                                                                                                                                                                                                      gggTAAGCCCCAGACCC‑CCC‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg        <  1:1075470/238‑1 (MQ=255)
                                                                                                                                                                                                                                      aGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCgccgcc                                                            >  2:1054638/1‑187 (MQ=255)
                                                                                                                                                                                                                                      aGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCgccgcc                                                            <  1:1054638/187‑1 (MQ=255)
                                                                                                                                                                                                                                           aGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGCCATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGt    >  1:542455/1‑238 (MQ=255)
                                                                                                                                                                                                                                             cACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgc  >  2:486184/1‑238 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
TAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAG‑TTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCG‑ACTGGTTGTCGG‑TACAGGTAAGCACCAGATCG‑CCA‑AGCCCCGCCATGCCCATAAAGGTGGCAGGGT‑CGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGC  >  NC_000913/3782678‑3783149

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑