Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,791,798 A→C 21.8% T176P (ACC→CCC)  gabD → succinate‑semialdehyde dehydrogenase I, NADP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,791,7980AC21.8% 100.0 / 2.9 69T176P (ACC→CCC) gabDsuccinate‑semialdehyde dehydrogenase I, NADP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (28/23);  new base C (7/11);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAA  >  NC_000913/2791563‑2792025
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                    
ccgggtgaaaCCTCTGGCCGAAGCGAAAGGCGAAAACAGCTCCGCCGCCTCCTTGTTTGGGTGGTTTACCGAAGAAGGCAAACGCTTTTAGGGCGACACCATCCCTGGTCATCACGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCa                                                                                                                                                                                                                                   <  2:134727/232‑1 (MQ=255)
    ggTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATggg                                                                                                                                                                                                                               >  1:84460/1‑238 (MQ=255)
    ggTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGt                                                                                                                                                                                                                               >  1:682535/1‑239 (MQ=255)
    ggTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCACCAGCCTGCACCATGGt                                                                                                                                                                                                                               >  1:733838/1‑239 (MQ=255)
    ggTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCGGGCCGCAGGCGACACCATgtc                                                                                                                                                                                                                               >  1:478375/1‑237 (MQ=255)
    ggTAAACCACGGGCCGAAGCGAAAGGCGAAATCAGCTCCGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCa                                                                                                                                                                                                                                   <  1:365064/235‑1 (MQ=255)
     ttAAACCCCGGGCCGATACTGAAGGTGAAATCAGCCTCGCAGCCTCGTTATCTGAGTGGGTGGCCGAAAAAGGCATACGCATTTATGGGGACACCTTTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCGCCGCGGCTATCATGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATg                                                                                                                                                                                                                                 <  2:347547/235‑1 (MQ=255)
     gTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGGTGATTAACCGCAAAGACGGTCCGGCCCTGGCACCAGGCCGCCCCACggg                                                                                                                                                                                                                               >  2:480837/1‑237 (MQ=255)
         aCCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGa                                                                                                                                                                                                                          >  1:229290/1‑239 (MQ=255)
              tGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGCTTGGCGTCACCGCGGCTATCCCGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAACCGGTCCGGCCCTGGCAGCCGGCCGCCCCAGGGGGccaaaacc                                                                                                                                                                                                                      >  2:325491/1‑231 (MQ=255)
              tGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGGTGCACCATGGTGGTGAAGcc                                                                                                                                                                                                                      >  1:443620/1‑238 (MQ=255)
              tGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGcc                                                                                                                                                                                                                      >  1:203771/1‑238 (MQ=255)
              tGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGACGCAGGCTGCACCATGGTGCTGAAccc                                                                                                                                                                                                                      >  1:439766/1‑238 (MQ=255)
              tGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGCCAGCAGGCCGCCCCATGGGGGTGaggggg                                                                                                                                                                                                                     >  2:365064/1‑234 (MQ=255)
              tGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCGGGCAGCAGGCTGCCCCATGGGGCTGAAGcc                                                                                                                                                                                                                      >  2:513012/1‑238 (MQ=255)
              tGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGAATACCCGCAAAACCGGTCCGGCGCTGGCAGCAGGCTGCCCCATGGCGCTTAAccc                                                                                                                                                                                                                      >  1:755201/1‑238 (MQ=255)
              tGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCGGATTGTTATCAAGCAGCCGATTGGCCTCACCGCGGCTATCACGCCGTGGAACTTCCCGGGGGCGAGGATTACCCGAAAAGCCGGCCCGGCCCGGGCAGCGGCTGGAACCAGGTGGTGaaaacc                                                                                                                                                                                                                     >  2:703321/1‑234 (MQ=255)
                    aaGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCCCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGt                                                                                                                                                                                                               >  1:630802/1‑239 (MQ=255)
                        gAAAGGCGAAATCGGCTACGCCGCCTCCTTTATGGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTAATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAg                                                                                                                                                                                                            <  1:185944/238‑1 (MQ=255)
                                            ccccttcctTTATTGGGTGGTTTGCCGAAGAAGGCAAACGCTTTATGGGCGACACCATTCCTGGTCATCAGGCCGATAACCGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTgcgctggcg                                                                                                                                                                                       <  1:703321/234‑1 (MQ=255)
                                                     tttttGGGTTGTTTGCCAAAGAAGGCAACCGCTTTTTTGGCGACCCCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAg                                                                                                                                                                              <  2:14139/236‑1 (MQ=255)
                                                          gAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGCCGCCGTTCTCTGCCCCTGCGCTGGCGGAACTGGc                                                                                                                                                                         >  1:453668/1‑239 (MQ=255)
                                                           aGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGTGGCTGAATCCCCCCAGTCAGACCCCGTTCCCTCCGCTGGCGCTGGGGGGGCTGGcc                                                                                                                                                                        >  1:892867/1‑238 (MQ=255)
                                                           aGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGGGCTGAAGCCCGCCAGTCAGACGCCGTTCTCCGCGCTGGCCCTGGGGGAggcgggg                                                                                                                                                                        >  2:295768/1‑233 (MQ=255)
                                                             tGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCTTCACCCCGGCTCTCCCGCCGTGGAACTTCCCGGCGGCGATGATTCCCCGCCAAGCCGAGCCGGCGCTGGCAGCCGCCTGCACCATGGTGCTCAAGCCCCCCACTCAGCCGCCCTTCCTTGCGCTGGCCCTCTCCGAGGTGGccg                                                                                                                                                                       >  2:90836/1‑236 (MQ=255)
                                                             tGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCGGGCGCTGCCAGCAGGGCGCCACCAGGGCCTGAAGCCCGCCAGCCCAACCCCGCTTCCTGCGCTGGCCCCGGCGGGGCCgggca                                                                                                                                                                       >  2:68905/1‑235 (MQ=255)
                                                             tGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCCCCCTGGGGCTCAAGCCCGCCAGTCAGACCCCGTTCTCTGCGCTGGCCCGGCCGGAGCGGGCgg                                                                                                                                                                       >  1:846854/1‑237 (MQ=255)
                                                             tGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCCGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGGGCTggggg                                                                                                                                                                       >  2:594012/1‑235 (MQ=255)
                                                             tGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCCGTCAGACGCCGTTCTCTGCGCTGGCGATGGCGGAGCTGGCGa                                                                                                                                                                       >  1:702911/1‑238 (MQ=255)
                                                             tGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGCGGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGGCAGACGCCGTTCTCTGGGCTGGGGCggggggggggggcgg                                                                                                                                                                       >  2:430903/1‑229 (MQ=255)
                                                             tGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGGCCGGCGCTGGAACCAGGCTGGACCATGGAGCTGAAGGCAGCCCAGCAAAACCCATTCTCTCCCCTGGGGGAGGCGGAGCCGGGGag                                                                                                                                                                      >  2:155974/1‑238 (MQ=255)
                                                              tgTTTGCCGAAGAAGGCAATCGCTTTTATGGGCCCACTTTTCCTGGTAATCAGGCCGATAAACCCCTGGTGTTTACCAAGCCGCCAATTGGCTCACCCGGGGCTATCACGGCGGGGACCTTCCCGGTGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCCCCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGAt                                                                                                                                                                      <  2:453668/237‑1 (MQ=255)
                                                                 ttGCCGAAGAAGGAAAACGAATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCCCGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCcg                                                                                                                                                                   <  1:513012/238‑1 (MQ=255)
                                                                  tGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGCTCTCTGCGCTGGCGCTGGGGGAGCCGGCGAACcggc                                                                                                                                                                 >  2:314586/1‑237 (MQ=255)
                                                                     cGAAGAAGCAAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCgg                                                                                                                                                              <  1:117483/239‑1 (MQ=255)
                                                                                 agcaaTTTAGGGCGACCCCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCCCGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTg                                                                                                                                                  <  1:480837/235‑1 (MQ=255)
                                                                                  ccggtttttgggCGACACCTTCCCTGTTCATCCGGCCAAAAAACGCCTGTTTGTTACCAAGCGGGCGGTGGGGGCCGCCGGGGTTATCCCGCCGTGGAGCTCCCCGGGGGCGGTGATTACCCGGAACGCGGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTg                                                                                                                                                  <  2:210146/229‑1 (MQ=255)
                                                                                   ccttttATGGCCACCTCTTTCCGGGCCACGAGGGCGAAAAACGCGTGTTTTTTAACAACCGGCCGATGGGCGTACCCGCGGCTATCCCGCCGTGGAACTTCCCGGCGGCGGTGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTggg                                                                                                                                                <  2:682535/236‑1 (MQ=255)
                                                                                    cctttttttgggACCCCCTTCCTGGTTCGCTGGCCGCAAAACCGCCTGTTTTTTATCAGGCGCCCGATGTGCGTCCCCGTGGCTCTCACGCCGTGAAACTTCCCGGCGGCGATGTTTCCCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTgg                                                                                                                                                 <  2:704955/228‑1 (MQ=255)
                                                                                            ggCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATAAAACAGCCGATTGGCGTCCCCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCt                                                                                                                                                   <  1:845315/228‑1 (MQ=255)
                                                                                                acacCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGCGGCTGAAGCCCGCCAGTCAGACGCCGGTCTCTGCGCTGGCCCCGGCGGGGCTGGCGATCCGCGCGGGCGGTCCGGCCGGGGGGTTTaaaggg                                                                                                                                    >  2:786069/1‑235 (MQ=255)
                                                                                                      ttCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAACGCCGGTCCGGCGCTGGCAGCAGGCTGCCCCATGGGGCTGAAGCCCGCCAGCCAGACGCCGTTCTCTCCGCTGGCGCTGGCGGGGCTGGCGTTCCGCGCGGGGCTGCCGGCGGGGGTGTTTTActggacctta                                                                                                                              >  2:845315/1‑230 (MQ=255)
                                                                                                              catcaGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGGATTTAACGGGGTCACCGGTTcggc                                                                                                                      >  1:170129/1‑239 (MQ=255)
                                                                                                              catcaGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAACCCGCCAGGCAGACGCCGTTCTCTGCGCTGGCCCCGGCGGGACTTGCGCCCCGCCCGGGCGTTCCGGCTGGGGGGTTTAAAGTGGGCCCCGGTGcgcg                                                                                                                      >  2:614002/1‑237 (MQ=255)
                                                                                                                          tAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCCGGCGCTGGCGGAGCTGGCGATCCGCCCGGGCGTTCCGGCCGGGGTGTTTAACGGGGTCACCGGTTCggcggggggggggg                                                                                                           >  2:10284/1‑231 (MQ=255)
                                                                                                                                    aTTGTTAGCAAGCGGCCGTTTGGGCTCACCGCGGCTATCCCGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACt                                                                                                 <  1:153485/239‑1 (MQ=255)
                                                                                                                                     tgtttATCAAGCAGCCGATTGGCGTCCCCGCGGCTATCACGCCGTGGAACTTCCCGGGGGCGATGATTCCCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACt                                                                                                 <  1:594012/235‑1 (MQ=255)
                                                                                                                                         taccAAGCAGCCGATTGGCGTCCCCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGAcc                                                                                             <  1:158479/235‑1 (MQ=255)
                                                                                                                                          agcAAGCCGCCGATGGGCGTCCCCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTCCCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAg                                                                                           <  1:786069/237‑1 (MQ=255)
                                                                                                                                                     gATTGGCGTCACCGGGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCCCCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCAc                                                                                                                               <  1:798264/192‑1 (MQ=255)
                                                                                                                                                                            aCGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGGGCGGTCGGGAACGAACCAACCCGGAACCCCCCGGTGGGCGAACCGGCGGTTACCgg                                                         >  1:553557/1‑239 (MQ=255)
                                                                                                                                                                             gccgggggACCTCCCCGGCGGAGATGATGACCCGCAAACGCGGCCCGGCGCTGGCAGCAGGCTGCCCCTTGGTGCTGACGCCCGCCAGTCAGACGCCGTTCTCTGCGCGGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGt                                                        <  2:619638/233‑1 (MQ=255)
                                                                                                                                                                             gccggTGGCACTTCCCGGCGGGATGGTTTGCCCGCAAAGCGGTGCCGGCGTGGGCGGCGGGTTCCACGATGGTGCTGAGGCCCGGCAGTCAGGCGGCGGTCGCTGCGGTGGCGCTGGTGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGt                                                        <  2:173593/235‑1 (MQ=255)
                                                                                                                                                                               ccGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAAACGTCGGTTACCGGtt                                                       >  1:515884/1‑238 (MQ=255)
                                                                                                                                                                                          cccGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGGAACGAACTGCCCAGTAAACCGCTGGGGCGGCAAATGTCGCTTTACCGTTGGACCGGAAATg                                           >  1:377864/1‑239 (MQ=255)
                                                                                                                                                                                            cggcgggcttgatTCCCCGAAGCGCCGGCCCGGCGCGGGCCGCGGGCGGCCCCTTGGTGCTGAAGCCCGCCAGTCAGCCGCCGTTCGCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGgc                                         <  2:583138/230‑1 (MQ=255)
                                                                                                                                                                                               cggcgATGCTTACCCGCAAAGCCGGTCCGGGCGTGGCAGCGGGCTGCCCCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGgccg                                       <  1:430903/238‑1 (MQ=255)
                                                                                                                                                                                                  tgatgatTACCCGAAAGGCCGGTCCGGGGCTGGACGCCGGCTGCACCTTGGTGCGGAACCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCa                                    <  1:325491/237‑1 (MQ=255)
                                                                                                                                                                                                  cgATGATTCCCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCCCCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCa                                    <  1:122602/238‑1 (MQ=255)
                                                                                                                                                                                                   cgatgttaccccGAAAGGCCGGTCCGGCACGGGCAGCGGGCTCCACACTGGTGCGGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCa                                    <  1:68905/229‑1 (MQ=255)
                                                                                                                                                                                                       gatTACCCGAAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTaa                               <  1:110030/239‑1 (MQ=255)
                                                                                                                                                                                                            cccGAAAGGGCGGGCCGGCGCTGGCAGGCGGCTGCCACAGGGGTCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGaa                          <  1:10284/239‑1 (MQ=255)
                                                                                                                                                                                                                cAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGGATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAGCGCACTGCCACGTAACCCGCTGGTGCGGCAACCGGCGGTTACCGGTGTGACCCAAACTGACCGCCAGTTAAAGGaagcag                      >  2:572776/1‑235 (MQ=255)
                                                                                                                                                                                                                      cGGTCCGGCGCGGGCCGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGa                <  1:295768/239‑1 (MQ=255)
                                                                                                                                                                                                                         tCCGGCGCGGGCGGCCGGCTGCCCCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGaaa              <  1:614002/238‑1 (MQ=255)
                                                                                                                                                                                                                         gccGGCGCTGGCAGCAGGCTGCCCCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGaaa              <  1:645215/237‑1 (MQ=255)
                                                                                                                                                                                                                                 gggCAGCAGGCTGCACCAGGGTGCGGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAg     <  1:572776/238‑1 (MQ=255)
                                                                                                                                                                                                                                   gctgtagGCTGCCCCATGGTGGTGAGTCCCGCCAGTCAGACCCCGTTCCCGGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCATTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGa    <  1:77340/233‑1 (MQ=255)
                                                                                                                                                                                                                                      cgagggCTGCCCCATGGGGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGaaa  <  1:756153/234‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                    
ACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGGCGACACCATTCCTGGTCATCAGGCCGATAAACGCCTGATTGTTATCAAGCAGCCGATTGGCGTCACCGCGGCTATCACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAA  >  NC_000913/2791563‑2792025

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: