Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,938,640 C→A 5.1% D188Y (GAT→TAT)  tcdA ← tRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,938,6400CA 5.1% 1064.0 / 19.8 273D188Y (GAT→TAT) tcdAtRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA
Reads supporting (aligned to +/- strand):  ref base C (134/124);  new base A (11/4);  total (145/128)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 274 total aligned reads displayed.

TCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCAATC  >  CP009273/2938491‑2938788
                                                                                                                                                     |                                                                                                                                                    
   aGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCt                                                                                                                                                   >  2:197254/1‑150 (MQ=255)
    gATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCtt                                                                                                                                                  <  1:3635578/150‑1 (MQ=255)
     aTGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCttt                                                                                                                                                 <  1:1963749/150‑1 (MQ=255)
       gCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTc                                                                                                                                               <  2:2705529/150‑1 (MQ=255)
        cacaATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCa                                                                                                                                              >  2:1116323/1‑150 (MQ=255)
         acaATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAg                                                                                                                                             <  2:56418/150‑1 (MQ=255)
          caATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAgg                                                                                                                                            >  2:3557821/1‑150 (MQ=255)
          caATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAgg                                                                                                                                            >  2:1362857/1‑150 (MQ=255)
          caATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCGTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTGTCAgg                                                                                                                                            >  2:1789312/1‑150 (MQ=255)
           aaTCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGgc                                                                                                                                           >  1:63157/1‑150 (MQ=255)
              ccatccGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCt                                                                                                                                        <  2:4005713/150‑1 (MQ=255)
               catccGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTc                                                                                                                                       <  1:1997820/150‑1 (MQ=255)
                   cGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCgcgc                                                                                                                                   >  2:621665/1‑150 (MQ=255)
                       tCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACt                                                                                                                               >  1:1676372/1‑150 (MQ=255)
                        cGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACtt                                                                                                                              <  1:1048136/150‑1 (MQ=255)
                         ggCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTc                                                                                                                             <  1:3215291/150‑1 (MQ=255)
                         ggCCCTTCTGCCGTGGCTTTCATCGCACACACCGGACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTc                                                                                                                             <  2:2087862/150‑1 (MQ=255)
                          gCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTcg                                                                                                                            <  1:4397246/150‑1 (MQ=255)
                            ccTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTcgcc                                                                                                                          >  2:2433133/1‑150 (MQ=255)
                             cTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTcgccg                                                                                                                         <  2:1484771/150‑1 (MQ=255)
                               tCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCgccgcc                                                                                                                       <  2:4151445/150‑1 (MQ=255)
                                 tGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAg                                                                                                                     <  2:2266352/150‑1 (MQ=255)
                                 tGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAg                                                                                                                     >  2:917429/1‑150 (MQ=255)
                                  gCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGc                                                                                                                    >  1:3911788/1‑150 (MQ=255)
                                       ggCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTc                                                                                                               <  2:3614373/150‑1 (MQ=255)
                                       ggCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTc                                                                                                               >  1:2984527/1‑150 (MQ=255)
                                        gCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCt                                                                                                             >  1:3960066/1‑150 (MQ=255)
                                         cTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCt                                                                                                             >  2:2169548/1‑150 (MQ=255)
                                            tCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGaa                                                                                                          <  1:256027/150‑1 (MQ=255)
                                              aTCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATc                                                                                                        >  1:4502131/1‑150 (MQ=255)
                                                 gCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGtt                                                                                                     >  1:465778/1‑150 (MQ=255)
                                                 gCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGtt                                                                                                     <  1:4146269/150‑1 (MQ=255)
                                                 gCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGtt                                                                                                     >  1:3864982/1‑150 (MQ=255)
                                                  cacacaCCGTACCGTCTGACTGCGGGTACACCATCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGttt                                                                                                    <  2:1062396/150‑1 (MQ=255)
                                                  cacacaCCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGttt                                                                                                    <  1:2207445/150‑1 (MQ=255)
                                                  cacacaCCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGttt                                                                                                    <  2:3783881/150‑1 (MQ=255)
                                                  cacacaCCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGttt                                                                                                    >  1:827605/1‑150 (MQ=255)
                                                   acacacCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGtttt                                                                                                   <  2:916240/150‑1 (MQ=255)
                                                   acacacCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGtttt                                                                                                   <  2:1680926/150‑1 (MQ=255)
                                                    cacacCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAATTTCGCCGCCAGCGGGTCCTGAATCGTTTTc                                                                                                  <  2:760041/150‑1 (MQ=255)
                                                    cacacCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTc                                                                                                  <  1:2812011/150‑1 (MQ=255)
                                                        ccgtaccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCa                                                                                              <  2:1445803/150‑1 (MQ=255)
                                                         cgtaccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAg                                                                                             >  1:916684/1‑150 (MQ=255)
                                                           taccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGAt                                                                                           >  1:4525402/1‑150 (MQ=255)
                                                                    aCTGCGGGTACACCAGCGCTTCCGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCt                                                                                  <  1:95267/150‑1 (MQ=255)
                                                                       gCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGaa                                                                               <  1:3337977/150‑1 (MQ=255)
                                                                            tACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGc                                                                          >  1:571568/1‑150 (MQ=255)
                                                                             acacCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCg                                                                         <  1:3683285/150‑1 (MQ=255)
                                                                              cacCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGt                                                                        <  1:3776849/150‑1 (MQ=255)
                                                                                 cAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCgg                                                                     <  1:136405/150‑1 (MQ=255)
                                                                                 cAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCgg                                                                     >  1:2080014/1‑150 (MQ=255)
                                                                                     gcTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCa                                                                 >  2:848363/1‑150 (MQ=255)
                                                                                     gcTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCa                                                                 <  1:2062190/150‑1 (MQ=255)
                                                                                     gcTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCa                                                                 <  1:1484934/150‑1 (MQ=255)
                                                                                      cTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCaa                                                                <  1:3481531/150‑1 (MQ=255)
                                                                                       ttCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAAt                                                               >  2:1868200/1‑150 (MQ=255)
                                                                                        tCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATc                                                              >  1:4455880/1‑150 (MQ=255)
                                                                                           gTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGc                                                           >  2:655804/1‑150 (MQ=255)
                                                                                             aGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGccc                                                         >  2:834261/1‑150 (MQ=255)
                                                                                             aGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGccc                                                         <  2:206979/150‑1 (MQ=255)
                                                                                                 aaCACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCAccc                                                     <  1:3421029/150‑1 (MQ=255)
                                                                                                      gCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGcc                                                <  1:3023183/150‑1 (MQ=255)
                                                                                                      gCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGcc                                                <  2:4639957/150‑1 (MQ=255)
                                                                                                          tCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCg                                            >  1:1784044/1‑150 (MQ=255)
                                                                                                            cACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGt                                          >  1:2482913/1‑150 (MQ=255)
                                                                                                               gCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTtgt                                       <  2:3340125/150‑1 (MQ=255)
                                                                                                                 cGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGa                                     <  1:2421268/150‑1 (MQ=255)
                                                                                                                  gAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGACGGATCAATCTGCCCACCCGCGCCACCGGTTGTGAc                                    >  1:2626708/1‑150 (MQ=255)
                                                                                                                   aGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGAcc                                   >  2:1575359/1‑150 (MQ=255)
                                                                                                                    gCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCa                                  >  1:3619190/1‑150 (MQ=255)
                                                                                                                     cTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAg                                 >  1:1616236/1‑150 (MQ=255)
                                                                                                                      ttACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGt                                <  2:1270321/150‑1 (MQ=255)
                                                                                                                      ttACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGt                                >  2:3602170/1‑150 (MQ=255)
                                                                                                                        aCCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTgg                              >  1:3527649/1‑150 (MQ=255)
                                                                                                                        aCCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTgg                              >  2:238310/1‑150 (MQ=255)
                                                                                                                         ccTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTggg                             >  1:1358952/1‑150 (MQ=255)
                                                                                                                          cTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGa                            >  1:2342670/1‑150 (MQ=255)
                                                                                                                            ttACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGAtt                          >  2:4321226/1‑150 (MQ=255)
                                                                                                                              aCTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGAtttt                        >  1:3197782/1‑150 (MQ=255)
                                                                                                                              aCTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGAtttt                        <  2:3671241/150‑1 (MQ=255)
                                                                                                                              aCTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGAtttt                        <  2:3848293/150‑1 (MQ=255)
                                                                                                                               cTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTa                       <  2:3298081/150‑1 (MQ=255)
                                                                                                                                  tttttCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTg                    >  1:1421642/1‑150 (MQ=255)
                                                                                                                                     ttCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTgcgc                 <  1:3830849/150‑1 (MQ=255)
                                                                                                                                     ttCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTgcgc                 >  2:3709331/1‑150 (MQ=255)
                                                                                                                                      tCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGcc                <  1:1091795/150‑1 (MQ=255)
                                                                                                                                      tCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGcc                >  2:3098425/1‑150 (MQ=255)
                                                                                                                                       cACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCg               >  2:116720/1‑150 (MQ=255)
                                                                                                                                        aCTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGa              <  2:3020424/150‑1 (MQ=255)
                                                                                                                                         cTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGAc             <  2:4543495/150‑1 (MQ=255)
                                                                                                                                         cTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGAc             <  1:3935765/150‑1 (MQ=255)
                                                                                                                                          tACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACa            >  1:2826471/1‑150 (MQ=255)
                                                                                                                                           aCGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACaa           >  2:3988566/1‑150 (MQ=255)
                                                                                                                                              ccAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATaa        >  2:2415596/1‑150 (MQ=255)
                                                                                                                                               cAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAg       >  2:3625818/1‑150 (MQ=255)
                                                                                                                                                aaaaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGc      >  2:4468057/1‑150 (MQ=255)
                                                                                                                                                aaaaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGc      >  2:4468092/1‑150 (MQ=255)
                                                                                                                                                aaaaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGc      >  2:3332156/1‑150 (MQ=255)
                                                                                                                                                aaaaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGc      <  2:2871157/150‑1 (MQ=255)
                                                                                                                                                  aaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCaa    >  2:4102864/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCAATC  >  CP009273/2938491‑2938788

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: