Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,924,898 T→A 12.7% D36V (GAT→GTT)  yebG ← DNA damage‑inducible protein regulated by LexA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,924,8980TA12.7% 811.7 / 89.3 251D36V (GAT→GTT) yebGDNA damage‑inducible protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (95/124);  new base A (19/13);  total (114/137)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

AGCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATATTT  >  CP009273/1924750‑1925046
                                                                                                                                                    |                                                                                                                                                    
 gCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATc                                                                                                                                                     <  1:4561982/149‑1 (MQ=255)
 gCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATc                                                                                                                                                     <  1:1414046/149‑1 (MQ=255)
  cACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATcc                                                                                                                                                    <  2:3672246/149‑1 (MQ=255)
       acgacCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGt                                                                                                                                               >  1:3152426/1‑149 (MQ=255)
       acgacCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGt                                                                                                                                               <  1:975693/149‑1 (MQ=255)
          acCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATc                                                                                                                                            <  1:1972817/149‑1 (MQ=255)
            cTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCaa                                                                                                                                          <  2:551964/149‑1 (MQ=255)
            cTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCaa                                                                                                                                          >  1:3801716/1‑149 (MQ=255)
               tgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCa                                                                                                                                       <  1:1491103/149‑1 (MQ=255)
               tgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCa                                                                                                                                       <  1:3567059/149‑1 (MQ=255)
                ggggACAGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCAt                                                                                                                                      >  2:2720439/1‑149 (MQ=255)
                  ggACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCAttt                                                                                                                                    <  1:3191891/149‑1 (MQ=255)
                  ggACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCAttt                                                                                                                                    >  1:604112/1‑149 (MQ=255)
                    aCGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTg                                                                                                                                  >  2:722435/1‑149 (MQ=255)
                    aCGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTg                                                                                                                                  >  2:4013020/1‑149 (MQ=255)
                       gTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCa                                                                                                                               <  1:4076540/149‑1 (MQ=255)
                          aTTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATaa                                                                                                                            <  2:169441/149‑1 (MQ=255)
                          aTTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATaa                                                                                                                            >  2:4098504/1‑149 (MQ=255)
                           tttGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAg                                                                                                                           >  2:382925/1‑149 (MQ=255)
                           tttGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAg                                                                                                                           >  2:382717/1‑149 (MQ=255)
                            ttGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGc                                                                                                                          <  2:4156638/149‑1 (MQ=255)
                               ccGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATc                                                                                                                       <  1:680255/149‑1 (MQ=255)
                               ccGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATc                                                                                                                       >  1:868426/1‑149 (MQ=255)
                                 ggTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCgg                                                                                                                     >  1:3543520/1‑149 (MQ=255)
                                  gTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGc                                                                                                                    <  2:1608464/149‑1 (MQ=255)
                                  gTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGc                                                                                                                    >  2:3104870/1‑149 (MQ=255)
                                  gTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGc                                                                                                                    >  1:2126631/1‑149 (MQ=255)
                                    tttCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCtt                                                                                                                  <  2:1874009/149‑1 (MQ=255)
                                    tttCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCtt                                                                                                                  <  1:383660/149‑1 (MQ=255)
                                     ttCAGAATGGTGCTCAGCTCATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTc                                                                                                                 >  2:4378619/1‑149 (MQ=255)
                                     ttCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTc                                                                                                                 >  2:2877328/1‑149 (MQ=255)
                                         gAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttt                                                                                                             >  1:1743139/1‑149 (MQ=255)
                                         gAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttt                                                                                                             <  1:4027565/149‑1 (MQ=255)
                                         gAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttt                                                                                                             <  2:4339171/149‑1 (MQ=255)
                                            tGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTg                                                                                                          >  1:4439455/1‑149 (MQ=255)
                                             ggTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGc                                                                                                         <  2:1027134/149‑1 (MQ=255)
                                             ggTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGc                                                                                                         <  2:3611974/149‑1 (MQ=255)
                                               tGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTg                                                                                                       <  1:2768673/149‑1 (MQ=255)
                                               tGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTg                                                                                                       <  1:426623/149‑1 (MQ=255)
                                                  tCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTa                                                                                                    >  1:2028960/1‑149 (MQ=255)
                                                    aGCACATCCTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTaaa                                                                                                  >  1:3725056/1‑149 (MQ=255)
                                                     gCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAAc                                                                                                 >  2:321055/1‑149 (MQ=255)
                                                       acaTCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGa                                                                                               >  1:1855248/1‑149 (MQ=255)
                                                         aTCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACa                                                                                             >  1:3775907/1‑149 (MQ=255)
                                                            tttttGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAttt                                                                                          <  2:3437707/149‑1 (MQ=255)
                                                            tttttGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAttt                                                                                          <  1:1559697/149‑1 (MQ=255)
                                                               ttGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTtct                                                                                       >  1:3610605/1‑149 (MQ=255)
                                                                  ttCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCtt                                                                                    <  2:2174787/149‑1 (MQ=255)
                                                                  ttCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCtt                                                                                    <  1:1787007/149‑1 (MQ=255)
                                                                       ccagccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGccc                                                                               <  2:2037614/149‑1 (MQ=255)
                                                                        cagccaTAGCGACAGGGCTTCACGTTGCTCGTCTTCCATGTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCt                                                                              <  1:2952558/149‑1 (MQ=255)
                                                                        cagccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCt                                                                              <  2:3137540/149‑1 (MQ=255)
                                                                         agccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTc                                                                             <  2:4510727/149‑1 (MQ=255)
                                                                          gccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCa                                                                            <  1:651078/149‑1 (MQ=255)
                                                                          gccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCa                                                                            <  1:2991724/149‑1 (MQ=255)
                                                                          gccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCa                                                                            <  1:2141498/149‑1 (MQ=255)
                                                                            caTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACg                                                                          >  1:1032830/1‑149 (MQ=255)
                                                                              tAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGaa                                                                        >  2:4480291/1‑149 (MQ=255)
                                                                              tAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGaa                                                                        >  2:1534473/1‑149 (MQ=255)
                                                                               aGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAAt                                                                       >  1:146157/1‑149 (MQ=255)
                                                                               aGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAAt                                                                       >  1:1733177/1‑149 (MQ=255)
                                                                                   aaaGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACt                                                                   >  1:1181900/1‑149 (MQ=255)
                                                                                   aaaGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACt                                                                   >  2:4628128/1‑149 (MQ=255)
                                                                                   aaaGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACt                                                                   >  2:2245983/1‑149 (MQ=255)
                                                                                    aaGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTa                                                                  >  2:803035/1‑149 (MQ=255)
                                                                                     aGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTAc                                                                 >  2:3357617/1‑149 (MQ=255)
                                                                                         cTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTAt                                                             <  2:938639/149‑1 (MQ=255)
                                                                                         cTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTAt                                                             >  1:572791/1‑149 (MQ=255)
                                                                                               gTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACt                                                       >  2:392888/1‑149 (MQ=255)
                                                                                                   cTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCaa                                                   <  2:3612000/149‑1 (MQ=255)
                                                                                                        cTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGcc                                              <  1:357912/149‑1 (MQ=255)
                                                                                                        cTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGcc                                              >  2:3047496/1‑149 (MQ=255)
                                                                                                          tCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCAt                                            <  2:1401419/149‑1 (MQ=255)
                                                                                                            cATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCAttt                                          <  2:366905/149‑1 (MQ=255)
                                                                                                             aTTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTg                                         >  1:2446156/1‑149 (MQ=255)
                                                                                                              tttGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGc                                        >  2:3226248/1‑149 (MQ=255)
                                                                                                              tttGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGc                                        >  2:3898020/1‑149 (MQ=255)
                                                                                                                 gCACTGGAGAATTTGTCAGTCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGc                                     <  1:2252550/149‑1 (MQ=255)
                                                                                                                 gCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGc                                     <  1:1028668/149‑1 (MQ=255)
                                                                                                                   aCTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCct                                   >  2:243799/1‑149 (MQ=255)
                                                                                                                    cTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCctc                                  <  2:3610483/149‑1 (MQ=255)
                                                                                                                     tGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCctct                                 <  1:3604431/149‑1 (MQ=255)
                                                                                                                      ggAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCtctc                                <  2:3554242/149‑1 (MQ=255)
                                                                                                                       gagaATTTGTCAGCCAGGTGTCAGGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCa                               <  1:3354740/149‑1 (MQ=255)
                                                                                                                          aaTTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCataa                            >  1:810989/1‑149 (MQ=255)
                                                                                                                           aTTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCataat                           >  1:3630768/1‑149 (MQ=255)
                                                                                                                              tGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAAc                        <  2:1957671/149‑1 (MQ=255)
                                                                                                                                   gCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGa                   >  1:1236798/1‑149 (MQ=255)
                                                                                                                                    ccAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGAt                  >  1:3496109/1‑149 (MQ=255)
                                                                                                                                     cAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGAtt                 <  1:3990316/149‑1 (MQ=255)
                                                                                                                                       ggTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGAtttt               >  2:2126935/1‑149 (MQ=255)
                                                                                                                                        gtgtCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTta              <  2:2803921/149‑1 (MQ=255)
                                                                                                                                         tgtCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTtat             <  1:911405/149‑1 (MQ=255)
                                                                                                                                          gtCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTtata            <  2:1300494/149‑1 (MQ=255)
                                                                                                                                           tCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATAc           <  1:4602141/149‑1 (MQ=255)
                                                                                                                                            cAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACa          <  2:2277494/149‑1 (MQ=255)
                                                                                                                                               gCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGta       <  1:2875346/149‑1 (MQ=255)
                                                                                                                                                 aGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGtata     >  1:3243608/1‑149 (MQ=255)
                                                                                                                                                  gATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTatat    <  2:1248729/149‑1 (MQ=255)
                                                                                                                                                   aTCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATAtt   <  1:1411366/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATATTT  >  CP009273/1924750‑1925046

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: