Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,711,397 C→T 79.7% R28Q (CGA→CAA)  tyrS ← tyrosyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,711,3970CT79.7% 784.6 / 206.1 345R28Q (CGA→CAA) tyrStyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (32/38);  new base T (129/145);  total (161/184)
Fisher's exact test for biased strand distribution p-value = 8.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

Only 100 of 345 total aligned reads displayed.

CCGTCGCGCCGCCTACCAGCGCAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAA  >  CP009273/1711251‑1711542
                                                                                                                                                  |                                                                                                                                                 
 cGTCGCGCCGCCTACCAGCGCAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGctc                                                                                                                                               >  1:2980856/1‑150 (MQ=255)
  gTCGCGCCGCCTACCAGCGCAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGctct                                                                                                                                              >  1:737338/1‑150 (MQ=255)
        ccgccTACCAGCGCAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAAc                                                                                                                                        <  2:1520606/150‑1 (MQ=255)
        ccgccTACCAGCGCAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAAc                                                                                                                                        >  1:2529439/1‑150 (MQ=255)
         cgccTACCAGCGCAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACg                                                                                                                                       >  1:1591554/1‑150 (MQ=255)
              aCCAGCGCAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTcc                                                                                                                                  >  2:1465523/1‑150 (MQ=255)
               ccAGCGCAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCt                                                                                                                                 >  2:5643052/1‑150 (MQ=255)
                cAGCGCAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTc                                                                                                                                >  1:4157614/1‑150 (MQ=255)
                   cgcAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTcgtc                                                                                                                             <  2:2198245/150‑1 (MQ=255)
                    gcAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTcgtcc                                                                                                                            >  2:5503942/1‑150 (MQ=255)
                     cAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTcgtccg                                                                                                                           >  2:78604/1‑150 (MQ=255)
                        ccGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCa                                                                                                                        <  1:4419793/150‑1 (MQ=255)
                         cGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCAc                                                                                                                       >  1:4495124/1‑150 (MQ=255)
                          ggCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGAGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCAcc                                                                                                                      <  2:5196261/150‑1 (MQ=255)
                           gCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCt                                                                                                                     <  1:581209/150‑1 (MQ=255)
                           gCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCt                                                                                                                     <  1:4870404/150‑1 (MQ=255)
                            cTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTg                                                                                                                    >  1:1020184/1‑150 (MQ=255)
                             ttGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTgg                                                                                                                   <  1:1331869/150‑1 (MQ=255)
                             ttGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTgg                                                                                                                   <  2:5191705/150‑1 (MQ=255)
                              tgtgGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTggg                                                                                                                  >  1:53856/1‑150 (MQ=255)
                               gtgGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGc                                                                                                                 >  2:1390870/1‑150 (MQ=255)
                                 ggCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTa                                                                                                               >  1:4923645/1‑150 (MQ=255)
                                 ggCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTa                                                                                                               <  1:5570283/150‑1 (MQ=255)
                                        ctgctgGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCATCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGccc                                                                                                        >  1:3912072/1‑150 (MQ=255)
                                        ctgctgGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGccc                                                                                                        >  2:3275377/1‑150 (MQ=255)
                                          gctgGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGccccc                                                                                                      <  2:1640674/150‑1 (MQ=255)
                                           ctgGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCg                                                                                                     <  1:4006176/150‑1 (MQ=255)
                                           ctgGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCg                                                                                                     <  1:2833653/150‑1 (MQ=255)
                                             ggAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGct                                                                                                   >  2:1091853/1‑150 (MQ=255)
                                                aGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCtt                                                                                                <  2:4844112/150‑1 (MQ=255)
                                                aGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCtt                                                                                                <  1:3136846/150‑1 (MQ=255)
                                                aGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCtt                                                                                                <  1:3450824/150‑1 (MQ=255)
                                                   gTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCa                                                                                             <  1:2898342/150‑1 (MQ=255)
                                                    tttCAGGCATAACAATGGAACAAGATGCCCCTAATGCATGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGGCGGTCACCTGGGCTACCAGCCCCCGCTCTTGCaa                                                                                            <  2:3683819/150‑1 (MQ=255)
                                                    tttCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCaa                                                                                            <  1:2194470/150‑1 (MQ=255)
                                                      tCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAAtt                                                                                          >  2:707896/1‑150 (MQ=255)
                                                         ggCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGtt                                                                                       <  1:1368195/150‑1 (MQ=255)
                                                         ggCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGtt                                                                                       <  1:2649784/150‑1 (MQ=255)
                                                           cATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTa                                                                                     >  2:2504629/1‑150 (MQ=255)
                                                            aTAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTaa                                                                                    <  1:2584421/150‑1 (MQ=255)
                                                             tAACCATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAAt                                                                                   >  2:4381113/1‑150 (MQ=255)
                                                              aaCAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATc                                                                                  >  1:1882112/1‑150 (MQ=255)
                                                               aCAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCa                                                                                 <  1:2126629/150‑1 (MQ=255)
                                                               aCAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCa                                                                                 >  2:4090398/1‑150 (MQ=255)
                                                                cAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCaa                                                                                >  2:2132649/1‑150 (MQ=255)
                                                                cAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCaa                                                                                >  2:4417125/1‑150 (MQ=255)
                                                                 aaTGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAg                                                                               >  2:3600554/1‑150 (MQ=255)
                                                                    ggAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTa                                                                            >  1:1686397/1‑150 (MQ=255)
                                                                     gAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTAc                                                                           >  2:4272240/1‑150 (MQ=255)
                                                                      aaCAAGAGGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACt                                                                          <  2:2792286/150‑1 (MQ=255)
                                                                         aaGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCt                                                                       >  2:3149121/1‑150 (MQ=255)
                                                                         aaGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCt                                                                       <  2:536746/150‑1 (MQ=255)
                                                                          aGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCtt                                                                      <  1:4288356/150‑1 (MQ=255)
                                                                          aGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCtt                                                                      >  2:835071/1‑150 (MQ=255)
                                                                             tGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGcc                                                                   <  1:1129994/150‑1 (MQ=255)
                                                                              gCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCa                                                                  >  2:4177942/1‑150 (MQ=255)
                                                                               ccccAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGGCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTTATCAAGTTACTGCTTGCCAt                                                                 <  2:4972756/150‑1 (MQ=255)
                                                                               ccccAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCAt                                                                 >  2:1008963/1‑150 (MQ=255)
                                                                                 ccAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCGCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCa                                                               >  2:1831001/1‑150 (MQ=255)
                                                                                  cAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCaa                                                              <  1:4095025/150‑1 (MQ=255)
                                                                                   aaaTGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCaaa                                                             >  2:2585016/1‑150 (MQ=255)
                                                                                      tGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAAtc                                                          >  2:4650005/1‑150 (MQ=255)
                                                                                      tGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAAtc                                                          >  2:3874342/1‑150 (MQ=255)
                                                                                       gCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAAtct                                                         <  2:4613076/150‑1 (MQ=255)
                                                                                          aGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCa                                                      >  2:4072573/1‑150 (MQ=255)
                                                                                           gCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCAt                                                     <  1:3564675/150‑1 (MQ=255)
                                                                                              gTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGta                                                  >  2:1446282/1‑150 (MQ=255)
                                                                                               tCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATAGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGtat                                                 >  2:3002646/1‑150 (MQ=255)
                                                                                                 aGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGtatat                                               <  2:3069534/150‑1 (MQ=255)
                                                                                                 aGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGtatat                                               >  1:1567488/1‑150 (MQ=255)
                                                                                                  gCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGtatata                                              <  2:5518084/150‑1 (MQ=255)
                                                                                                    ggTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAAc                                            >  1:3430262/1‑150 (MQ=255)
                                                                                                    ggTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAAc                                            <  1:1206884/150‑1 (MQ=255)
                                                                                                    ggTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAAc                                            <  2:4495124/150‑1 (MQ=255)
                                                                                                       aGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGAc                                         <  1:3126726/150‑1 (MQ=255)
                                                                                                        ggATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACt                                        <  2:500102/150‑1 (MQ=255)
                                                                                                              aaGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCt                                  >  1:2789253/1‑150 (MQ=255)
                                                                                                              aaGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCt                                  >  1:2668910/1‑150 (MQ=255)
                                                                                                              aaGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCGTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCt                                  <  2:2736330/150‑1 (MQ=255)
                                                                                                              aaGCCGCAAAATAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACGAGCCCCCGCTCTTGCAAATGTTTAATCAAGTTACTGCTTGCCATCAAAATATCCATGTATATAACGACTGCAcat                                  >  1:4745896/1‑148 (MQ=255)
                                                                                                               aGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCtt                                 <  1:3218072/150‑1 (MQ=255)
                                                                                                               aGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCtt                                 <  2:4610380/150‑1 (MQ=255)
                                                                                                                gCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCttt                                >  2:2082474/1‑150 (MQ=255)
                                                                                                                gCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCttt                                <  1:3845976/150‑1 (MQ=255)
                                                                                                                   gCAATAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGcc                             <  1:4964191/150‑1 (MQ=255)
                                                                                                                     aaTAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCgg                           >  2:4700374/1‑150 (MQ=255)
                                                                                                                      aTAGAGCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGt                          >  2:4009860/1‑150 (MQ=255)
                                                                                                                          agCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTacac                      <  2:3912072/150‑1 (MQ=255)
                                                                                                                          agCGCGATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTacac                      <  2:3914756/150‑1 (MQ=255)
                                                                                                                              cgATCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACt                  <  1:4045059/150‑1 (MQ=255)
                                                                                                                               gATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACtt                 >  2:96896/1‑150 (MQ=255)
                                                                                                                                aTCGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTGCATGTATATAACGACTGCACCTTTGCCGGTACACGACttt                <  2:5051252/150‑1 (MQ=255)
                                                                                                                                aTCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACttt                <  1:781037/150‑1 (MQ=255)
                                                                                                                                  cGGGCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTc              <  2:1341720/150‑1 (MQ=255)
                                                                                                                                    ggCCTTGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGc            >  2:507603/1‑150 (MQ=255)
                                                                                                                                         tGCGCCAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGAt       <  1:4703781/150‑1 (MQ=255)
                                                                                                                                          gcgcCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATg      <  2:2110261/150‑1 (MQ=255)
                                                                                                                                             ccAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGa   <  1:196750/150‑1 (MQ=255)
                                                                                                                                              cAGTTGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGaa  >  2:899945/1‑150 (MQ=255)
                                                                                                                                              cAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGaa  <  2:53856/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCGTCGCGCCGCCTACCAGCGCAACCGGCTTGTGGCCCGCCTGCTGGAAGCGTTTCAGGCATAACAATGGAACAAGATGCCCCAAATGCAAGCTGTCAGCGGTAGGATCGAAGCCGCAATAGAGCGCGATCGGGCCTTGCGCCAGTCGCTCTGCTAACGCTTCCTCGTCCGTCACCTGGGCTACCAGCCCCCGCTCTTGCAATTGTTTAATCAAGTTACTGCTTGCCATCAAAATCTCCATGTATATAACGACTGCACCTTTGCCGGTACACGACTTTTCGCCAGATGCGAA  >  CP009273/1711251‑1711542

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: ‑