Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 264,231 A→C 6.1% T328P (ACT→CCT) ‡ gabD‑I → succinate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947264,2310AC 6.1% 431.7 / 14.0 132T328P (ACT→CCT) ‡gabD‑Isuccinate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (44/80);  new base C (5/3);  total (49/83)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 141 total aligned reads displayed.

TCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGTTGACGTACCGAAGACC  >  NC_002947/264082‑264380
                                                                                                                                                     |                                                                                                                                                     
tcatcTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGcaccacca                                                                                                                                                       >  1:240584/1‑150 (MQ=255)
  atcTCCAAGTACCGTAACAACGGCCAGACCTGCGTCCGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACt                                                                                                                                                     >  2:3618388/1‑150 (MQ=255)
   tctcCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTg                                                                                                                                                    <  2:2726515/150‑1 (MQ=255)
   actcCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTg                                                                                                                                                    <  1:2398132/149‑1 (MQ=255)
      ccAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGcc                                                                                                                                                 <  1:4025179/150‑1 (MQ=255)
      ccAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGcc                                                                                                                                                 >  1:4112019/1‑150 (MQ=255)
      ccAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGcc                                                                                                                                                 <  2:123899/150‑1 (MQ=255)
       cAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGccc                                                                                                                                                >  2:5053341/1‑150 (MQ=255)
          ggaCCGTACCCCCGGCCAGACCTGCGTCTGCGCCACCCGGCTCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGACGGCACCACCACTGGCCCGCt                                                                                                                                             <  2:3854256/148‑1 (MQ=255)
             ccGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGAt                                                                                                                                          >  2:1735544/1‑150 (MQ=255)
              cGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATc                                                                                                                                         <  2:1138145/150‑1 (MQ=255)
              cGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATc                                                                                                                                         >  2:2349112/1‑150 (MQ=255)
               gcaacaacGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCg                                                                                                                                        >  1:2793585/3‑150 (MQ=255)
                 aacaacGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGAt                                                                                                                                      >  2:4151154/1‑150 (MQ=255)
                  acaacGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGATg                                                                                                                                     <  1:1079911/150‑1 (MQ=255)
                   caccGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGATga                                                                                                                                    <  2:4844797/147‑2 (MQ=255)
                   caacGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATgg                                                                                                                                    >  2:1853090/1‑150 (MQ=255)
                    aacGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGATgag                                                                                                                                   <  1:1545036/150‑3 (MQ=255)
                       ggCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAg                                                                                                                                <  1:3508096/150‑1 (MQ=255)
                        gCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAgg                                                                                                                               >  1:301816/1‑150 (MQ=255)
                         ccAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGc                                                                                                                              >  1:1355983/1‑150 (MQ=255)
                         ccAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGc                                                                                                                              <  2:68800/150‑1 (MQ=255)
                         ccAAACCTGCGTCTGCGCCAACCGTCTCTATGTACAGGACGGCGTCTATGATGCCTTTGCCGAGAAGCTCTCCAAGGCGGTCGGCTCGCTCAAGACCGGCAACGGCTTCGACGAGGGCATCAATCTCGGCCCGCTGATCGATgaggcggc                                                                                                                              >  1:2053030/1‑143 (MQ=255)
                          cAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGTACGGCGTATACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGTAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCt                                                                                                                             <  1:337368/150‑1 (MQ=255)
                          cAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCt                                                                                                                             >  1:1944774/1‑150 (MQ=255)
                          cAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCt                                                                                                                             <  2:919628/150‑1 (MQ=255)
                             aCCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTc                                                                                                                          <  1:1105533/150‑1 (MQ=255)
                             aCCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTc                                                                                                                          <  2:3557138/150‑1 (MQ=255)
                                  cGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCaa                                                                                                                     <  2:1957408/150‑1 (MQ=255)
                                   gTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAg                                                                                                                    <  1:897791/150‑1 (MQ=255)
                                      tGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTc                                                                                                                 >  2:5006577/1‑150 (MQ=255)
                                      tGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTc                                                                                                                 <  1:1866792/150‑1 (MQ=255)
                                       gcgcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTcc                                                                                                                <  2:4008619/150‑1 (MQ=255)
                                       gcgcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCa                                                                                                               <  2:1234828/150‑1 (MQ=255)
                                        cgcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCa                                                                                                               <  2:1631468/150‑1 (MQ=255)
                                        cgcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCa                                                                                                               <  2:1631206/150‑1 (MQ=255)
                                        cgcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCa                                                                                                               <  2:1487734/150‑1 (MQ=255)
                                        cgcCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCa                                                                                                               <  2:1234806/150‑1 (MQ=255)
                                          ccAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAgg                                                                                                             >  2:4933152/1‑150 (MQ=255)
                                              ccGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAaca                                                                                                         >  1:4932311/1‑150 (MQ=255)
                                                    cTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGa                                                                                                   <  1:4102834/150‑1 (MQ=255)
                                                    cTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGa                                                                                                   <  1:283016/150‑1 (MQ=255)
                                                    cTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGa                                                                                                   <  1:2894941/150‑1 (MQ=255)
                                                     tACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAg                                                                                                  <  1:4917792/150‑1 (MQ=255)
                                                      aCGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAgg                                                                                                 <  1:3740733/150‑1 (MQ=255)
                                                       cGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGa                                                                                                <  1:4240025/150‑1 (MQ=255)
                                                        gTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGAc                                                                                               >  1:5037878/1‑150 (MQ=255)
                                                           cAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGcc                                                                                            <  2:4762839/150‑1 (MQ=255)
                                                              gtcGGCGTCTCCGCCGCGGTCGCCGAGAAGCTGGCCGCTGCAGGTGCCAAGCTGAAGCTCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTc                                                                                         <  2:4223861/148‑1 (MQ=255)
                                                                 ggCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGc                                                                                      <  2:3740462/150‑1 (MQ=255)
                                                                      cTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAgg                                                                                 <  1:4529660/150‑1 (MQ=255)
                                                                       tACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGgc                                                                                <  1:3566310/150‑1 (MQ=255)
                                                                         cgacgCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGgcgc                                                                              <  2:3553991/150‑1 (MQ=255)
                                                                            cgcgTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGAt                                                                                                                                      >  2:4852241/1‑91 (MQ=255)
                                                                            cgcgTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGAt                                                                                                                                      <  1:4852241/91‑1 (MQ=255)
                                                                             gcgTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCaaa                                                                          <  2:4113681/150‑1 (MQ=255)
                                                                                ttCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCACCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTg                                                                       <  2:2734884/150‑1 (MQ=255)
                                                                                ttCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCACCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTg                                                                       <  2:2275999/150‑1 (MQ=255)
                                                                                   gCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTg                                                                    >  2:2552215/1‑150 (MQ=255)
                                                                                    ccGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGt                                                                   <  1:4919324/150‑1 (MQ=255)
                                                                                     cGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTc                                                                  >  2:4994041/1‑150 (MQ=255)
                                                                                     cGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTc                                                                  >  2:2423794/1‑150 (MQ=255)
                                                                                         aaGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGt                                                              <  1:1263880/150‑1 (MQ=255)
                                                                                              ggCCGCTGCAGTTGCCAAGCTGAAGATCGGTATCGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCATAGGCGCCAAAGTGCTGTCCGGTGGCaa                                                         <  2:3483581/150‑1 (MQ=255)
                                                                                              ggCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCaa                                                         <  2:1420307/150‑1 (MQ=255)
                                                                                                 cGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCt                                                      <  1:960376/150‑1 (MQ=255)
                                                                                                  gCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTg                                                     >  2:2107925/1‑150 (MQ=255)
                                                                                                    tGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGAt                                                   <  1:4467793/150‑1 (MQ=255)
                                                                                                    tGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGAt                                                   >  1:738719/1‑150 (MQ=255)
                                                                                                    tGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGAt                                                   <  1:1661637/150‑1 (MQ=255)
                                                                                                        gTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGaa                                               <  1:1111484/150‑1 (MQ=255)
                                                                                                         ttGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAg                                              <  2:4032403/150‑1 (MQ=255)
                                                                                                         ttGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAACAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAg                                              <  2:466658/150‑1 (MQ=255)
                                                                                                          tGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAgg                                             <  2:4239681/150‑1 (MQ=255)
                                                                                                            ccAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCCAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCa                                           >  2:2574200/1‑150 (MQ=255)
                                                                                                             cAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCaa                                          >  1:672586/1‑150 (MQ=255)
                                                                                                                   gAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAActtctt                                    <  2:1294355/150‑1 (MQ=255)
                                                                                                                   gAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAActtctt                                    <  2:1294339/150‑1 (MQ=255)
                                                                                                                    aaGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACttcttc                                   <  1:94472/150‑1 (MQ=255)
                                                                                                                    aaGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACttcttc                                   <  1:3573236/150‑1 (MQ=255)
                                                                                                                    aaGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACttcttc                                   <  1:2712089/150‑1 (MQ=255)
                                                                                                                      gATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGa                                 >  1:4896966/1‑150 (MQ=255)
                                                                                                                      gATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGa                                 >  1:4896942/1‑150 (MQ=255)
                                                                                                                               cGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCAt                        <  2:1627380/150‑1 (MQ=255)
                                                                                                                               cGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCAt                        >  2:866126/1‑150 (MQ=255)
                                                                                                                                 gCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATcc                      >  2:1005701/1‑150 (MQ=255)
                                                                                                                                    tGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTgg                   <  2:4296898/150‑1 (MQ=255)
                                                                                                                                    tGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTgg                   <  2:833095/150‑1 (MQ=255)
                                                                                                                                     ggAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGt                  <  2:3948916/150‑1 (MQ=255)
                                                                                                                                      gaagaaGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGtt                 <  2:4630004/150‑1 (MQ=255)
                                                                                                                                        agaagGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGTTGa               <  1:3708538/150‑1 (MQ=255)
                                                                                                                                         gaagGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACCTCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCCTCCTGGTTGAc              >  2:4940823/1‑150 (MQ=255)
                                                                                                                                         gaagGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGTTGAc              <  2:2453393/150‑1 (MQ=255)
                                                                                                                                                 caccacTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGTTGACGTACCGaa      >  1:3164816/1‑150 (MQ=255)
                                                                                                                                                 caccacTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGTTGACGTACCGaa      >  1:608274/1‑150 (MQ=255)
                                                                                                                                                 caccacTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGTTGACGTACCGaa      >  1:3164825/1‑150 (MQ=255)
                                                                                                                                                 caccacTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGTTGACGTACCGaa      <  2:3474875/150‑1 (MQ=255)
                                                                                                                                                   ccacTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGTTGACGTACCGAAGa    <  1:1212008/150‑1 (MQ=255)
                                                                                                                                                   ccacTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGTTGACGTACCGAAGa    >  1:4357476/1‑150 (MQ=255)
                                                                                                                                                     acTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCATGGTTGACGTACCGAAGAcc  >  1:3645808/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGACGCGTTCGCCGAGAAGCTGGCCGCTGCAGTTGCCAAGCTGAAGATCGGTAACGGCCTGGAAGAAGGCACCACCACTGGCCCGCTGATCGATGGCAAGGCTGTCGCCAAGGTCCAGGAACACATCGAGGACGCCGTCAGCAAAGGCGCCAAAGTGCTGTCCGGTGGCAAGCTGATCGAAGGCAACTTCTTCGAGCCGACCATCCTGGTTGACGTACCGAAGACC  >  NC_002947/264082‑264380

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