Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,422,841 C→A 81.6% D293Y (GAT→TAT)  purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,422,8410CA81.6% 745.7 / 193.9 333D293Y (GAT→TAT) purFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/35);  new base A (120/152);  total (146/187)
Fisher's exact test for biased strand distribution p-value = 8.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 334 total aligned reads displayed.

CAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCACAT  >  CP009273/2422692‑2422989
                                                                                                                                                     |                                                                                                                                                    
  gctgctGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATAtt                                                                                                                                                    <  2:1794516/150‑1 (MQ=255)
   ctgctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTc                                                                                                                                                   <  2:4585010/150‑1 (MQ=255)
    tgctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTcc                                                                                                                                                  <  1:2495869/150‑1 (MQ=255)
    tgctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTGAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTcc                                                                                                                                                  <  2:3723408/150‑1 (MQ=255)
     gctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTccc                                                                                                                                                 <  1:5447891/150‑1 (MQ=255)
       tgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCAt                                                                                                                                               >  1:3673603/1‑150 (MQ=255)
        ggCCCGGCATGATAAAGGTGCGGCCAACATAGAGGTTTTTAACGAAGCCATGGCGGTACGGATTCCCCACAATACGCGCAATTTCCAGCGCGATATCCCACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCaat                                                                                                                                              >  2:5013120/1‑148 (MQ=255)
         gCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTc                                                                                                                                             <  1:3461669/150‑1 (MQ=255)
         gCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTc                                                                                                                                             <  2:4866396/150‑1 (MQ=255)
         gCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTc                                                                                                                                             >  1:3611076/1‑150 (MQ=255)
             ggCATGATAAAGTTGCGGCCAACATAGCGGATTTTAACGAAGCCCTGGCGGTACGGTATGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCgg                                                                                                                                         >  2:2069484/1‑150 (MQ=255)
             ggCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCgg                                                                                                                                         >  2:2516768/1‑150 (MQ=255)
             ggCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCgg                                                                                                                                         <  1:5642027/150‑1 (MQ=255)
                aTGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCa                                                                                                                                      >  2:4925527/1‑150 (MQ=255)
                    tAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAAttt                                                                                                                                  <  2:4090097/150‑1 (MQ=255)
                      aaGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTtc                                                                                                                                >  2:5555454/1‑150 (MQ=255)
                       aGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtct                                                                                                                               >  1:2813632/1‑150 (MQ=255)
                        ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtctc                                                                                                                              <  1:1335767/150‑1 (MQ=255)
                         gTGCGGCCAACATAGCGGTTTTTAACGAAGCCCGGGCGGTACGGTTTGCCCAGAATCCGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCg                                                                                                                             <  1:4366337/150‑1 (MQ=255)
                             ggCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCa                                                                                                                         <  2:1429371/150‑1 (MQ=255)
                              gCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAg                                                                                                                        <  1:2233361/150‑1 (MQ=255)
                                 aaCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGttt                                                                                                                     >  1:5214420/1‑150 (MQ=255)
                                  aCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTc                                                                                                                    >  1:613061/1‑150 (MQ=255)
                                  aCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTc                                                                                                                    <  1:3506860/150‑1 (MQ=255)
                                  aCATAGCGGTTGTGAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTc                                                                                                                    <  2:4024714/150‑1 (MQ=255)
                                    aTAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGt                                                                                                                  >  2:2911647/1‑150 (MQ=255)
                                    aTAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGt                                                                                                                  >  1:1063994/1‑150 (MQ=255)
                                     tAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTg                                                                                                                 >  1:5379653/1‑150 (MQ=255)
                                        cGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGccc                                                                                                              >  2:2319924/1‑150 (MQ=255)
                                         ggTTTTGAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGCTAGTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGGGCCCa                                                                                                             <  2:2276207/150‑1 (MQ=255)
                                             tttAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATAtt                                                                                                         >  1:735956/1‑150 (MQ=255)
                                               tAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTca                                                                                                       <  1:956034/150‑1 (MQ=255)
                                                  cGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCacac                                                                                                    <  1:3180373/150‑1 (MQ=255)
                                                   gAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcg                                                                                                   >  1:501055/1‑150 (MQ=255)
                                                    aaGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcgc                                                                                                  <  2:5228188/150‑1 (MQ=255)
                                                     aGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCATATATTCCCATTCGCGGTCAATTTTCTCGCACAGTTACGTGCCCATATTCACAcgct                                                                                                 >  2:3480425/1‑149 (MQ=255)
                                                      gCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcgcgc                                                                                                <  2:3140056/150‑1 (MQ=255)
                                                       cccTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcgcgcg                                                                                               >  2:2187513/1‑150 (MQ=255)
                                                           ggCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGt                                                                                           <  1:3598678/150‑1 (MQ=255)
                                                            gCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTa                                                                                          <  1:2985167/150‑1 (MQ=255)
                                                             cGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTaa                                                                                         <  2:1642339/150‑1 (MQ=255)
                                                                tACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACg                                                                                      <  2:778934/150‑1 (MQ=255)
                                                                 aCGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACgg                                                                                     >  2:2028317/1‑150 (MQ=255)
                                                                   ggTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaa                                                                                   <  2:1916719/150‑1 (MQ=255)
                                                                   ggTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaa                                                                                   >  1:3536987/1‑150 (MQ=255)
                                                                    gTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaaa                                                                                  >  1:1204884/1‑150 (MQ=255)
                                                                    gTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaaa                                                                                  >  1:5537145/1‑150 (MQ=255)
                                                                       tGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAAttt                                                                               >  2:5648368/1‑150 (MQ=255)
                                                                         cccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTg                                                                             <  1:2730347/150‑1 (MQ=255)
                                                                          ccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGt                                                                            >  1:212625/1‑150 (MQ=255)
                                                                             gAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGa                                                                         >  1:1806376/1‑150 (MQ=255)
                                                                              aaTACGAGCATTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGAt                                                                        <  1:3879102/150‑1 (MQ=255)
                                                                               aTACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATa                                                                       <  1:3328018/150‑1 (MQ=255)
                                                                                    aGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGa                                                                  >  2:5174526/1‑150 (MQ=255)
                                                                                     gCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCGGTAAACGGAAATTTTGTCGATAAACGAg                                                                 <  2:844198/150‑1 (MQ=255)
                                                                                      cAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGt                                                                <  1:4611659/150‑1 (MQ=255)
                                                                                       aaTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTc                                                               <  2:1226214/150‑1 (MQ=255)
                                                                                        aTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGTAAATTTTGTCGATAAACGAGTca                                                              >  2:4528392/1‑149 (MQ=255)
                                                                                             cAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcgggcg                                                         >  2:4345345/1‑150 (MQ=255)
                                                                                              aGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcgggcgg                                                        >  2:5126788/1‑150 (MQ=255)
                                                                                                 gcgATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCa                                                     <  1:1180408/150‑1 (MQ=255)
                                                                                                 gcgATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCa                                                     <  1:2096190/150‑1 (MQ=255)
                                                                                                   gATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCaaa                                                   <  2:5407473/150‑1 (MQ=255)
                                                                                                    atatCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAg                                                  <  1:2797175/150‑1 (MQ=255)
                                                                                                      atCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGta                                                <  1:1456298/150‑1 (MQ=255)
                                                                                                      atCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGta                                                >  2:4069481/1‑150 (MQ=255)
                                                                                                       tCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGACATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGtat                                               >  1:4064063/1‑150 (MQ=255)
                                                                                                         acacGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatac                                             <  1:2541295/150‑1 (MQ=255)
                                                                                                         acacGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatac                                             >  2:2753990/1‑150 (MQ=255)
                                                                                                          cacGAGGTTTCTGGGATCGAGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTataca                                            >  2:1170411/1‑150 (MQ=255)
                                                                                                            cGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatacata                                          >  1:4329793/1‑150 (MQ=255)
                                                                                                             gAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatacatac                                         >  2:5587716/1‑150 (MQ=255)
                                                                                                               ggTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTc                                       <  1:1622657/150‑1 (MQ=255)
                                                                                                                gTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCa                                      <  2:4742593/150‑1 (MQ=255)
                                                                                                                      gggatcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAgg                                >  1:3823006/1‑150 (MQ=255)
                                                                                                                       ggatcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGc                               <  2:2583314/150‑1 (MQ=255)
                                                                                                                        gatcgggatcACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCa                              >  1:4716499/1‑150 (MQ=255)
                                                                                                                        gatcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCa                              >  1:3398970/1‑150 (MQ=255)
                                                                                                                         atcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCAc                             >  1:863261/1‑150 (MQ=255)
                                                                                                                            gggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGa                          >  1:1515187/1‑150 (MQ=255)
                                                                                                                             ggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATCTTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGAt                         >  2:1326998/1‑150 (MQ=255)
                                                                                                                              gatcACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGAtt                        <  1:4925917/150‑1 (MQ=255)
                                                                                                                                tcACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGc                      >  2:1864302/1‑150 (MQ=255)
                                                                                                                                tcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGc                      <  1:632937/150‑1 (MQ=255)
                                                                                                                                 caccacGTAGATATCCATATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCt                     <  2:4900639/150‑1 (MQ=255)
                                                                                                                                  accacGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTg                    >  1:1872177/1‑150 (MQ=255)
                                                                                                                                     acGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGAcc                 <  2:3820500/150‑1 (MQ=255)
                                                                                                                                     acGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGAcc                 >  1:2423829/1‑150 (MQ=255)
                                                                                                                                         cGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGAt             >  2:3286514/1‑150 (MQ=255)
                                                                                                                                          gATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGAtt            <  1:5137092/150‑1 (MQ=255)
                                                                                                                                          gATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGAtt            <  2:2820181/150‑1 (MQ=255)
                                                                                                                                            tatCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGt          >  2:3666306/1‑150 (MQ=255)
                                                                                                                                            tatCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGt          <  1:411859/150‑1 (MQ=255)
                                                                                                                                             atCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTc         >  1:4385652/1‑150 (MQ=255)
                                                                                                                                             atCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTc         <  1:5585707/150‑1 (MQ=255)
                                                                                                                                              tCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCa        <  1:4528409/150‑1 (MQ=255)
                                                                                                                                                 aGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGca     <  2:3084782/150‑1 (MQ=255)
                                                                                                                                                  gATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGcac    >  1:4521224/1‑150 (MQ=255)
                                                                                                                                                  gATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGcac    <  1:5667801/150‑1 (MQ=255)
                                                                                                                                                    tCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCACAt  >  1:3112271/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCACAT  >  CP009273/2422692‑2422989

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: