Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,938,640 C→A 5.6% 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.6% 608.8 / 10.4 159D188Y (GAT→TAT) tcdAtRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA
Reads supporting (aligned to +/- strand):  ref base C (77/73);  new base A (5/4);  total (82/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 164 total aligned reads displayed.

CAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCAAT  >  CP009273/2938493‑2938787
                                                                                                                                                   |                                                                                                                                                   
cAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATc                                                                                                                                                     >  2:890796/1‑148 (MQ=255)
 aGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCg                                                                                                                                                    <  1:2752777/148‑1 (MQ=255)
   aTGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCt                                                                                                                                                  >  2:3021987/1‑148 (MQ=255)
   aTGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCt                                                                                                                                                  <  1:1239968/148‑1 (MQ=255)
   aTGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCt                                                                                                                                                  <  1:310036/148‑1 (MQ=255)
    tGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCtt                                                                                                                                                 <  2:2599254/148‑1 (MQ=255)
      cacaATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCtttt                                                                                                                                               >  2:2467408/1‑148 (MQ=255)
       acaATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTc                                                                                                                                              <  1:1152939/148‑1 (MQ=255)
         aaTCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAg                                                                                                                                            >  1:3175393/1‑148 (MQ=255)
            ccatccGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGgcg                                                                                                                                         <  1:1336028/148‑1 (MQ=255)
                ccGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTcg                                                                                                                                     >  1:571950/1‑148 (MQ=255)
                ccGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTcg                                                                                                                                     >  1:2185304/1‑148 (MQ=255)
                      cGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAAc                                                                                                                               <  1:890796/148‑1 (MQ=255)
                          ccTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTcg                                                                                                                           >  2:1633640/1‑148 (MQ=255)
                           cTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTcgc                                                                                                                          >  1:1758774/1‑148 (MQ=255)
                             tCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTcgccg                                                                                                                        >  1:306425/1‑148 (MQ=255)
                              cTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTcgccgc                                                                                                                       <  1:3037271/148‑1 (MQ=255)
                                     ggCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCggg                                                                                                                >  2:1242264/1‑148 (MQ=255)
                                     ggCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCggg                                                                                                                >  2:109179/1‑148 (MQ=255)
                                         ttCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCt                                                                                                            >  1:1917444/1‑148 (MQ=255)
                                           cATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGa                                                                                                          <  1:2746247/148‑1 (MQ=255)
                                            aTCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGaa                                                                                                         >  1:1354892/1‑148 (MQ=255)
                                            aTCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGaa                                                                                                         >  2:1923686/1‑148 (MQ=255)
                                             tCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAAt                                                                                                        <  2:1488960/148‑1 (MQ=255)
                                               gCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCg                                                                                                      <  1:2372723/148‑1 (MQ=255)
                                                  cacacCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGttt                                                                                                   <  1:2101718/148‑1 (MQ=255)
                                                  cacacCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGttt                                                                                                   <  1:1149682/148‑1 (MQ=255)
                                                  cacacCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGttt                                                                                                   <  1:1021315/148‑1 (MQ=255)
                                                    cacCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTc                                                                                                 >  2:1445635/1‑148 (MQ=255)
                                                     accgtaccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCg                                                                                                <  2:2489732/148‑1 (MQ=255)
                                                       cgtaccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGcc                                                                                              >  1:284337/1‑148 (MQ=255)
                                                       cgtaccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGcc                                                                                              >  1:284437/1‑148 (MQ=255)
                                                          accgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGa                                                                                           >  1:641347/1‑148 (MQ=255)
                                                          accgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGa                                                                                           >  2:2082966/1‑148 (MQ=255)
                                                           ccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGAt                                                                                          >  1:8992/1‑148 (MQ=255)
                                                           ccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGAt                                                                                          >  1:2812548/1‑148 (MQ=255)
                                                           ccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGAt                                                                                          >  2:1207551/1‑148 (MQ=255)
                                                                    tGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCt                                                                                 >  2:3011233/1‑148 (MQ=255)
                                                                        ggTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTTGCCAGATCGGTAACCTGAAt                                                                             <  2:2083586/148‑1 (MQ=255)
                                                                        ggTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTTGCCAGATCGGTAACCTGAAt                                                                             <  2:2083246/148‑1 (MQ=255)
                                                                           acacCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTg                                                                          <  1:2918834/148‑1 (MQ=255)
                                                                               cAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTc                                                                      >  1:771192/1‑148 (MQ=255)
                                                                               cAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTc                                                                      >  2:182573/1‑148 (MQ=255)
                                                                               cAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTc                                                                      >  2:182579/1‑148 (MQ=255)
                                                                               aaGCTCTTCAGTAGAAAACACGCAATCCACGCCGAGCTAACATTTACTCTTTATCACTACGACAAAATCGCTTTACAGTCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTAGCCATAACGGTAACCTGAATCTGCGTc                                                                      >  2:2207564/2‑148 (MQ=255)
                                                                                aGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCg                                                                     >  1:1389531/1‑148 (MQ=255)
                                                                                aGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCg                                                                     >  1:1389495/1‑148 (MQ=255)
                                                                                  cgcTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGa                                                                   <  1:2106309/148‑1 (MQ=255)
                                                                                    cTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATc                                                                 >  2:1263475/1‑148 (MQ=255)
                                                                                       cAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAAt                                                              <  2:1870047/148‑1 (MQ=255)
                                                                                          tAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTg                                                           >  2:669281/1‑148 (MQ=255)
                                                                                          tAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTg                                                           >  1:1422265/1‑148 (MQ=255)
                                                                                           aGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGc                                                          <  1:2891347/148‑1 (MQ=255)
                                                                                             aaaaCACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGACTCTGCGTCGGATCAATCTGccc                                                        >  1:1599931/1‑148 (MQ=255)
                                                                                             aaaaCACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGccc                                                        >  1:1599909/1‑148 (MQ=255)
                                                                                             aaaaCACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGccc                                                        >  2:122679/1‑148 (MQ=255)
                                                                                               aaCACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCAc                                                      >  2:189779/1‑148 (MQ=255)
                                                                                               aaCACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCAc                                                      <  2:615665/148‑1 (MQ=255)
                                                                                               aaCACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCAc                                                      <  2:615543/148‑1 (MQ=255)
                                                                                                       aTCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCa                                              <  1:1273077/148‑1 (MQ=255)
                                                                                                           aCGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCgg                                          <  1:2027122/148‑1 (MQ=255)
                                                                                                                 aGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGa                                    <  2:752472/148‑1 (MQ=255)
                                                                                                                 aGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGa                                    <  2:752462/148‑1 (MQ=255)
                                                                                                                 aGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGa                                    >  2:1352763/1‑148 (MQ=255)
                                                                                                                  gCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGAc                                   <  1:2209587/148‑1 (MQ=255)
                                                                                                                  gCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAAATTCGCCGCCAGCGGGTCCTGAATCGTTTTCGACAGATCGGTAACATGAATCTGAGTCGGATCAATCTGCCCACCCGCGGCACCGGTTGTGAc                                   <  1:2207564/148‑1 (MQ=255)
                                                                                                                   cTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGAcc                                  >  1:2921624/1‑148 (MQ=255)
                                                                                                                   cTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGAcc                                  <  1:418614/148‑1 (MQ=255)
                                                                                                                   cTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGAcc                                  >  2:3017670/1‑148 (MQ=255)
                                                                                                                   cTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGAcc                                  <  1:2083814/148‑1 (MQ=255)
                                                                                                                   cTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGAcc                                  <  1:2083668/148‑1 (MQ=255)
                                                                                                                    ttACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCa                                 >  1:425553/1‑148 (MQ=255)
                                                                                                                      aCCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGt                               <  1:1657484/148‑1 (MQ=255)
                                                                                                                      aCCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGt                               >  1:96863/1‑148 (MQ=255)
                                                                                                                      aCCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGt                               >  2:2640630/1‑148 (MQ=255)
                                                                                                                        cTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTgg                             >  2:2168264/1‑148 (MQ=255)
                                                                                                                           tACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGAt                          <  2:1972760/148‑1 (MQ=255)
                                                                                                                            aCTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGAtt                         <  1:592170/148‑1 (MQ=255)
                                                                                                                              tGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGAtttt                       <  1:1929488/148‑1 (MQ=255)
                                                                                                                                   ttCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTgc                  >  2:3044245/1‑148 (MQ=255)
                                                                                                                                    tCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTgcg                 <  2:1354892/148‑1 (MQ=255)
                                                                                                                                     cACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTgcgc                >  2:1061005/1‑148 (MQ=255)
                                                                                                                                       cTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCg              <  1:104187/148‑1 (MQ=255)
                                                                                                                                       cTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCg              >  1:523293/1‑148 (MQ=255)
                                                                                                                                          cGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACa           <  1:3021987/148‑1 (MQ=255)
                                                                                                                                           gCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACaa          <  2:2206827/148‑1 (MQ=255)
                                                                                                                                              aaaaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATaa       <  1:35004/148‑1 (MQ=255)
                                                                                                                                               aaaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAg      >  1:1900341/1‑148 (MQ=255)
                                                                                                                                                aaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGc     >  1:2381947/1‑148 (MQ=255)
                                                                                                                                                aaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGc     >  2:310592/1‑148 (MQ=255)
                                                                                                                                                aaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGc     >  2:310608/1‑148 (MQ=255)
                                                                                                                                                aaTAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGc     >  2:2967362/1‑148 (MQ=255)
                                                                                                                                                 aTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCa    <  1:1633640/148‑1 (MQ=255)
                                                                                                                                                 aTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCa    <  1:2946257/148‑1 (MQ=255)
                                                                                                                                                  tCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCaa   <  2:399086/148‑1 (MQ=255)
                                                                                                                                                  tCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCaa   <  1:3066861/148‑1 (MQ=255)
                                                                                                                                                  tCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCaa   <  2:2063317/148‑1 (MQ=255)
                                                                                                                                                   cGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCAAt  >  1:2598193/1‑148 (MQ=255)
                                                                                                                                                   cGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCAAt  <  1:3187771/148‑1 (MQ=255)
                                                                                                                                                   cGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCAAt  <  1:799812/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCAAT  >  CP009273/2938493‑2938787

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: