Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,938,640 C→A 16.9% 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,6400CA16.9% 772.1 / 101.3 252D188Y (GAT→TAT) tcdAtRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA
Reads supporting (aligned to +/- strand):  ref base C (118/91);  new base A (18/25);  total (136/116)
Fisher's exact test for biased strand distribution p-value = 9.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

TCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCAATCA  >  CP009273/2938491‑2938789
                                                                                                                                                     |                                                                                                                                                     
tCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATc                                                                                                                                                       >  2:367227/1‑150 (MQ=255)
tCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATc                                                                                                                                                       >  2:617693/1‑150 (MQ=255)
tCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATc                                                                                                                                                       >  1:3072375/1‑150 (MQ=255)
   aGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCt                                                                                                                                                    <  2:3930873/150‑1 (MQ=255)
      tGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCtttt                                                                                                                                                 >  2:2631723/1‑150 (MQ=255)
          caATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAgg                                                                                                                                             <  1:611275/150‑1 (MQ=255)
           aaTCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGgc                                                                                                                                            >  1:4492689/1‑150 (MQ=255)
            atccatccGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGgcg                                                                                                                                           >  2:665343/1‑150 (MQ=255)
              ccatccGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCt                                                                                                                                         >  2:3716983/1‑150 (MQ=255)
                atccGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTcg                                                                                                                                       >  2:1093882/1‑150 (MQ=255)
                atccGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTcg                                                                                                                                       >  1:3513866/1‑150 (MQ=255)
                atccGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTcg                                                                                                                                       <  1:3971770/150‑1 (MQ=255)
                    gCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCa                                                                                                                                   <  1:784804/150‑1 (MQ=255)
                    gCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCa                                                                                                                                   >  2:5079436/1‑150 (MQ=255)
                      ttCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAAc                                                                                                                                 <  1:1329544/150‑1 (MQ=255)
                         ggCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCACAATAGCTTTTCAGGCGCTCGCGCAACTTc                                                                                                                              >  2:4166194/1‑150 (MQ=255)
                          gCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTcg                                                                                                                             <  2:2451547/150‑1 (MQ=255)
                            ccTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTcgcc                                                                                                                           <  2:2121182/150‑1 (MQ=255)
                            ccTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTcgcc                                                                                                                           <  1:2395994/150‑1 (MQ=255)
                              ttCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTcgccgc                                                                                                                         >  1:775676/1‑150 (MQ=255)
                                 tNCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAg                                                                                                                      >  1:326038/1‑150 (MQ=255)
                                  gCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGc                                                                                                                     <  1:1173258/150‑1 (MQ=255)
                                    cGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCgg                                                                                                                   <  2:1907369/150‑1 (MQ=255)
                                      tGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGt                                                                                                                 <  1:1451308/150‑1 (MQ=255)
                                         cTTTCATCGCACCCACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCt                                                                                                              <  1:2956775/150‑1 (MQ=255)
                                         cTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCt                                                                                                              <  1:4328391/150‑1 (MQ=255)
                                          gttCCTCGCACACACCGGACCGTCTGACTGCGGGTACACCAGCGCGTCAGTAGAAAACCCGCAATCCCCGCCGAGCTTACCTTTACTGTTTTTCCCTCCGCCAACATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTg                                                                                                             <  1:1253650/149‑1 (MQ=255)
                                            tCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGaa                                                                                                           >  2:4392431/1‑150 (MQ=255)
                                              aTCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATc                                                                                                         <  1:1959493/150‑1 (MQ=255)
                                               tCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCg                                                                                                        >  1:3141992/1‑150 (MQ=255)
                                               tCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCg                                                                                                        >  1:4567801/1‑150 (MQ=255)
                                                cGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGt                                                                                                       >  1:2804585/1‑150 (MQ=255)
                                                    cacacCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTc                                                                                                   <  2:3141992/150‑1 (MQ=255)
                                                       accgtaccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGcc                                                                                                >  2:4251594/1‑150 (MQ=255)
                                                          gtaccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGa                                                                                             >  2:1140455/1‑150 (MQ=255)
                                                           taccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGAt                                                                                            <  2:188488/150‑1 (MQ=255)
                                                           taccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGAt                                                                                            <  1:1129402/150‑1 (MQ=255)
                                                            accgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATc                                                                                           <  1:5079436/150‑1 (MQ=255)
                                                            accgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATc                                                                                           <  1:769609/150‑1 (MQ=255)
                                                             ccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCg                                                                                          >  1:4152443/1‑150 (MQ=255)
                                                             ccgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCCGCGGGTCCTGAATCGTTTTCGCCAGATCg                                                                                          >  2:1479017/1‑150 (MQ=255)
                                                              cgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCgg                                                                                         >  2:2371762/1‑150 (MQ=255)
                                                              cgtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCgg                                                                                         <  2:2253685/150‑1 (MQ=255)
                                                               gtCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGt                                                                                        >  1:5026067/1‑150 (MQ=255)
                                                                 cTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTCACTGTTTTTCACTCCGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGTTaa                                                                                      <  2:3607130/150‑1 (MQ=255)
                                                                  tgccTGCGGGTTCACCAGCGCTTCAGTAGAAAACACGCAATCCCCGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAACATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAAc                                                                                     <  2:2429578/147‑1 (MQ=255)
                                                                       gCGGGTACACCAGCGCTTCAGTAGAAATCACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAACTCGCTTTTCTGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTTAATCGTTTTCGCCAGATCGGTAACCTGaa                                                                                <  1:4125783/150‑1 (MQ=255)
                                                                         gggTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATc                                                                              <  2:5104484/150‑1 (MQ=255)
                                                                         gggTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATc                                                                              >  1:4897847/1‑150 (MQ=255)
                                                                          ggTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCt                                                                             >  1:31993/1‑150 (MQ=255)
                                                                             acacCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCg                                                                          >  1:268148/1‑150 (MQ=255)
                                                                             acacCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCg                                                                          <  2:381934/150‑1 (MQ=255)
                                                                               acCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTc                                                                        <  1:2620890/150‑1 (MQ=255)
                                                                                 cAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCgg                                                                      >  1:3547319/1‑150 (MQ=255)
                                                                                    ctcTTCAGTAGAAAACACTCAATCCACGACGAGCTTACCTTTACTGTTTTTCAATACGCAAAAATCGATTTTCAGGCGCTCGCGCAACTTCGACGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCTGTAACATTAATCTGCGTAGGAta                                                                   <  2:2224573/148‑2 (MQ=255)
                                                                                    cgcTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCg                                                                       <  2:3292288/146‑1 (MQ=255)
                                                                                    cgcTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATc                                                                   <  2:2008609/150‑1 (MQ=255)
                                                                                       ttCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAAt                                                                >  1:3913516/1‑150 (MQ=255)
                                                                                            tAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGcc                                                           >  2:2722777/1‑150 (MQ=255)
                                                                                               aaaaCACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCAc                                                        >  1:1836403/1‑150 (MQ=255)
                                                                                                aaaCACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCAcc                                                       >  2:4861936/1‑150 (MQ=255)
                                                                                                aaaCACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCAcc                                                       <  1:4937630/150‑1 (MQ=255)
                                                                                                  acacGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCcg                                                     >  1:2414975/1‑150 (MQ=255)
                                                                                                    acGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCcgcg                                                   >  2:5103491/1‑150 (MQ=255)
                                                                                                     cGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCgcgc                                                  >  1:1484449/1‑150 (MQ=255)
                                                                                                      gCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGcc                                                 >  1:3465934/1‑150 (MQ=255)
                                                                                                         aTCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCAcc                                              <  2:800500/150‑1 (MQ=255)
                                                                                                          tCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCg                                             <  1:2531024/150‑1 (MQ=255)
                                                                                                          tCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATAGTTTTCGCCAGATCGGTACCCTGAATCTTCGTAGGATCAATCTGCCCACCCGCGCCACCg                                             >  2:3979282/1‑150 (MQ=255)
                                                                                                          tCCACGCCGAGCTTACCTCTACTGATTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGACTCTGCGTCGGATCAATCTGCCCACCCGCGCCACCg                                             <  2:67339/150‑1 (MQ=255)
                                                                                                          tCCACGCCGAGCGTACCTGGACTGCTTTTCACTACGCCAAACTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGATTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCg                                             <  2:67432/150‑1 (MQ=255)
                                                                                                                ccGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTtgtg                                       <  2:2551619/150‑1 (MQ=255)
                                                                                                                ccGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTtgtg                                       >  1:2042335/1‑150 (MQ=255)
                                                                                                                ccGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGATTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTtgtg                                       >  2:2546559/1‑150 (MQ=255)
                                                                                                                 cGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGa                                      <  1:578531/150‑1 (MQ=255)
                                                                                                                 cGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGa                                      >  2:4456445/1‑150 (MQ=255)
                                                                                                                    gCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCa                                   >  1:2837392/1‑150 (MQ=255)
                                                                                                                       tACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTg                                <  1:2720600/150‑1 (MQ=255)
                                                                                                                         ccTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTggg                              >  2:5047079/1‑150 (MQ=255)
                                                                                                                          cTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGa                             <  1:3716983/150‑1 (MQ=255)
                                                                                                                          cTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGa                             <  1:4360537/150‑1 (MQ=255)
                                                                                                                           cttACTGTTTTTCACTACGCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCGGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGAt                            <  1:3585964/149‑1 (MQ=255)
                                                                                                                             tACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCGCACCCGCGCCACCGGTTGTGACCAGTGGGAttt                          <  2:3401147/150‑1 (MQ=255)
                                                                                                                             tACTGGTTTTCACTACGCCAACATCGCTTGTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGAttt                          <  2:2775135/150‑1 (MQ=255)
                                                                                                                                 gTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTAtt                      >  1:1902136/1‑150 (MQ=255)
                                                                                                                                   ttttCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTgc                    <  1:3837739/150‑1 (MQ=255)
                                                                                                                                   ttttCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTgc                    >  2:3028585/1‑150 (MQ=255)
                                                                                                                                     ttCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTgcgc                  <  1:4003165/150‑1 (MQ=255)
                                                                                                                                         cTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGAc              >  1:4287756/1‑150 (MQ=255)
                                                                                                                                         cTACGCCAAAAGAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCGGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGAc              <  1:2199919/150‑1 (MQ=255)
                                                                                                                                          tACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACa             >  1:2256995/1‑150 (MQ=255)
                                                                                                                                          tACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTGTATTGCGCCGACa             >  1:1032673/1‑150 (MQ=255)
                                                                                                                                            cGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAAt           <  1:3965484/150‑1 (MQ=255)
                                                                                                                                             gCCAAAATAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATa          <  2:4580889/150‑1 (MQ=255)
                                                                                                                                              ccAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATaa         <  2:1333540/150‑1 (MQ=255)
                                                                                                                                               cAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAg        >  1:4847557/1‑150 (MQ=255)
                                                                                                                                                 aaaTAGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCa      <  2:3473620/150‑1 (MQ=255)
                                                                                                                                                  aaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCaa     <  1:1236369/150‑1 (MQ=255)
                                                                                                                                                  aaTCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCaa     >  2:1691817/1‑150 (MQ=255)
                                                                                                                                                  aaTCGCTTTTCAGGCGCTCGCGAAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCaa     <  1:667415/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTTTTTCACTACGCCAAAATCGCTTTTCAGGCGCTCGCGCAACTTCGCCGCCAGCGGGTCCTGAATCGTTTTCGCCAGATCGGTAACCTGAATCTGCGTCGGATCAATCTGCCCACCCGCGCCACCGGTTGTGACCAGTGGGATTTTATTGCGCCGACAATAAGCAATCA  >  CP009273/2938491‑2938789

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