Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,067,329 A→G 18.3% intergenic (‑279/+6) epd ← / ← yggC D‑erythrose 4‑phosphate dehydrogenase/putative PanK family P‑loop kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,067,3290AG18.3% 580.9 / 119.8 237intergenic (‑279/+6)epd/yggCD‑erythrose 4‑phosphate dehydrogenase/putative PanK family P‑loop kinase
Reads supporting (aligned to +/- strand):  ref base A (88/105);  new base G (23/20);  total (111/126)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 247 total aligned reads displayed.

GCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAATCCTCCCG  >  CP009273/3067180‑3067477
                                                                                                                                                     |                                                                                                                                                    
gCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATa                                                                                                                                                      >  2:4496400/1‑150 (MQ=255)
 cTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATaa                                                                                                                                                     >  2:3028755/1‑150 (MQ=255)
  ttCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATaaa                                                                                                                                                    >  2:2421346/1‑150 (MQ=255)
     aGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTg                                                                                                                                                 <  2:2650715/150‑1 (MQ=255)
       cTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTc                                                                                                                                               <  2:2602309/150‑1 (MQ=255)
        tAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCa                                                                                                                                              <  1:992539/150‑1 (MQ=255)
        tAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCa                                                                                                                                              <  1:4573085/150‑1 (MQ=255)
        tAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCa                                                                                                                                              <  1:2005059/150‑1 (MQ=255)
           gATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCt                                                                                                                                           <  1:1127027/150‑1 (MQ=255)
            aTAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCtt                                                                                                                                          <  1:1382024/150‑1 (MQ=255)
               aGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTtgtg                                                                                                                                       >  2:150558/1‑150 (MQ=255)
               aGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTtgtg                                                                                                                                       >  2:4339601/1‑150 (MQ=255)
               aGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTtgtg                                                                                                                                       <  1:3255190/150‑1 (MQ=255)
                gCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGa                                                                                                                                      <  1:2489518/150‑1 (MQ=255)
                    aaCGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTg                                                                                                                                  <  2:3731813/150‑1 (MQ=255)
                     aCGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGa                                                                                                                                 <  2:1194820/150‑1 (MQ=255)
                       gTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATa                                                                                                                               <  2:704280/150‑1 (MQ=255)
                        tGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATAc                                                                                                                              >  2:584266/1‑150 (MQ=255)
                           aaTAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCgg                                                                                                                           >  1:1645433/1‑150 (MQ=255)
                             tAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGcc                                                                                                                         >  2:4389580/1‑150 (MQ=255)
                              aaaaGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCct                                                                                                                        >  2:1922514/1‑150 (MQ=255)
                              aaaaGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCct                                                                                                                        <  1:3594159/150‑1 (MQ=255)
                               aaaGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCctt                                                                                                                       >  1:4789876/1‑150 (MQ=255)
                                     aTGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAg                                                                                                                 <  1:1248541/150‑1 (MQ=255)
                                        tttGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCa                                                                                                              >  2:3583140/1‑150 (MQ=255)
                                          tGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCAtc                                                                                                            <  2:3871348/150‑1 (MQ=255)
                                            tCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCAtctc                                                                                                          <  2:920044/150‑1 (MQ=255)
                                            tCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCAtctc                                                                                                          <  1:3191543/150‑1 (MQ=255)
                                            tCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCAtatc                                                                                                          <  1:3457989/150‑1 (MQ=255)
                                                   gaagaaGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATc                                                                                                   >  1:4426522/1‑150 (MQ=255)
                                                   gaagaaGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATc                                                                                                   <  2:421407/150‑1 (MQ=255)
                                                     agaagaCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAg                                                                                                 <  1:2198920/150‑1 (MQ=255)
                                                      gaagaCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGa                                                                                                >  2:4778078/1‑150 (MQ=255)
                                                       aagaCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGAt                                                                                               <  1:1279899/150‑1 (MQ=255)
                                                        agaCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGAtt                                                                                              <  2:4653931/150‑1 (MQ=255)
                                                         gaCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTg                                                                                             <  2:4568042/150‑1 (MQ=255)
                                                         gaCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTg                                                                                             <  2:3368996/150‑1 (MQ=255)
                                                          aCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTgg                                                                                            <  1:1362579/150‑1 (MQ=255)
                                                          aCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTgg                                                                                            <  1:3827560/150‑1 (MQ=255)
                                                          aCATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTgg                                                                                            >  2:2127824/1‑150 (MQ=255)
                                                              ttATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTg                                                                                        >  2:3164993/1‑150 (MQ=255)
                                                               tATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGa                                                                                       <  2:3222216/150‑1 (MQ=255)
                                                               tATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGa                                                                                       <  2:2795217/150‑1 (MQ=255)
                                                                aTCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGAc                                                                                      <  1:3912713/150‑1 (MQ=255)
                                                                aTCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGAc                                                                                      <  1:1191048/150‑1 (MQ=255)
                                                                aTCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGAc                                                                                      <  2:4070970/150‑1 (MQ=255)
                                                                  cTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACgg                                                                                    >  1:101327/1‑150 (MQ=255)
                                                                  cTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACgg                                                                                    >  2:2510032/1‑150 (MQ=255)
                                                                       tcacatcacaCTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTAtt                                                                               <  1:4339601/150‑1 (MQ=255)
                                                                        cacatcacaCTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATtc                                                                              >  1:2791166/1‑150 (MQ=255)
                                                                          catcacaCTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATtcat                                                                            >  1:4684395/1‑150 (MQ=255)
                                                                          catcacaCTTTTATCCCCTGTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATtcat                                                                            >  2:2993377/1‑150 (MQ=255)
                                                                            tcacaCTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATtcatca                                                                          <  1:3175516/150‑1 (MQ=255)
                                                                                acTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGcac                                                                      >  1:4034794/1‑150 (MQ=255)
                                                                                 cTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGcaca                                                                     <  1:2972713/150‑1 (MQ=255)
                                                                                 cTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGcaca                                                                     <  2:4426522/150‑1 (MQ=255)
                                                                                  ttttATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCacac                                                                    <  1:1180855/150‑1 (MQ=255)
                                                                                   tttATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACg                                                                   >  1:1508575/1‑150 (MQ=255)
                                                                                   tttATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAAGAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACg                                                                   <  2:1385023/150‑1 (MQ=255)
                                                                                     tATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCt                                                                 <  2:1214915/150‑1 (MQ=255)
                                                                                     tATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCt                                                                 >  2:3643582/1‑150 (MQ=255)
                                                                                      aTCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTc                                                                <  2:3381202/150‑1 (MQ=255)
                                                                                      aTCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTc                                                                <  1:2500890/150‑1 (MQ=255)
                                                                                       tCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCg                                                               <  2:838883/150‑1 (MQ=255)
                                                                                       tCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCg                                                               <  2:838877/150‑1 (MQ=255)
                                                                                          ccTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACa                                                            >  2:1774483/1‑150 (MQ=255)
                                                                                             tttGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTg                                                         <  1:2257137/150‑1 (MQ=255)
                                                                                             tttGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTg                                                         <  1:2257425/150‑1 (MQ=255)
                                                                                                   ggAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCa                                                   <  1:2063226/150‑1 (MQ=255)
                                                                                                      aGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCg                                                <  1:4000983/150‑1 (MQ=255)
                                                                                                           tATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCg                                           >  2:1453680/1‑150 (MQ=255)
                                                                                                            aTTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCgg                                          <  1:886714/150‑1 (MQ=255)
                                                                                                                cAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCAt                                      >  1:2955479/1‑150 (MQ=255)
                                                                                                                    cTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATaaaa                                  <  1:595065/150‑1 (MQ=255)
                                                                                                                     tGGCGTAATAATAACCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAAt                                 >  2:872966/1‑150 (MQ=255)
                                                                                                                      ggCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATg                                <  1:712164/150‑1 (MQ=255)
                                                                                                                      ggCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATg                                <  1:1819456/150‑1 (MQ=255)
                                                                                                                      ggCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATg                                <  1:4677487/150‑1 (MQ=255)
                                                                                                                      ggCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATg                                >  1:561742/1‑150 (MQ=255)
                                                                                                                       gCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGc                               <  1:2688016/150‑1 (MQ=255)
                                                                                                                        cGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCt                              <  1:368939/150‑1 (MQ=255)
                                                                                                                         gTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCtt                             <  1:3709303/150‑1 (MQ=255)
                                                                                                                             taataaCCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGc                         <  2:2447264/150‑1 (MQ=255)
                                                                                                                              aataaCCCTACAAGAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCt                        <  2:2234032/150‑1 (MQ=255)
                                                                                                                                 aaCCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACt                     >  2:658969/1‑150 (MQ=255)
                                                                                                                                  aCCCTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACtt                    <  2:492984/150‑1 (MQ=255)
                                                                                                                                  aCCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACtt                    <  2:3967696/150‑1 (MQ=255)
                                                                                                                                   cccTACAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTg                   >  1:2867400/1‑150 (MQ=255)
                                                                                                                                   cccTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTg                   >  1:3186261/1‑150 (MQ=255)
                                                                                                                                    ccTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGc                  <  1:3124340/150‑1 (MQ=255)
                                                                                                                                      tACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCcg                >  2:1093488/1‑150 (MQ=255)
                                                                                                                                       aCAATAACTGGAATGAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCcgc               >  1:3453021/1‑150 (MQ=255)
                                                                                                                                          aTAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGtt            <  2:1645433/150‑1 (MQ=255)
                                                                                                                                          aTAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGtt            >  1:487862/1‑150 (MQ=255)
                                                                                                                                           tAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTa           <  1:1019437/150‑1 (MQ=255)
                                                                                                                                            aaCTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTaa          >  1:2337389/1‑150 (MQ=255)
                                                                                                                                              cTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAAtc        >  2:2839423/1‑150 (MQ=255)
                                                                                                                                               tGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAAtcc       >  2:3548667/1‑150 (MQ=255)
                                                                                                                                                    tAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAATCCTCCCg  >  2:3781913/1‑150 (MQ=255)
                                                                                                                                                    tAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATATCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAATCCTCCCg  >  2:4454424/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTTCAGCTAGGATAAGCGAAACGTGGAATAAAAGGAATGTTTGTCCAGCCGAAGAAGACATTTATCTGACTCACATCACACTTTTATCCCCTTTTGTGGGAAGCTTTATTCCAGGCTGGCGTAATAATAACCCTACAATAACTGGAATAAATTGTCAGCTTGTGAAGTGATACCGGCCTTCTTCAGTCATCTCTACTATCAGATTGGCCTGACGGCTATTCATCAGCACACGCTCGACATTGGGGCCATCGGTGCGGGCATAAAATGCTTCGGCTACTTGCCGCGTTAATCCTCCCG  >  CP009273/3067180‑3067477

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: