Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,355,596 G→C 6.8% A292P (GCA→CCA)  gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,355,5960GC 6.8% 220.5 / 8.6 82A292P (GCA→CCA) gltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base G (47/27);  new base C (4/2);  total (53/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

GCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTCTGCGTCCGGCGCGCTACGTCATCACCAAAGATAAGCTCA  >  NC_000913/3355372‑3355827
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         
gcTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCCCCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATCCTGGAACTGCTGCTGACAGCCGggatggaa                                                                                                                                                                                                                             >  2:216607/1‑238 (MQ=255)
         cTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCCCTTACTCCAGTTCGATGGATAAAATGCTGGAACCGCTGCTGGCAGGCggggggggaaaaaccagg                                                                                                                                                                                                                    >  1:866024/1‑224 (MQ=255)
            acacCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACCGCTGCTGGCAGTCGGGATGGATATAATCCGCGcc                                                                                                                                                                                                                 >  2:581887/1‑239 (MQ=255)
              acCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCACGCGGGGTGGATATAAACCGCGCCaa                                                                                                                                                                                                               >  1:17951/1‑238 (MQ=255)
                 gTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACTCCGCACCGTTCGCCAACGAAACCGCCCCTGACTCCAGTTCGCAGGAAAACCTCAGGGAACTGCTCGTGGAGAGCgggagggaaaagagcccgtcaaaggg                                                                                                                                                                                                            >  1:317695/1‑221 (MQ=255)
                 gTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGCTAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGc                                                                                                                                                                                                             >  2:853349/1‑238 (MQ=255)
                 gTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGCTAACATGCTGGAACTCCTGCTGGCAGACGGGATGAATTCCCTCCTCGCCCTga                                                                                                                                                                                                             >  2:152221/1‑237 (MQ=255)
                 gTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGggatggat                                                                                                                                                                                                                             >  2:545342/1‑222 (MQ=255)
                 gTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATACTGGAACTGCTGCTGCCAGCGGGGAGGGATAACAACACCCCCATGc                                                                                                                                                                                                             >  1:223955/1‑238 (MQ=255)
                 gTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCCCTGACCCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGGgggagggaaaacaacccccccaagg                                                                                                                                                                                                             >  1:712377/1‑221 (MQ=255)
                  tACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGCACTGCTGCTGCCAGGCGGGATGGGAATCATCTGGGCAATGc                                                                                                                                                                                                             >  2:704625/1‑237 (MQ=255)
                  tACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGAGGGATatcat                                                                                                                                                                                                                        >  2:310483/1‑226 (MQ=255)
                         ctggccgctggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCCCCAACCCCAGTTCGCTGGATAACATGCTGGACCTGCTGCTCGCAGGCGGGATGGATATCACCCGCCCCGTGGGGCTATTaa                                                                                                                                                                                                    >  1:162705/1‑238 (MQ=255)
                          tggccgctggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACGGCTGCTGGCAGGCGGGATGGCGTGCATCCGCGCCATGCGCCTAtaggt                                                                                                                                                                                                   >  2:62830/1‑235 (MQ=255)
                            gccgctggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGTTAACATGCTGGAACTGCTTCTGGCAGGCGGGATGGATAGCATCCGCGCCATGCGTCTATTAGTa                                                                                                                                                                                                  >  2:804402/1‑238 (MQ=255)
                            gccgctggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGCTAACATGCTGGAACTGCGGATGGCACGCGGGATTGATATCATCCGCCCCCTGCGTCTATTAGTa                                                                                                                                                                                                  >  2:993474/1‑238 (MQ=255)
                                 tggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGTCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTaccacc                                                                                                                                                                                             >  2:894241/1‑238 (MQ=255)
                                             tCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGCTGGATACCATGCTGCAACTGCTGCTGGCCGGCGGGTTGGCTACCACCCGAGCCCTTCGGCTATTAATGCCACCCGccgtggcgaaa                                                                                                                                                                                >  1:963870/1‑230 (MQ=255)
                                              ccGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGCAACCGGCTCTGACTCCAGTTCGATGGATAACAAGCTGGAACTGCTGCTGTCAGGGGGGATGGATATCTTCCGCGCCACGAGCCCAATAGAAACACCCGCACGGCAGaaca                                                                                                                                                                               >  1:284505/1‑239 (MQ=255)
                                                 ttATCCGGGGCAAAACGTGTAAAACAACCCCACCCCCGGGAACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGCCAgcagcag                                                                                                                                                                             <  2:712377/238‑7 (MQ=255)
                                                  tattggccgccatagtgcaaaaTCACCCCCACCCCCGGTTACCCGCCAATGGGCGCGGGCGCGTACTTATAAATTCCAGACCCCGCTTATCCCTGACCGGCCCGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAAcc                                                                                                                                                                            <  1:806194/220‑1 (MQ=255)
                                                    tCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGCTATCATCCGCGCCATCCGCCTATTAGTACCACCCCCCTGGCAGAACAACCCgg                                                                                                                                                                         >  2:272893/1‑239 (MQ=255)
                                                            aTAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGGTATAAACCCCGCCATGCGTCTTTTAGTAACACCCCCCGGGCAGAAAAACCCGGATATGGa                                                                                                                                                                  >  1:446917/1‑238 (MQ=255)
                                                              aaCGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGTTTGAAATCATCCGCGCCATGGGGCTTTTATTACAACCCCCCGGGCAGAACAACCCGGAAATGGAccc                                                                                                                                                               >  1:622332/1‑239 (MQ=255)
                                                              aaCGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGAccc                                                                                                                                                               >  2:83811/1‑239 (MQ=255)
                                                              aaCGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCGGGCCGAACAACCCgg                                                                                                                                                                         >  2:293655/1‑229 (MQ=255)
                                                                      tgagtaagggacaggccggaaagccGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAAc                                                                                                                                                          <  2:356762/215‑1 (MQ=255)
                                                                           ccccacccccacttcaaCCGCCAATGGGGGCGGCCGGTGATTATAAAATTCCGGACCCCGCTTGTCCCTGGCCCGGACGAGGCGCGACGTTTCGCAAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGc                                                                                                                                                   <  1:806356/224‑1 (MQ=255)
                                                                             acCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCCCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCctt                                                                                                                                                >  2:267118/1‑239 (MQ=255)
                                                                                 tCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCttt                                                                                                                                            >  2:600690/1‑239 (MQ=255)
                                                                                 tCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCCCCTGGCAGAAAAACCCGGCTATGGCCCCGGAAATGCGGGCCTTCttt                                                                                                                                            >  1:958463/1‑239 (MQ=255)
                                                                                 cctcgcgtttcactgCCAATGGGCGTGTGCGCGTATTTATAAATTCCAGACACCGCTTATCCCGGACCTGCACGAGGCCGCTCCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCcttct                                                                                                                                              <  1:316767/225‑1 (MQ=255)
                                                                                  cccccgggcaCCCCCAATGGGCGCGTGGGCTTACTTATAAATTCCAGACCCCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCttt                                                                                                                                            <  1:894241/230‑1 (MQ=255)
                                                                                         gcgcccGCCAAGGGGCGAGTTCGCTGACTTAAAACTCCCAAACCCCGCGTAGCCCTGACCTGCACGACGGCGCACCGTTCCGGAGCGAAACCGGGTCGGGCTCCAGTTCGGTGGAGAAAATGCTGGAAGTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATAGGGGCCCGGAACTGCGTGCCTTCTTTGGCttt                                                                                                                                      <  1:744066/235‑1 (MQ=255)
                                                                                          aacccGCCAATGGGCGCCTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACttt                                                                                                                                      <  2:709335/235‑1 (MQ=255)
                                                                                                caatggggCGCTGCGCGTACTTATAAATTCCAGCCACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCAt                                                                                                                              <  2:962657/234‑1 (MQ=255)
                                                                                                 cAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCAt                                                                                                                              <  2:17951/239‑1 (MQ=255)
                                                                                                        cgcgtgcgcgtACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATTCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGGGCCTTATTTGACTTTAACTCCATGCCTATg                                                                                                                       >  2:715882/1‑239 (MQ=255)
                                                                                                         gcgtgcgcgtACTTATAAATTCCAGACCCCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATg                                                                                                                       <  2:540277/238‑1 (MQ=255)
                                                                                                              cgcgtACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGaac                                                                                                                                                                                <  2:874469/176‑1 (MQ=255)
                                                                                                              cgcgtACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGaaa                                                                                                                                                                                >  1:874469/1‑175 (MQ=255)
                                                                                                                cgtACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAAATGCGTGCCTTCTTTGACTTTAACTCCATGCCTATGGGGCCGTg                                                                                                               >  2:305909/1‑239 (MQ=255)
                                                                                                                     ttATAAATTCCAGGCACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTTACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTaa                                                                                                                                    <  2:345770/213‑1 (MQ=255)
                                                                                                                       ataAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCCGGCAGAACAACCCGGATATGGACCCGGAAATGCGTGCCTTCTTTGGCTTTAACTCCCTCCATTTGGAGCCGgggggaggc                                                                                                        >  1:945634/1‑234 (MQ=255)
                                                                                                                             ttCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCcggc                                                                                                  <  2:284505/239‑1 (MQ=255)
                                                                                                                                cAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCCAATGGAGCCGTGGGATGGCCcggcgg                                                                                                >  2:1006961/1‑238 (MQ=255)
                                                                                                                                cAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAAACCGGATACGGAACCGGACCATCCGGCCCTCCTTGAACTTACATCCAGACAAATGGAAACGGGGGATTGCCCgggggg                                                                                               >  1:919673/1‑239 (MQ=255)
                                                                                                                                  gACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCa                                                                                             <  2:945634/239‑1 (MQ=255)
                                                                                                                                       cGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATTGAGCCGTGGGATGGCCCGGCGGGGCTCGTg                                                                                        >  2:520096/1‑239 (MQ=255)
                                                                                                                                             tCCCTGGCCTGCACGACCCCCCACCTTTCGTCCACTAACCCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGGTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTcc                                                                                  <  1:520096/239‑1 (MQ=255)
                                                                                                                                                   aCCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGCCGGt                                                                            >  2:763236/1‑239 (MQ=255)
                                                                                                                                                   aCCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGGGCCTTCTTTTCCTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCgggaaaccggagtgccggcggg                                                                            >  1:289780/1‑220 (MQ=255)
                                                                                                                                                   aCCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCCCCTGGCAGAACAACCCGGATATGGACCCGGAAATTCGTGCCTTCTTTGCCTTTAACTCCCGTCATTTGGAGCCCGGGGGTGCCCCGGGGGGCGTCgggtagccccggggg                                                                            >  1:857852/1‑225 (MQ=255)
                                                                                                                                                        cACGACGCCGCCCCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGc                                                                                                                                                   <  1:981899/164‑1 (MQ=255)
                                                                                                                                                        cACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGc                                                                                                                                                   >  2:981899/1‑164 (MQ=255)
                                                                                                                                                            acgCCGCACCGTTCGGCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTgcc                                                                   <  2:137241/239‑1 (MQ=255)
                                                                                                                                                                cgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGCTGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGGTATCATCCGCGCCATGCGTCTCTTAGTCCCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTCCTTTTACTTTTAACCCGTGCAAAGGGAGCCGTGGGTGGGCCCGGGGGGCCTCCTGAGGTCCGCCGAGCGTTTTgcccccc                                                               >  1:444876/1‑238 (MQ=255)
                                                                                                                                                                cgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCg                                                                                                                                                                                                            <  2:314759/99‑1 (MQ=255)
                                                                                                                                                                cgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCg                                                                                                                                                                                                            >  1:314759/1‑99 (MQ=255)
                                                                                                                                                                  cACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGGGATGTCCGCCGGGCGTTTTGCCGCCTgg                                                             >  2:247278/1‑238 (MQ=255)
                                                                                                                                                                       ttCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCGGAACAACCCGGCAATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCAATTGGAGCCCTGGGATGGCCACGGTGGGATCGTGATGTCCCACGGTGGTTTTTCCCCCTGTaaac                                                         >  1:347465/1‑236 (MQ=255)
                                                                                                                                                                       ttCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGGGATGTCCGACGGTCGTTTTGCCGCCTGTAAcc                                                         >  2:356871/1‑238 (MQ=255)
                                                                                                                                                                        tCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTc                                                       <  2:958463/239‑1 (MQ=255)
                                                                                                                                                                        tCGTCAACGAAACCGGCTCTGACCCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTc                                                       <  2:866024/239‑1 (MQ=255)
                                                                                                                                                                          gTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGGGATGGCCGACGGTCGTTTTTCCGCCTGGAAACTCg                                                      >  2:37031/1‑238 (MQ=255)
                                                                                                                                                                                 gaaCCGGCTCGTGCCCCAGGTCGAGGGAAAACCTGCTGGAAATGCGGCGGGCGGGCGGGGTGGATATCACCCGCGCCAGGCGTCTATTAGGCCCCCCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAAc                                              <  1:658956/238‑1 (MQ=255)
                                                                                                                                                                                                  aGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTTCCGCCTGTAACCTGCACCGTAACGGTCTGCGGCCGgcgc                              >  2:48527/1‑238 (MQ=255)
                                                                                                                                                                                                           cgggcccatatgcgggttgttctgCCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAATACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTCTGCGTCCGGCGCGCTACGtcat                    <  2:279110/218‑1 (MQ=255)
                                                                                                                                                                                                            aatcaaTGCTGGACCTGTTGCTGGCAGGGGGTTTGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTCTGCGTCCGGCGCGCTACGtcat                    <  1:356871/233‑1 (MQ=255)
                                                                                                                                                                                                              aaCATGGGGGATCTTCTGCTGGCAGGGGGGATGGATATCATCCGCGCCAGGCGTCTATTAGTACCCCCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTCCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTCTGCGTCCGGCGCGCTACGTCATCAc                <  1:272893/239‑1 (MQ=255)
                                                                                                                                                                                                               aCATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTCTGCGTCCGGCGCGCTACGTCATCAcc               <  2:446917/239‑1 (MQ=255)
                                                                                                                                                                                                                   gCTGGACCTGCTGCGGGCAGGGGGGATGGATATCATCCGCGCCATGCGTTTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTCTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTCTGCGTCCGGCGCGCTACGTCATCACCAAAg           <  2:162705/239‑1 (MQ=255)
                                                                                                                                                                                                                   gCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTCTGCGTCCGGCGCGCTCCGTCATCACCAAAg           >  2:18199/1‑239 (MQ=255)
                                                                                                                                                                                                                     tGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAAAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGGGATGTCCGACGGGCGTTTTTCCGCCCGTAACCCGCACCGGAAAGGGCCGCGGCCGgcggggaaagtcacaacacaaaaa          >  1:418119/1‑218 (MQ=255)
                                                                                                                                                                                                                      ggctctggttctgGCAGGGGGGTTGGATATCACCCGCGCCATGCGTTTATTAGTCCCCCCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTCTGCGTCCGGCGCGCTACGTCATCACCAAAGATa        <  1:1006961/229‑1 (MQ=255)
                                                                                                                                                                                                                         ggttgttctgCCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCCGCGGGCCTCGTGTAGTCCGACGGGCGTTTTTCCGCCGGTAACCCCGCCCCAAACGGGGTGGGGCCGGGGGGCTACGGCGTCCACAAAAAAAAGc     >  1:279110/8‑239 (MQ=255)
                                                                                                                                                                                                                          ctgctgctgGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTCTGCGTCCGGCGCGCTACGGCATCACCAAAGATAAGc     >  2:21282/1‑238 (MQ=255)
                                                                                                                                                                                                                            gctgctgGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCCCCCGGAACCCCCACCCGAACGGTCTGGTTCCGGCGcgccacgccacaacccaaaaaaaacaca   >  1:20852/1‑221 (MQ=255)
                                                                                                                                                                                                                             gttctgCCAGGCGGGATGGATATCATCCGCGTCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCAt                                                                                                                              <  1:593517/113‑1 (MQ=255)
                                                                                                                                                                                                                             gttctgCCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCAt                                                                                                                              >  2:593517/4‑116 (MQ=255)
                                                                                                                                                                                                                             ctgctgGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTTCCGCCTGTAACCTCGACCGTAACGGACTGGGTCCGGGGCGCTACGCATCCCCCAAAGCTAAGCtcc  >  2:756510/1‑237 (MQ=255)
                                                                                                                                                                                                                             ctgctgGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTCTGCGTCCGGCGCGCTACGTACTCCCCCAAGATAAGCcca  >  2:546435/1‑238 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         
GCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTCTGCGTCCGGCGCGCTACGTCATCACCAAAGATAAGCTCA  >  NC_000913/3355372‑3355827

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