Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,151,733 A→T 8.8% D242V (GAT→GTT)  yqhG → DUF3828 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,151,7330AT 8.8% 769.5 / 9.1 262D242V (GAT→GTT) yqhGDUF3828 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (93/142);  new base T (22/2);  total (118/144)
Fisher's exact test for biased strand distribution p-value = 8.50e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 274 total aligned reads displayed.

CGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGC  >  CP009273/3151587‑3151881
                                                                                                                                                  |                                                                                                                                                    
cGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTtgatg                                                                                                                                                    <  1:1014546/149‑1 (MQ=255)
 gACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAATATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTtgatga                                                                                                                                                   >  2:599977/1‑149 (MQ=255)
 gACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTtgatga                                                                                                                                                   <  1:3508977/149‑1 (MQ=255)
  aCTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGaa                                                                                                                                                  <  1:3624015/149‑1 (MQ=255)
  aCTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGaa                                                                                                                                                  <  2:2302003/149‑1 (MQ=255)
    tctTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAAt                                                                                                                                                <  2:3137382/149‑1 (MQ=255)
     ctTACCGCCAGGCGGTTAATGCCTTCTATACAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAAtt                                                                                                                                               <  1:2779524/149‑1 (MQ=255)
      ttACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTGACGCGTTGATGAAAATc                                                                                                                                              >  2:814089/1‑149 (MQ=255)
      ttACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTc                                                                                                                                              <  2:4444307/149‑1 (MQ=255)
      ttACCGAGAGGCGGTTATTTACTTCTATAAATCCACATTCTTGACGTCAAATAATCACGATGAATATTTGCAATTATAACGGCAATTTATTTACTACAATCCGATCTGTTTTGGCGAAGAGTCTCTATTTTCACGAGTTGTTTAAATTc                                                                                                                                              <  2:650357/149‑1 (MQ=255)
       tACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAAGTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCg                                                                                                                                             >  2:4679457/1‑149 (MQ=255)
        aCCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCGATTTTCACGCGTTGATGAAATTCGa                                                                                                                                            >  2:434197/1‑149 (MQ=255)
         ccGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAg                                                                                                                                           >  2:1835896/1‑149 (MQ=255)
         ccGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAg                                                                                                                                           <  1:1327920/149‑1 (MQ=255)
         ccGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAg                                                                                                                                           <  1:3273912/149‑1 (MQ=255)
         ccGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGGTGATGAAATTCGAg                                                                                                                                           >  1:2263461/1‑149 (MQ=255)
            ccAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTcc                                                                                                                                        <  2:3230046/149‑1 (MQ=255)
             cAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCt                                                                                                                                       <  2:545055/149‑1 (MQ=255)
               ggCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTAACGCGTTGATGTAAATAGAGTCCTgg                                                                                                                                     >  1:986094/1‑149 (MQ=255)
                gCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACCCGTTGATGAAATTAGAGTCCAGga                                                                                                                                    >  1:3429018/1‑149 (MQ=255)
                   gTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGaa                                                                                                                                 <  1:85539/149‑1 (MQ=255)
                    ttAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTAACGCGTTGATGAAATTAGAGTACAGGGGaaa                                                                                                                                >  1:4167477/1‑149 (MQ=255)
                    ttAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCGATTTTGACGCGGTGGTGAAAATAGAGTCCAGGAGaaa                                                                                                                                >  1:162878/1‑149 (MQ=255)
                       aTGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAAcc                                                                                                                             >  2:1480649/1‑149 (MQ=255)
                         gCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAAccgc                                                                                                                           <  1:1074812/149‑1 (MQ=255)
                         gCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAAccgc                                                                                                                           <  1:1074841/149‑1 (MQ=255)
                         gCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTAACGCGTGTATGAAATTCGAGTCCGGGGGAAAAacgc                                                                                                                           >  2:2365891/1‑149 (MQ=255)
                          ccTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATTTAATTAGAGTCCAGGAGAAAAacgcc                                                                                                                          >  1:3088145/1‑149 (MQ=255)
                            ttCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTAGAGTCCAGGATAAAACCGCCGa                                                                                                                        >  2:2933446/1‑149 (MQ=255)
                              cTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGAtt                                                                                                                      <  2:656784/149‑1 (MQ=255)
                                ataAAACCACCTTCAAGACGGCCAATAATCCCGATGAATATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACCCGTTGATGAAAATCGAGTCCTGGAGAAAAACGCCCAttac                                                                                                                    >  1:3473259/1‑147 (MQ=255)
                                ataAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTcc                                                                                                                    <  1:481297/149‑1 (MQ=255)
                                ataAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTcc                                                                                                                    <  2:2511529/149‑1 (MQ=255)
                                  aaaaCCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTccgc                                                                                                                  <  2:4225223/149‑1 (MQ=255)
                                      ccaccTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTccgcccgc                                                                                                              <  2:4553234/149‑1 (MQ=255)
                                          cTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTc                                                                                                          >  2:3817326/1‑149 (MQ=255)
                                           ttCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCa                                                                                                         <  1:2096718/149‑1 (MQ=255)
                                           ttCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTAACGCGGTGATGAAATTAGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCAATCa                                                                                                         >  2:4098464/1‑149 (MQ=255)
                                              aaGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTAACGCGGGGTTGAAAATAGAGTCCTGGAGTAAACCCCCGAGTCCCCCCGCATTCATGt                                                                                                      >  1:3566758/1‑149 (MQ=255)
                                               aGACGGCCAATAATCCCGATGAATATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGtt                                                                                                     <  1:4082175/149‑1 (MQ=255)
                                               aGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTAACGCGGTGATTAAATTAGATTCCAGGCGAAAAACCCCCCTTACGCCCCCAATCAAGtt                                                                                                     >  1:89565/1‑149 (MQ=255)
                                                gACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTc                                                                                                    <  2:2537125/149‑1 (MQ=255)
                                                gACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTc                                                                                                    <  2:195000/149‑1 (MQ=255)
                                                  cGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGt                                                                                                  <  2:2072928/149‑1 (MQ=255)
                                                   ggCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCAGGGGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGtt                                                                                                 <  2:2984783/149‑1 (MQ=255)
                                                   ggCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGGTGTTGAAATTAGAGTCCATGAGAAAACCCCCGATTCCGCCCGCATTAATGTTCggt                                                                                                 >  2:3551917/1‑147 (MQ=255)
                                                    gCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGttt                                                                                                <  2:1727964/149‑1 (MQ=255)
                                                      cAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCCCCGATTCCGCCCGCATTCAAGTTAGTTTTa                                                                                              >  2:3397494/1‑149 (MQ=255)
                                                       aataatCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGGTGATGAAAATCGAGTCCGGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTAc                                                                                             >  2:246022/1‑149 (MQ=255)
                                                        ataatCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATTAAAATAGAGTCCGGGAGAAAACCGCCGATTACGCCCGCAATCATGTTAGTTTTACg                                                                                            >  2:3606399/1‑149 (MQ=255)
                                                              ccGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTTTTTTCACGCGGTGTTGAAATTAGAGTCCTGGAGAAAACCGCCCATTACGCCCCCAATAAAGTTCGTTTTACGCTGaac                                                                                      >  1:4588393/1‑147 (MQ=255)
                                                              ccGATGAAGAGATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGAcc                                                                                      <  1:2636248/149‑1 (MQ=255)
                                                                gATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCaa                                                                                    <  1:4556019/149‑1 (MQ=255)
                                                                  tGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGTTGAAAATCGAGTACTGGAGAAAAACCCCGAATACGCCCGCATTAAAGTTAGTTTTAAGCAGACCAAtt                                                                                  >  1:341239/1‑148 (MQ=255)
                                                                    aaGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGc                                                                                <  2:2101551/149‑1 (MQ=255)
                                                                    aaGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTAGAGTCCAGGAGAAAACCGCCGATTCCGCCCGCATTAATGTTCGTTTTACGCTGACCAAAGGc                                                                                >  2:748243/1‑149 (MQ=255)
                                                                     aGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTGACGCGTTGTTGAAATTAGAGTCCTGGGGAAAACCGCCCAATCCACCCCCAATAAAGTTAGTTTTTAGCTGTCCAAAGGca                                                                               >  1:1447181/1‑149 (MQ=255)
                                                                        tatGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGcaaca                                                                            >  2:1060603/1‑149 (MQ=255)
                                                                          tGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCaacaac                                                                          <  1:2809138/149‑1 (MQ=255)
                                                                           gCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAAATCGAGTCCTGGAGAAAAACGCCGATTCCCCCCGCCTTCCTGTTCGTTTTTCGCTGACCAATGGCAACAACg                                                                         >  2:2816214/1‑149 (MQ=255)
                                                                             aaaTAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGaa                                                                       <  1:1534655/149‑1 (MQ=255)
                                                                                 aGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCCGTGATGAAATTCGAGTCCAGGAGAAAAACCCCCAGTCCCCCCCCATTCAAGGTCGTTTTACGCTGACCACTGGCCACAACCAACAGc                                                                   >  2:3322621/1‑149 (MQ=255)
                                                                                   aaCGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGTTGAAATTCGAGTCCTGGAGTAAACCGCCGATTCCGCCCGCAATCATGTTCGTTTTAAGCTGACCAATGGCAACCACGAAGAGCaa                                                                 >  2:1087559/1‑149 (MQ=255)
                                                                                   aaCGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCaa                                                                 <  1:552980/149‑1 (MQ=255)
                                                                                         aaTTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTAACGCGGTGATTAAAATAGAGTCCTGGCGAAAACCCCCCAATCCCCCCCCAATAATGTTCGTTTTAAGCAGACCAAATGCAACAACCAAGAGCAAGAAcag                                                           >  1:4130913/1‑147 (MQ=255)
                                                                                         aaTTTATTTACGACAATCCGATCTGTCTTGGCGAAGAGTCGCTTTTTTAACGCGGTGATTAAAATAGAGTCCAGGTGAAAAACGCCGAATACGCCCGCAATCATGTTCGTTTTAAGCGGACCAAATGCAACAACGAAGAGCCAGAAcag                                                           >  1:1888134/1‑147 (MQ=255)
                                                                                          atttatttACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTgg                                                          <  2:2425046/149‑1 (MQ=255)
                                                                                             tatttaCGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGttt                                                       <  1:1480649/149‑1 (MQ=255)
                                                                                              atttaCGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGtttt                                                      <  2:4240143/149‑1 (MQ=255)
                                                                                              atttaCGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTAACGCGGTGGTGAAAATATAGTCCTGGAGAAAAACCCCCATTACGCCCGCCAACATGTTCGTTTTTAGCAGACCAATGGCAACAACAAAGAGCAAGAAAAGGtttt                                                      >  1:2613351/1‑149 (MQ=255)
                                                                                                ttaCGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTAc                                                    <  1:442200/149‑1 (MQ=255)
                                                                                                ttaCGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTAc                                                    <  1:1226995/149‑1 (MQ=255)
                                                                                                 taCGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGGTGATGAAATTAGAGTCCTGGAGAAAACCGCCGCTTACGCCCGCATTCAAGTTCGTTTTACGCTGACCAAAGGCAACAACCAAGAGCAAGAACAGGTTTTACa                                                   >  2:2631155/1‑149 (MQ=255)
                                                                                                     aCAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCAGGAGAAAACCGCCGATTACGCCCCCATTCATGTTAGTTTTAAGCTGACCCAAGTCAACAACCAAGAGCAAGAACTGGTTTTAAAgcgg                                               >  2:149062/1‑149 (MQ=255)
                                                                                                        aTCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGgcgc                                            <  2:4054162/149‑1 (MQ=255)
                                                                                                        aTCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGgcgc                                            <  2:4280543/149‑1 (MQ=255)
                                                                                                               cTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCa                                     <  2:4588428/149‑1 (MQ=255)
                                                                                                                 gTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCCCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAAAGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGACGGCAAg                                   >  2:3418883/1‑149 (MQ=255)
                                                                                                                     tGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTggg                               <  1:4260492/149‑1 (MQ=255)
                                                                                                                     tGGCGAAGAGTCGCTATTTTCACGCGGTGATGAAATTAGAGTCCTGGAGAAAACCGCCGATTACGCCCGCATTCAAGTTCGTTTTACGCTGCACAATTGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTggg                               >  2:128332/1‑149 (MQ=255)
                                                                                                                     tGGCGAAGAGTCGCTATTTTCACGCCTGGTTGAAAATCGACTCCTGGAGAAAAACCCCCATTACCCCCGCAATAAAGTTAGGTTTAACCAGACCAAATGCAACAAACAAGATCAAGAACTGGTTTTACACCGGCGCGAACGCCAGTggg                               >  2:1325858/1‑149 (MQ=255)
                                                                                                                      ggCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTAGAGTCCTGGAGAAAACCGCCGCTTCCGCCCGCAATCATGTTCGTTTTAAGCTGACCAAATGCAACAACGAAGAGCAAGAACAGGTTTTAAAGCGGCGCGAAGGCCAGTGGGa                              >  2:4632107/1‑149 (MQ=255)
                                                                                                                      ggCGAAGAGTCGCTATTTTCACGCGGTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCAATCATGTTCGTTTTACGCAGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGa                              >  2:2721148/1‑149 (MQ=255)
                                                                                                                         gAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAGCCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAAt                           <  2:3284569/149‑1 (MQ=255)
                                                                                                                         gAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAAt                           <  1:2399955/149‑1 (MQ=255)
                                                                                                                         gAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAAt                           <  2:172978/149‑1 (MQ=255)
                                                                                                                              gTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTg                      >  2:3770760/1‑149 (MQ=255)
                                                                                                                              gTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTg                      >  1:611797/1‑149 (MQ=255)
                                                                                                                                   tATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGAtttt                 >  2:3832617/1‑149 (MQ=255)
                                                                                                                                   tATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGAtttt                 <  2:1347905/149‑1 (MQ=255)
                                                                                                                                   tATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGAtttt                 <  2:1460072/149‑1 (MQ=255)
                                                                                                                                    aTTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTa                <  1:714165/149‑1 (MQ=255)
                                                                                                                                    aTTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTa                <  2:1629159/149‑1 (MQ=255)
                                                                                                                                     ttttCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTAt               >  2:2267246/1‑149 (MQ=255)
                                                                                                                                      tttCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATc              <  1:3976791/149‑1 (MQ=255)
                                                                                                                                         cACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGAACAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGc           >  1:1234869/1‑149 (MQ=255)
                                                                                                                                          aCGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCACCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGcc          <  1:3011043/149‑1 (MQ=255)
                                                                                                                                                 gatgaAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAg   >  1:94084/1‑149 (MQ=255)
                                                                                                                                                  atgaAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGc  <  2:3429018/149‑1 (MQ=255)
                                                                                                                                                  atgaAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGc  >  1:3932511/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGACTCTTACCGCCAGGCGGTTAATGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCCTGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGC  >  CP009273/3151587‑3151881

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑