Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,198,010 T→A 25.7% intergenic (‑550/‑65) purH ← / → rrsE fused IMP cyclohydrolase/phosphoribosylaminoimidazolecarboxamide formyltransferase/16S ribosomal RNA of rrnE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,198,0100TA25.7% 248.8 / 108.3 136intergenic (‑550/‑65)purH/rrsEfused IMP cyclohydrolase/phosphoribosylaminoimidazolecarboxamide formyltransferase/16S ribosomal RNA of rrnE operon
Reads supporting (aligned to +/- strand):  ref base T (46/55);  new base A (19/16);  total (65/71)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 139 total aligned reads displayed.

AGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGC  >  CP009273/4197861‑4198157
                                                                                                                                                     |                                                                                                                                                   
aGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTACAGCGAGTCTCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAt                                                                                                                                                     <  1:3481154/150‑1 (MQ=255)
aGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAt                                                                                                                                                     <  1:3481170/150‑1 (MQ=40)
  ggTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATAc                                                                                                                                                   >  1:1966025/1‑150 (MQ=40)
  ggTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATAc                                                                                                                                                   >  1:1966018/1‑150 (MQ=40)
  ggTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATAc                                                                                                                                                   >  1:1965959/1‑150 (MQ=40)
    ttGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCa                                                                                                                                                 >  1:2441265/1‑150 (MQ=40)
     tGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCaa                                                                                                                                                <  1:2018182/150‑1 (MQ=255)
     tGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCaa                                                                                                                                                >  2:2841232/1‑150 (MQ=255)
     tGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAAATTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACTTCGCAAGACGAAAAATGAAAACaaa                                                                                                                                                >  2:2841243/1‑150 (MQ=255)
      gACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAg                                                                                                                                               >  2:1251851/1‑150 (MQ=40)
       aCTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGt                                                                                                                                              <  1:2779014/150‑1 (MQ=40)
       aCTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGt                                                                                                                                              >  1:1883719/1‑150 (MQ=255)
        ctctGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtc                                                                                                                                             <  1:2508680/150‑1 (MQ=40)
        ctctGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtc                                                                                                                                             <  1:1437443/150‑1 (MQ=40)
         tctGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtct                                                                                                                                            <  2:342879/150‑1 (MQ=40)
         tctGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtct                                                                                                                                            <  1:4095437/150‑1 (MQ=40)
         tctGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtct                                                                                                                                            >  1:604426/1‑150 (MQ=40)
          ctGAAAGATGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtctc                                                                                                                                           >  1:2508146/1‑150 (MQ=255)
          ctGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtctc                                                                                                                                           <  2:1073304/150‑1 (MQ=40)
             aaaGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAag                                                                                                                                        >  2:1929939/1‑150 (MQ=40)
             aaaGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAag                                                                                                                                        >  2:2125817/1‑150 (MQ=40)
              aaGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCTCTAAAGCGCGTCGCAACTGCTCTTAAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACTAAAAATTAATACAAAGTCTCAaga                                                                                                                                       >  1:4110157/1‑150 (MQ=14)
              aaGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAaga                                                                                                                                       <  1:911220/150‑1 (MQ=255)
              aaGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAaga                                                                                                                                       <  1:911219/150‑1 (MQ=255)
               agagGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAagag                                                                                                                                      >  2:3399758/1‑150 (MQ=40)
               agagGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAagag                                                                                                                                      >  2:2061405/1‑150 (MQ=40)
               agagGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAagag                                                                                                                                      >  2:3399521/1‑150 (MQ=40)
                 agGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTg                                                                                                                                    <  2:3659011/150‑1 (MQ=40)
                 agGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTg                                                                                                                                    <  2:3658698/150‑1 (MQ=40)
                  ggAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGa                                                                                                                                   >  1:3966124/1‑150 (MQ=40)
                  ggAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGa                                                                                                                                   >  1:2312893/1‑150 (MQ=40)
                   gAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGaa                                                                                                                                  >  2:2550677/1‑150 (MQ=255)
                    aaaGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAac                                                                                                                                 >  1:2990278/1‑150 (MQ=255)
                     aaGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAaca                                                                                                                                <  2:1585191/150‑1 (MQ=40)
                     aaGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAaca                                                                                                                                >  1:1195081/1‑150 (MQ=40)
                     aaGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGACGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAaca                                                                                                                                >  1:1195117/1‑150 (MQ=255)
                      aGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAacac                                                                                                                               <  1:4094753/150‑1 (MQ=40)
                       gCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACg                                                                                                                              <  2:1464299/150‑1 (MQ=40)
                       gCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACg                                                                                                                              <  1:316573/150‑1 (MQ=40)
                       gCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACg                                                                                                                              >  2:255299/1‑150 (MQ=255)
                       gCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACg                                                                                                                              <  1:241575/150‑1 (MQ=255)
                        cGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGt                                                                                                                             >  2:2572683/1‑150 (MQ=40)
                        cGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGt                                                                                                                             <  2:428795/150‑1 (MQ=255)
                        cGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGt                                                                                                                             <  2:428838/150‑1 (MQ=255)
                        cGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGt                                                                                                                             <  2:428704/150‑1 (MQ=255)
                        cGTAATATACGCAACCTCGCGACAGTGAGCTAAATCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAAGCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGt                                                                                                                             >  2:489398/1‑150 (MQ=255)
                         gTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTa                                                                                                                            <  2:724662/150‑1 (MQ=40)
                          tAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTaa                                                                                                                           <  1:1697120/150‑1 (MQ=40)
                           aaTATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCTAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAAt                                                                                                                          >  2:2832895/1‑150 (MQ=255)
                           aaTATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAAt                                                                                                                          <  2:549597/150‑1 (MQ=40)
                           aaTATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAAt                                                                                                                          >  2:3434376/1‑150 (MQ=40)
                           aaTATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAAt                                                                                                                          <  2:212977/150‑1 (MQ=40)
                           aaTATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAAt                                                                                                                          <  2:3902374/150‑1 (MQ=40)
                           aaTATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGTAAt                                                                                                                          >  1:1776620/1‑150 (MQ=255)
                           aaTATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAAAAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAAt                                                                                                                          <  1:1034755/150‑1 (MQ=255)
                             tataCGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTc                                                                                                                        <  2:1608218/150‑1 (MQ=40)
                             tataCGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGTAATTc                                                                                                                        <  2:3181710/150‑1 (MQ=255)
                             tataCGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGTAATTc                                                                                                                        <  2:3181660/150‑1 (MQ=255)
                                aCGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCAtt                                                                                                                     <  1:1136964/150‑1 (MQ=40)
                                 cGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa                                                                                                                    <  1:1442219/150‑1 (MQ=40)
                                  gCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTAc                                                                                                                   <  1:3944092/150‑1 (MQ=40)
                                   ccACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg                                                                                                                  <  1:2730916/150‑1 (MQ=40)
                                   ccACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGTAATTCATTACg                                                                                                                  >  2:121219/1‑150 (MQ=255)
                                    cACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGa                                                                                                                 >  1:3747020/1‑150 (MQ=40)
                                     aCCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGaa                                                                                                                <  2:1842021/150‑1 (MQ=40)
                                     aCCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGaa                                                                                                                <  1:1632220/150‑1 (MQ=40)
                                     aCCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGaa                                                                                                                <  2:3613582/150‑1 (MQ=40)
                                      ccTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAg                                                                                                               <  2:1703674/150‑1 (MQ=40)
                                       cTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt                                                                                                              <  1:3917847/150‑1 (MQ=40)
                                       cTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt                                                                                                              <  1:3917838/150‑1 (MQ=40)
                                       cTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt                                                                                                              >  2:1356020/1‑150 (MQ=255)
                                       cTCGCGACAGTGAGATAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACACTATGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGAAGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt                                                                                                              <  1:2046046/150‑1 (MQ=11)
                                         cgcgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGttt                                                                                                            >  1:287658/1‑150 (MQ=40)
                                         cgcgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGttt                                                                                                            <  1:548465/150‑1 (MQ=40)
                                         cgcgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGttt                                                                                                            <  1:1084132/150‑1 (MQ=40)
                                         cgcgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGttt                                                                                                            >  2:4061419/1‑150 (MQ=40)
                                         cgcgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGttt                                                                                                            <  1:2296897/150‑1 (MQ=40)
                                         cgcgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGttt                                                                                                            >  2:1247867/1‑150 (MQ=40)
                                          gcgACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTa                                                                                                           <  1:1425126/150‑1 (MQ=40)
                                            gACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAAt                                                                                                         >  1:1523781/1‑150 (MQ=40)
                                            gACAGTGCGCTAAAGCGCGTCGCAACTGCTCATTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAAt                                                                                                         <  2:1557943/150‑1 (MQ=255)
                                             aCAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAAtt                                                                                                        <  1:2939708/150‑1 (MQ=255)
                                               aGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCt                                                                                                      >  1:1071651/1‑150 (MQ=39)
                                               aGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCt                                                                                                      >  1:1464217/1‑150 (MQ=39)
                                               aGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCt                                                                                                      >  1:2085050/1‑150 (MQ=255)
                                                 tGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCttt                                                                                                    >  1:1913662/1‑150 (MQ=38)
                                                   cgcTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGa                                                                                                  >  1:1787386/1‑150 (MQ=35)
                                                     cTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGc                                                                                                >  2:1366566/1‑150 (MQ=21)
                                                     cTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAAACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGc                                                                                                >  1:505904/1‑150 (MQ=21)
                                                      tAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCCGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCa                                                                                               <  2:1011127/150‑1 (MQ=18)
                                                      tAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCa                                                                                               >  2:1913042/1‑150 (MQ=33)
                                                       aaaGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCAt                                                                                              <  1:3932207/150‑1 (MQ=32)
                                                         aGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCa                                                                                            >  1:1589676/1‑150 (MQ=31)
                                                         aGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCa                                                                                            <  1:3542183/150‑1 (MQ=31)
                                                         aGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCa                                                                                            <  1:3542207/150‑1 (MQ=31)
                                                         aGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCa                                                                                            <  1:3542221/150‑1 (MQ=31)
                                                                                                                                                aTGAAAACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGgaaga     <  1:2837458/150‑1 (MQ=11)
                                                                                                                                                 tGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGgaagaa    >  2:2516960/1‑150 (MQ=31)
                                                                                                                                                 tGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAAAGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGgaagaa    >  2:2517125/1‑150 (MQ=14)
                                                                                                                                                   aaTACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgc  >  1:3701836/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                   
AGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGC  >  CP009273/4197861‑4198157

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: