Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,990,572 C→A 45.4% I10I (ATC→ATA yigB → 5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase; pyrimidine phosphatase; riboflavin synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,990,5720CA45.4% 218.3 / 315.3 250I10I (ATC→ATAyigB5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase; pyrimidine phosphatase; riboflavin synthesis
Reads supporting (aligned to +/- strand):  ref base C (68/68);  new base A (46/68);  total (114/136)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

GTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGC  >  CP009273/3990424‑3990721
                                                                                                                                                    |                                                                                                                                                     
gtgCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCAtct                                                                                                                                                      <  1:1836922/150‑1 (MQ=255)
    aGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGgc                                                                                                                                                  <  1:417389/150‑1 (MQ=255)
    aGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGgc                                                                                                                                                  >  2:1012615/1‑150 (MQ=255)
    aGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGgc                                                                                                                                                  >  2:2743828/1‑150 (MQ=255)
    aGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGgc                                                                                                                                                  <  2:964099/150‑1 (MQ=255)
     ggAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGgcg                                                                                                                                                 <  1:3053266/150‑1 (MQ=255)
         ctgctgGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCa                                                                                                                                             <  1:525848/150‑1 (MQ=255)
          tgctgGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCAc                                                                                                                                            >  1:2310184/1‑150 (MQ=255)
           gctgGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCAcc                                                                                                                                           >  1:2301055/1‑150 (MQ=255)
            ctgGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCt                                                                                                                                          <  2:1564483/150‑1 (MQ=255)
               ggTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTg                                                                                                                                       <  2:4102213/150‑1 (MQ=255)
                 tCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGAc                                                                                                                                     <  1:33680/150‑1 (MQ=255)
                   aTGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCt                                                                                                                                   <  1:1739311/150‑1 (MQ=255)
                     gCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTgg                                                                                                                                 <  2:2284331/150‑1 (MQ=255)
                            tctcCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATAc                                                                                                                          >  1:4119939/1‑150 (MQ=255)
                            tctcCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATAc                                                                                                                          >  1:3028534/1‑150 (MQ=255)
                              tcCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATAccc                                                                                                                        >  2:3254795/1‑150 (MQ=255)
                                  ccacGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTa                                                                                                                    <  1:1155665/150‑1 (MQ=255)
                                        aaaTCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATaa                                                                                                              <  1:3406686/150‑1 (MQ=255)
                                          aTCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAAcc                                                                                                            >  1:2722606/1‑150 (MQ=255)
                                            cTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGt                                                                                                          >  2:2230186/1‑150 (MQ=255)
                                             tataCTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTc                                                                                                         >  1:2524579/1‑150 (MQ=255)
                                              ataCTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTcc                                                                                                        >  2:709124/1‑150 (MQ=255)
                                                    aTCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGAt                                                                                                  >  1:4146373/1‑150 (MQ=255)
                                                        tGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTg                                                                                              >  2:3862222/1‑150 (MQ=255)
                                                        tGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTg                                                                                              >  1:1243696/1‑150 (MQ=255)
                                                        tGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTg                                                                                              >  2:3541117/1‑150 (MQ=255)
                                                         gATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTgc                                                                                             >  2:2804328/1‑150 (MQ=255)
                                                          aTTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTgcg                                                                                            >  2:2023954/1‑150 (MQ=255)
                                                          aTTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTgcg                                                                                            >  2:212611/1‑150 (MQ=255)
                                                           ttttCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTgcgc                                                                                           >  2:2507135/1‑150 (MQ=255)
                                                            tttCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCa                                                                                          >  1:321819/1‑150 (MQ=255)
                                                               cAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACcg                                                                                       >  1:586337/1‑150 (MQ=255)
                                                                   acCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACcgagcg                                                                                   <  2:3000237/150‑1 (MQ=255)
                                                                     cTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACcgagcgag                                                                                 >  2:577721/1‑150 (MQ=255)
                                                                     cTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACcgagcgag                                                                                 <  1:2445519/150‑1 (MQ=255)
                                                                      ttGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCgaga                                                                                <  1:2202492/150‑1 (MQ=255)
                                                                      ttGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCgaga                                                                                <  1:121858/150‑1 (MQ=255)
                                                                         ccTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGgc                                                                             <  2:1438323/150‑1 (MQ=255)
                                                                          cTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGgcg                                                                            >  2:3588929/1‑150 (MQ=255)
                                                                          cTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGgcg                                                                            >  2:1714057/1‑150 (MQ=255)
                                                                             ggTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCtt                                                                         <  2:1980596/150‑1 (MQ=255)
                                                                                gtACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTAcc                                                                      >  1:2891959/1‑150 (MQ=255)
                                                                                  aCGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCtt                                                                    >  2:4155038/1‑150 (MQ=255)
                                                                                    gATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTtg                                                                  >  1:857538/1‑150 (MQ=255)
                                                                                       gcggcgCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGc                                                               <  2:1353048/150‑1 (MQ=255)
                                                                                        cggcgCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCa                                                              <  2:2321165/150‑1 (MQ=255)
                                                                                        cggcgCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCa                                                              <  2:1162185/150‑1 (MQ=255)
                                                                                           cgcATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCaaaa                                                           >  1:3698243/1‑150 (MQ=255)
                                                                                            gcATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATACCCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAAt                                                          <  2:2449889/150‑1 (MQ=255)
                                                                                            gcATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAAt                                                          <  2:2449840/150‑1 (MQ=255)
                                                                                             cATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAAtt                                                         >  1:602756/1‑150 (MQ=255)
                                                                                                ccACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTatc                                                      <  1:1658200/150‑1 (MQ=255)
                                                                                                ccACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTatc                                                      <  1:626344/150‑1 (MQ=255)
                                                                                                  aCGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTatcat                                                    <  1:2139532/150‑1 (MQ=255)
                                                                                                   cgcgCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTatcatc                                                   <  2:3644155/150‑1 (MQ=255)
                                                                                                    gcgcCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATcc                                                  >  1:3445293/1‑150 (MQ=255)
                                                                                                     cgcCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCg                                                 <  2:3214974/150‑1 (MQ=255)
                                                                                                     cgcCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCg                                                 >  2:2791175/1‑150 (MQ=255)
                                                                                                     cgcCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCg                                                 <  1:3578051/150‑1 (MQ=255)
                                                                                                        cAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGgcg                                              <  1:3152416/150‑1 (MQ=255)
                                                                                                          aaCGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCt                                            <  1:4124517/150‑1 (MQ=255)
                                                                                                          aaCGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCt                                            <  1:485471/150‑1 (MQ=255)
                                                                                                           aCGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTg                                           >  1:2774947/1‑150 (MQ=255)
                                                                                                           aCGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTg                                           >  1:3086989/1‑150 (MQ=255)
                                                                                                             gggggAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTgcg                                         <  1:2531748/150‑1 (MQ=255)
                                                                                                             gggggAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTgcg                                         >  1:797585/1‑150 (MQ=255)
                                                                                                              ggggAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTgcgc                                        >  1:2439619/1‑150 (MQ=255)
                                                                                                              ggggAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTgcgc                                        >  1:1288336/1‑150 (MQ=255)
                                                                                                              ggggAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTgcgc                                        >  2:2693561/1‑150 (MQ=255)
                                                                                                               gggAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCa                                       <  1:1472086/150‑1 (MQ=255)
                                                                                                               gggAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCa                                       >  1:4104525/1‑150 (MQ=255)
                                                                                                               gggAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCa                                       <  2:2243255/150‑1 (MQ=255)
                                                                                                                  aaaTAATGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCt                                    >  1:3606407/1‑150 (MQ=255)
                                                                                                                     taatGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTcc                                 >  2:3792868/1‑150 (MQ=255)
                                                                                                                     taatGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTcc                                 >  2:3318946/1‑150 (MQ=255)
                                                                                                                      aatGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCa                                <  2:48289/150‑1 (MQ=255)
                                                                                                                      aatGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCa                                <  2:3425652/150‑1 (MQ=255)
                                                                                                                      aatGCGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATAACCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCa                                >  2:2143095/1‑150 (MQ=255)
                                                                                                                          cGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAAt                            <  1:3424392/150‑1 (MQ=255)
                                                                                                                          cGTTTTTACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAAt                            <  2:333299/150‑1 (MQ=255)
                                                                                                                            tttttACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGa                          <  1:1624197/150‑1 (MQ=255)
                                                                                                                            tttttACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGa                          <  1:3917680/150‑1 (MQ=255)
                                                                                                                              tttACCGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAg                        <  1:2023952/150‑1 (MQ=255)
                                                                                                                                tACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATAACCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTGTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGTAGAt                      >  1:3871791/1‑150 (MQ=255)
                                                                                                                                  ccGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCt                    >  2:1548258/1‑150 (MQ=255)
                                                                                                                                  ccGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCt                    >  2:1885120/1‑150 (MQ=255)
                                                                                                                                   cGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTg                   >  1:1032082/1‑150 (MQ=255)
                                                                                                                                   cGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTg                   <  2:2914486/150‑1 (MQ=255)
                                                                                                                                   cGGCCTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTg                   <  1:1090352/150‑1 (MQ=255)
                                                                                                                                      ccTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCaa                >  2:353817/1‑150 (MQ=255)
                                                                                                                                       cTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAAc               <  1:3284008/150‑1 (MQ=255)
                                                                                                                                       cTTTGGGGCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAAc               <  2:2880720/150‑1 (MQ=255)
                                                                                                                                        tttGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACg              <  1:3781485/150‑1 (MQ=255)
                                                                                                                                          tGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGcc            <  1:2640758/150‑1 (MQ=255)
                                                                                                                                           ggggCGCATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCt           <  2:2565557/150‑1 (MQ=255)
                                                                                                                                                 cATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCa     <  1:2014938/150‑1 (MQ=255)
                                                                                                                                                 cATATCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATCACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCa     <  2:453066/150‑1 (MQ=255)
                                                                                                                                                    ctcGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGc  >  1:1516872/1‑150 (MQ=255)
                                                                                                                                                    ctcGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGc  <  1:1858135/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGCAGGAGCTGCTGGGTCATGCCAACCTCTCCACCACGCAAATCTATACTCATCTTGATTTTCAACACCTTGCCTCGGTGTACGATGCGGCGCATCCACGCGCCAAACGGGGGAAATAATGCGTTTTTACCGGCCTTTGGGGCGCATCTCGGCGCTCACCTTTGACCTGGATGATACCCTTTACGATAACCGTCCGGTGATTTTGCGCACCGAGCGAGAGGCGCTTACCTTTGTGCAAAATTATCATCCGGCGCTGCGCAGCTTCCAGAATGAAGATCTGCAACGCCTGCGCCAGGC  >  CP009273/3990424‑3990721

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: