Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,924,898 T→A 7.9% 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,8980TA 7.9% 1085.9 / 63.2 312D36V (GAT→GTT) yebGDNA damage‑inducible protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (140/147);  new base A (7/18);  total (147/165)
Fisher's exact test for biased strand distribution p-value = 5.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 315 total aligned reads displayed.

CAGCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATATTTC  >  CP009273/1924749‑1925047
                                                                                                                                                     |                                                                                                                                                     
   cACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATccg                                                                                                                                                    <  1:3202256/150‑1 (MQ=255)
     ccGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATccgcc                                                                                                                                                  <  2:2025787/150‑1 (MQ=255)
      cgacgaCCTGTGGGGACGGTAATTTGCCGGTTTTCAGAAGGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATccgccg                                                                                                                                                 <  2:3775491/150‑1 (MQ=255)
        acgacCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTa                                                                                                                                               >  2:2357252/1‑150 (MQ=255)
          gacCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATc                                                                                                                                             <  2:1032283/150‑1 (MQ=255)
              tgtgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGc                                                                                                                                         <  1:1251501/150‑1 (MQ=255)
              tgtgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGc                                                                                                                                         <  1:4697217/150‑1 (MQ=255)
               gtgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCa                                                                                                                                        <  1:707329/150‑1 (MQ=255)
                     aCGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGt                                                                                                                                  <  1:1020552/150‑1 (MQ=255)
                      cGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTc                                                                                                                                 <  1:1057094/150‑1 (MQ=255)
                       ggTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCa                                                                                                                                >  1:4377221/1‑150 (MQ=255)
                       ggTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCa                                                                                                                                >  1:2755372/1‑150 (MQ=255)
                       ggTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCa                                                                                                                                <  2:4756892/150‑1 (MQ=255)
                         tAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATa                                                                                                                              >  2:1147750/1‑150 (MQ=255)
                           aTTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAg                                                                                                                            >  1:363929/1‑150 (MQ=255)
                           aTTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTGTGTCATAAg                                                                                                                            >  2:2919544/1‑150 (MQ=255)
                            tttGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGc                                                                                                                           <  2:380330/150‑1 (MQ=255)
                            tttGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGc                                                                                                                           <  2:424419/150‑1 (MQ=255)
                              tGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCAt                                                                                                                         <  1:2531200/150‑1 (MQ=255)
                                   gTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCt                                                                                                                    >  1:4076498/1‑150 (MQ=255)
                                   gTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCt                                                                                                                    <  1:2900929/150‑1 (MQ=255)
                                     tttCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTc                                                                                                                  <  1:3492115/150‑1 (MQ=255)
                                     tttCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGAGAAAGGGCTTCACGTTGCTCGTATTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTc                                                                                                                  <  2:3492075/150‑1 (MQ=255)
                                          gAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCtttt                                                                                                             <  1:5088932/150‑1 (MQ=255)
                                          gAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATATTGTCATAAGCATCGGCTTCCtttt                                                                                                             >  1:3025290/1‑150 (MQ=255)
                                           aaTGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttttt                                                                                                            >  1:1366982/1‑150 (MQ=255)
                                           aaTGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttttt                                                                                                            <  2:3150395/150‑1 (MQ=255)
                                           aaTGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttttt                                                                                                            <  1:4218054/150‑1 (MQ=255)
                                                 gCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGt                                                                                                      >  1:921197/1‑150 (MQ=255)
                                                 gCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGt                                                                                                      >  1:3978589/1‑150 (MQ=255)
                                                   tCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTaa                                                                                                    <  1:3600488/150‑1 (MQ=255)
                                                   tCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTaa                                                                                                    >  2:1246950/1‑150 (MQ=255)
                                                     aGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAAc                                                                                                  >  1:2957880/1‑150 (MQ=255)
                                                     aGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAAc                                                                                                  >  1:3765753/1‑150 (MQ=255)
                                                     aGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAAc                                                                                                  >  2:2700846/1‑150 (MQ=255)
                                                      gCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACg                                                                                                 >  1:4748928/1‑150 (MQ=255)
                                                       cacaTCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGa                                                                                                <  1:2426593/150‑1 (MQ=255)
                                                          aTCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAt                                                                                             <  2:4119750/150‑1 (MQ=255)
                                                           tCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAtt                                                                                            <  1:4993590/150‑1 (MQ=255)
                                                            cTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCTCTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAttt                                                                                           <  2:721728/150‑1 (MQ=255)
                                                            cTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAgtt                                                                                           >  1:3580507/1‑150 (MQ=255)
                                                            cTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAttt                                                                                           <  1:2260035/150‑1 (MQ=255)
                                                                ttGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTtctc                                                                                       <  1:2233500/150‑1 (MQ=255)
                                                                    tCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTGTTGCTGGTAAACGACATTTTCTCTTCg                                                                                   >  1:2012651/1‑150 (MQ=255)
                                                                      tGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGcc                                                                                 <  2:1310622/150‑1 (MQ=255)
                                                                          agccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCa                                                                             <  1:86456/150‑1 (MQ=255)
                                                                          agccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCa                                                                             >  1:1296638/1‑150 (MQ=255)
                                                                          agccaTAGCGAAAGGGCTTCACGTTGCTAGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCa                                                                             <  2:3622474/150‑1 (MQ=255)
                                                                            ccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACg                                                                           >  1:2856191/1‑150 (MQ=255)
                                                                            ccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACg                                                                           >  1:4110161/1‑150 (MQ=255)
                                                                             caTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGa                                                                          >  1:777179/1‑150 (MQ=255)
                                                                             caTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGa                                                                          >  1:2940246/1‑150 (MQ=255)
                                                                             caTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGa                                                                          >  2:2066406/1‑150 (MQ=255)
                                                                                    aaaGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTa                                                                   <  2:3517483/150‑1 (MQ=255)
                                                                                     aaGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTAc                                                                  >  2:4175243/1‑150 (MQ=255)
                                                                                      aGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACg                                                                 <  1:4983509/150‑1 (MQ=255)
                                                                                        ggCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAATCGACATTTTCTCTTCGCCCTCACGAATGTCTAAGTa                                                               >  2:3669025/1‑150 (MQ=255)
                                                                                         gCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTAt                                                              <  2:630718/150‑1 (MQ=255)
                                                                                          cTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTAtt                                                             >  1:3922158/1‑150 (MQ=255)
                                                                                           ttCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCATTTTGCTGGTAAACGCCCTATTCTCTTCGCCCTCACGAATGACTACGTAttt                                                            >  2:2817542/1‑150 (MQ=255)
                                                                                             cACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTaa                                                          >  2:2176144/1‑150 (MQ=255)
                                                                                             cACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTaa                                                          >  2:464113/1‑150 (MQ=255)
                                                                                               cGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACt                                                        <  2:5131941/150‑1 (MQ=255)
                                                                                                 ttGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTc                                                      <  1:3369796/150‑1 (MQ=255)
                                                                                                 ttGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTc                                                      <  1:4884405/150‑1 (MQ=255)
                                                                                                  tGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCa                                                     >  2:1548346/1‑150 (MQ=255)
                                                                                                  tGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCa                                                     >  2:3065802/1‑150 (MQ=255)
                                                                                                    cTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAAc                                                   >  2:808245/1‑150 (MQ=255)
                                                                                                      cGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCg                                                 >  2:2511145/1‑150 (MQ=255)
                                                                                                          ttCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCAt                                             <  2:4420798/150‑1 (MQ=255)
                                                                                                           tCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCAtt                                            <  2:1226989/150‑1 (MQ=255)
                                                                                                            ccATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCAttt                                           <  2:4192859/150‑1 (MQ=255)
                                                                                                             cATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTg                                          >  1:4480936/1‑150 (MQ=255)
                                                                                                              aTTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGc                                         >  2:780475/1‑150 (MQ=255)
                                                                                                               tttGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCCTTTGCa                                        >  1:1400244/1‑150 (MQ=255)
                                                                                                                ttGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAg                                       <  2:4073023/150‑1 (MQ=255)
                                                                                                                ttGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTGTGTCATAAGCATCGGCGTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCACCCGCCATTTGCAg                                       <  2:1414249/150‑1 (MQ=255)
                                                                                                                  gCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGcc                                     >  1:4565444/1‑150 (MQ=255)
                                                                                                                    aCTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCGTACGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCctc                                   >  1:1649077/1‑150 (MQ=255)
                                                                                                                    aCTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCctc                                   >  1:4653489/1‑150 (MQ=255)
                                                                                                                     cTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCTTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCATCCGCCATTTGCAGCctct                                  <  2:1271294/150‑1 (MQ=255)
                                                                                                                     cTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCctct                                  >  1:2688617/1‑150 (MQ=255)
                                                                                                                        gagaATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCat                               <  2:2623540/150‑1 (MQ=255)
                                                                                                                         agaATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCata                              >  2:4697681/1‑150 (MQ=255)
                                                                                                                           aaTTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCataat                            <  2:5015706/150‑1 (MQ=255)
                                                                                                                           aaTTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCataat                            <  2:4748928/150‑1 (MQ=255)
                                                                                                                            aTTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCataata                           >  1:2575314/1‑150 (MQ=255)
                                                                                                                                     ccAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGAtt                  >  2:3334221/1‑150 (MQ=255)
                                                                                                                                        ggTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTta               <  2:4530659/150‑1 (MQ=255)
                                                                                                                                        ggTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTta               <  1:1900627/150‑1 (MQ=255)
                                                                                                                                        ggTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTta               <  1:1906581/150‑1 (MQ=255)
                                                                                                                                         gtgtCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTtat              >  1:1737746/1‑150 (MQ=255)
                                                                                                                                            tCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACa           >  1:1348452/1‑150 (MQ=255)
                                                                                                                                               aGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGta        <  2:4005248/150‑1 (MQ=255)
                                                                                                                                               aGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGta        >  2:5225365/1‑150 (MQ=255)
                                                                                                                                                 cAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGtata      <  2:676420/150‑1 (MQ=255)
                                                                                                                                                 cAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGtata      >  2:1970368/1‑150 (MQ=255)
                                                                                                                                                 cAGATCCGCAGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGtata      <  2:1865877/150‑1 (MQ=255)
                                                                                                                                                  aGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTatat     >  1:513571/1‑150 (MQ=255)
                                                                                                                                                    aTCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATAttt   >  1:57526/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATATTTC  >  CP009273/1924749‑1925047

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: