Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,584,756 A→T 46.7% *305K (TAA→AAA)  yjjM ← yjjMN operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,584,7560AT46.7% 68.8 / inf 313*305K (TAA→AAA) yjjMyjjMN operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (91/76);  new base T (76/70);  total (167/146)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 316 total aligned reads displayed.

GGCAGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATAT  >  CP009273/4584607‑4584905
                                                                                                                                                     |                                                                                                                                                     
ggCAGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTa                                                                                                                                                       <  1:1955308/150‑1 (MQ=255)
  cAGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAAt                                                                                                                                                     <  1:4253274/150‑1 (MQ=255)
   aGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATg                                                                                                                                                    <  2:2760175/150‑1 (MQ=255)
    ggTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGa                                                                                                                                                   >  1:4386481/1‑150 (MQ=255)
     gTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAg                                                                                                                                                  >  1:1075790/1‑150 (MQ=255)
      ttAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGc                                                                                                                                                 <  2:2541725/150‑1 (MQ=255)
        aaaGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCAt                                                                                                                                               <  1:2966434/150‑1 (MQ=255)
           gTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAAc                                                                                                                                            <  2:5075111/150‑1 (MQ=255)
            ttttAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAAcg                                                                                                                                           <  2:4975120/150‑1 (MQ=255)
            ttttAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAAcg                                                                                                                                           <  1:769501/150‑1 (MQ=255)
              ttAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAAcgcg                                                                                                                                         >  2:1578346/1‑150 (MQ=255)
              ttAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAAcgcg                                                                                                                                         >  2:201155/1‑150 (MQ=255)
               tAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCgt                                                                                                                                        <  1:1104244/150‑1 (MQ=255)
               tAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCgt                                                                                                                                        <  2:4585269/150‑1 (MQ=255)
                aGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCgtg                                                                                                                                       >  2:3678790/1‑150 (MQ=255)
                    cTTTACTACCGCTCACAAAAACATTTTTAGTCTACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTtct                                                                                                                                   >  2:1890762/1‑150 (MQ=255)
                     tttACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTtctc                                                                                                                                  <  2:3263980/150‑1 (MQ=255)
                        actacCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCCTAACGCGGGTTCTCAtt                                                                                                                               >  2:3703356/1‑150 (MQ=255)
                          tacCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTaa                                                                                                                             >  2:1366481/1‑150 (MQ=255)
                           acCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAAt                                                                                                                            <  1:3579225/150‑1 (MQ=255)
                              gCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCTGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGa                                                                                                                         <  2:232889/150‑1 (MQ=255)
                               cTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGAt                                                                                                                        <  1:4556135/150‑1 (MQ=255)
                                tCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATc                                                                                                                       <  1:2160428/150‑1 (MQ=255)
                                tCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATc                                                                                                                       <  1:4762515/150‑1 (MQ=255)
                                tCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATc                                                                                                                       <  1:4714541/150‑1 (MQ=255)
                                 cacaAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCg                                                                                                                      >  1:4835658/1‑150 (MQ=255)
                                  acaAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCgn                                                                                                                     <  1:3775244/150‑2 (MQ=255)
                                   caAAAACATTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAAt                                                                                                                   >  1:932699/1‑150 (MQ=255)
                                    aaaaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAAt                                                                                                                   <  1:856340/150‑1 (MQ=255)
                                    aaaaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAAt                                                                                                                   <  2:5031389/150‑1 (MQ=255)
                                     aaaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATc                                                                                                                  >  1:1752713/1‑150 (MQ=255)
                                     aaaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATc                                                                                                                  >  2:4465211/1‑150 (MQ=255)
                                      aaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCa                                                                                                                 >  2:3913719/1‑150 (MQ=255)
                                      aaaCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCa                                                                                                                 >  2:859827/1‑150 (MQ=255)
                                        aCATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCatt                                                                                                               >  1:3883206/1‑150 (MQ=255)
                                         cATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCattg                                                                                                              >  1:4096430/1‑150 (MQ=255)
                                            ttttAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCattgatt                                                                                                           <  1:3708841/150‑1 (MQ=255)
                                              ttAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGt                                                                                                         >  2:3779303/1‑150 (MQ=255)
                                              ttAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGt                                                                                                         >  2:2800915/1‑150 (MQ=255)
                                              ttAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGt                                                                                                         <  2:5121738/150‑1 (MQ=255)
                                                aGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGttt                                                                                                       >  2:1626319/1‑150 (MQ=255)
                                                 gTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTa                                                                                                      <  1:1283540/150‑1 (MQ=255)
                                                 gTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTa                                                                                                      <  1:3349069/150‑1 (MQ=255)
                                                    gACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCt                                                                                                   >  1:2773168/1‑150 (MQ=255)
                                                    gACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCt                                                                                                   >  1:1496067/1‑150 (MQ=255)
                                                    gACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCt                                                                                                   >  1:4895924/1‑150 (MQ=255)
                                                     aCTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTg                                                                                                  >  1:3951524/1‑150 (MQ=255)
                                                     aCTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTg                                                                                                  >  1:2617263/1‑150 (MQ=255)
                                                      cTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGa                                                                                                 <  1:2051034/150‑1 (MQ=255)
                                                       tGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGaa                                                                                                <  2:1293476/150‑1 (MQ=255)
                                                         aTTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAAtt                                                                                              <  1:1865460/150‑1 (MQ=255)
                                                          ttGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTc                                                                                             <  2:3029602/150‑1 (MQ=255)
                                                             tgAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCaaa                                                                                          >  2:1780105/1‑150 (MQ=255)
                                                               aGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATg                                                                                        <  2:1075539/150‑1 (MQ=255)
                                                               aGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATg                                                                                        <  2:2085630/150‑1 (MQ=255)
                                                               aGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATg                                                                                        >  1:5004598/1‑150 (MQ=255)
                                                                 cATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGt                                                                                      >  1:4481307/1‑150 (MQ=255)
                                                                  aTAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGtt                                                                                     <  2:4162011/150‑1 (MQ=255)
                                                                  aTAGGTAAGCGGGAATTACCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAAGGGATCGAATCAGAGATTGGTTAGCTGAATACAAATgggt                                                                                     >  2:3136593/1‑148 (MQ=255)
                                                                   tAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTa                                                                                    <  1:222075/150‑1 (MQ=255)
                                                                       tAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGa                                                                                >  2:1431851/1‑150 (MQ=255)
                                                                            gggAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCa                                                                           >  2:2940882/1‑150 (MQ=255)
                                                                              gAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAg                                                                         <  1:2722765/150‑1 (MQ=255)
                                                                               aaTTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAgg                                                                        >  2:4137422/1‑150 (MQ=255)
                                                                                 tttCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGt                                                                      <  1:1235912/150‑1 (MQ=255)
                                                                                  ttCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGTTTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTg                                                                     >  2:2224479/1‑150 (MQ=255)
                                                                                    cccGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGc                                                                   <  1:2416534/150‑1 (MQ=255)
                                                                                     ccGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCa                                                                  <  2:3523197/150‑1 (MQ=255)
                                                                                      cGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCAt                                                                 <  1:1645596/150‑1 (MQ=255)
                                                                                       ggTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATc                                                                <  2:2821204/150‑1 (MQ=255)
                                                                                               tCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCAc                                                        >  2:468004/1‑150 (MQ=255)
                                                                                                 ttttGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTg                                                      >  1:3131055/1‑150 (MQ=255)
                                                                                                         aaaCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAAt                                              <  2:810915/150‑1 (MQ=255)
                                                                                                         aaaCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAAt                                              <  1:1043363/150‑1 (MQ=255)
                                                                                                         aaaCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAAt                                              <  1:3678790/150‑1 (MQ=255)
                                                                                                          aaCCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATc                                             >  2:639970/1‑150 (MQ=255)
                                                                                                            ccGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCag                                           <  2:1951906/150‑1 (MQ=255)
                                                                                                              gTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCagag                                         >  1:57763/1‑150 (MQ=255)
                                                                                                              gTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCagag                                         >  1:3163949/1‑150 (MQ=255)
                                                                                                                  aCCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCt                                     <  1:725036/150‑1 (MQ=255)
                                                                                                                  aCCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCt                                     <  1:5057398/150‑1 (MQ=255)
                                                                                                                          aTAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTATTGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATgg                             >  1:1083913/1‑150 (MQ=255)
                                                                                                                          aTAGTTTGCTGTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATgg                             <  1:5124098/150‑1 (MQ=255)
                                                                                                                            aGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTc                           >  1:4610384/1‑150 (MQ=255)
                                                                                                                               ttGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCata                        >  1:3478031/1‑150 (MQ=255)
                                                                                                                               ttGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCata                        <  1:3070339/150‑1 (MQ=255)
                                                                                                                                tGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCatat                       <  2:658289/150‑1 (MQ=255)
                                                                                                                                   tGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCt                    <  1:859827/150‑1 (MQ=255)
                                                                                                                                    gTTTTTAATCCGTACTTTATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTc                   <  2:430548/150‑1 (MQ=255)
                                                                                                                                    gTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTc                   >  2:1512939/1‑150 (MQ=255)
                                                                                                                                    gTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTc                   <  2:2040953/150‑1 (MQ=255)
                                                                                                                                    gTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTc                   >  2:1127543/1‑150 (MQ=255)
                                                                                                                                      atttAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCAt                 <  2:4906727/149‑1 (MQ=255)
                                                                                                                                         tAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGa              >  1:1445849/1‑150 (MQ=255)
                                                                                                                                              cGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCa         <  2:3581321/150‑1 (MQ=255)
                                                                                                                                                  cTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGa     >  1:4308786/1‑150 (MQ=255)
                                                                                                                                                   ttAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGat    >  1:3382143/1‑150 (MQ=255)
                                                                                                                                                    tAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGata   >  1:2326771/1‑150 (MQ=255)
                                                                                                                                                    tAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGata   <  1:1019009/150‑1 (MQ=255)
                                                                                                                                                     aaTGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGatat  >  2:4009275/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCAGGTTAAAGTTTTAGGACTTTACTACCGCTCACAAAAACATTTTTAGTCGACTGATTGTGAGCATAGGTAAGCGGGAATTTCCCGGTATTTATCTTTTGATAAAACCGTGCACCCGCTAATAGTTTGCTGTTTTTAATCCGTACTTAATGAGCATAACGCGTGTTCTCATTAATGGATCGAATCATTGATTGTTTAGCTGAATTCAAATGGTTACGGAGCGCAAGGGTGGCATCAAGGTCACTGCGGCAAATCAGAGCGCTGAGAATGGTCATATGCTCATCGACGGCAATGATAT  >  CP009273/4584607‑4584905

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: