Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 212,726 T→G 5.4% S132R (AGT→AGG tilS → tRNA(Ile)‑lysidine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913212,7260TG 5.4% 353.1 / 10.2 112S132R (AGT→AGGtilStRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base T (52/54);  new base G (4/2);  total (56/56)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 112 total aligned reads displayed.

CGCAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCA  >  NC_000913/212578‑212872
                                                                                                                                                    |                                                                                                                                                  
cgcAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTg                                                                                                                                                   >  1:4352/1‑150 (MQ=255)
 gcAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTgg                                                                                                                                                  <  2:1408982/150‑1 (MQ=255)
  cAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGc                                                                                                                                                 <  2:1440235/150‑1 (MQ=255)
  cAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGc                                                                                                                                                 <  2:1977919/150‑1 (MQ=255)
   aagaagGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGcc                                                                                                                                                <  1:774710/150‑1 (MQ=255)
    agaagGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGccc                                                                                                                                               >  2:389681/1‑150 (MQ=255)
             gggCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGGGCCACTGCCgcgtt                                                                                                                                      <  2:1941111/150‑5 (MQ=255)
             gggCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGGGCCACTGCCgcgtt                                                                                                                                      <  1:968442/150‑5 (MQ=255)
                  ttGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCgg                                                                                                                                 <  1:292596/150‑1 (MQ=255)
                   tGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGc                                                                                                                                <  1:1934935/150‑1 (MQ=255)
                   tGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGc                                                                                                                                <  1:1354751/150‑1 (MQ=255)
                   tGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGc                                                                                                                                <  1:1444459/150‑1 (MQ=255)
                     aGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTa                                                                                                                              >  2:869606/1‑150 (MQ=255)
                       gCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTa                                                                                                                              <  1:1411125/148‑1 (MQ=255)
                       gCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTa                                                                                                                              >  2:1411125/1‑148 (MQ=255)
                         ccAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGc                                                                                                                          >  2:1877454/1‑150 (MQ=255)
                           aGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCgg                                                                                                                        <  2:118174/150‑1 (MQ=255)
                              cgcgGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAg                                                                                                                     <  2:1239301/150‑1 (MQ=255)
                                     ggcaCGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGa                                                                                                              >  2:16313/1‑150 (MQ=255)
                                      gcaCGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAg                                                                                                             <  1:436462/150‑1 (MQ=255)
                                        aCGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGtt                                                                                                           >  1:285650/1‑150 (MQ=255)
                                            tATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGcc                                                                                                       <  2:1825068/150‑1 (MQ=255)
                                             ttCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCg                                                                                                      <  2:331171/149‑1 (MQ=255)
                                               cAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCTTGCCGGGCCTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGa                                                                                                    <  1:1161695/150‑1 (MQ=255)
                                                 ggCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAAc                                                                                                  >  2:469250/1‑150 (MQ=255)
                                                          ccGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTAt                                                                                         >  1:1431708/1‑150 (MQ=255)
                                                               ccTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGcc                                                                                    <  2:209855/150‑1 (MQ=255)
                                                               ccTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGcc                                                                                    <  2:1111744/150‑1 (MQ=255)
                                                                cTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGccc                                                                                   >  1:1605718/1‑150 (MQ=255)
                                                                 ttgttgCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCg                                                                                  >  2:74557/1‑150 (MQ=255)
                                                                     tgCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCgg                                                                                              >  2:1055359/1‑134 (MQ=255)
                                                                     tgCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCgg                                                                                              <  1:1055359/134‑1 (MQ=255)
                                                                     tgCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGc                                                                              >  2:1353945/1‑150 (MQ=255)
                                                                      gCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCt                                                                             >  2:732470/1‑150 (MQ=255)
                                                                       cccGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTc                                                                            >  1:979275/1‑150 (MQ=255)
                                                                       cccGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTc                                                                            >  2:1729884/1‑150 (MQ=255)
                                                                          ggTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGcc                                                                         <  2:1136770/150‑1 (MQ=255)
                                                                            tGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTTGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGCCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCg                                                                       <  2:923294/150‑1 (MQ=255)
                                                                             gAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGc                                                                      >  2:1565111/1‑150 (MQ=255)
                                                                                  gCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCAcgcg                                                                 >  1:1052828/1‑150 (MQ=255)
                                                                                  gCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCAcgcg                                                                 >  2:359122/1‑150 (MQ=255)
                                                                                    tGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCggg                                                               <  1:389681/150‑1 (MQ=255)
                                                                                     ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggg                                                              >  2:1190786/1‑150 (MQ=255)
                                                                                     ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggg                                                              >  2:1146216/1‑150 (MQ=255)
                                                                                       tCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggggg                                                            >  1:653463/1‑150 (MQ=255)
                                                                                        cACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGa                                                           >  1:552536/1‑150 (MQ=255)
                                                                                          ccGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAAc                                                         <  2:806400/150‑1 (MQ=255)
                                                                                               cAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTg                                                    >  1:1214290/1‑150 (MQ=255)
                                                                                                 aCATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCa                                                  <  2:451672/150‑1 (MQ=255)
                                                                                                   aTCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGt                                                <  2:1322733/150‑1 (MQ=255)
                                                                                                   aTCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGt                                                >  1:1159652/1‑150 (MQ=255)
                                                                                                    tctcGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTg                                               >  1:1936289/1‑150 (MQ=255)
                                                                                                     ctcGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTgg                                              >  2:1332162/1‑150 (MQ=255)
                                                                                                     ctcGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTgg                                              <  1:1353945/150‑1 (MQ=255)
                                                                                                       cgacgaTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGgc                                            >  2:1177342/1‑150 (MQ=255)
                                                                                                       cgacgaTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGgc                                            >  1:1454237/1‑150 (MQ=255)
                                                                                                       cgacgaTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGgc                                            <  2:1472090/150‑1 (MQ=255)
                                                                                                          cgaTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGcgcg                                         <  2:1431708/150‑1 (MQ=255)
                                                                                                          cgaTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGcgcg                                         >  1:663904/1‑150 (MQ=255)
                                                                                                             tCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCa                                      >  1:1018248/1‑150 (MQ=255)
                                                                                                             tCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCa                                      >  1:1011477/1‑150 (MQ=255)
                                                                                                                  gtgAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATg                                 <  1:1474686/150‑1 (MQ=255)
                                                                                                                  gtgAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATg                                 >  2:1811221/1‑150 (MQ=255)
                                                                                                                  gtgAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATg                                 <  2:1649650/150‑1 (MQ=255)
                                                                                                                    gAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGTGGAACTGGTGCAGTGGGCGCGTCAGTATGAt                               <  2:1404617/150‑1 (MQ=255)
                                                                                                                     aaaCCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGAtt                              <  2:386543/150‑1 (MQ=255)
                                                                                                                        ccTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTAc                           >  2:803793/1‑150 (MQ=255)
                                                                                                                        ccTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTAc                           <  1:732470/150‑1 (MQ=255)
                                                                                                                         cTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACg                          <  1:546692/150‑1 (MQ=255)
                                                                                                                         cTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACg                          >  1:631168/1‑150 (MQ=255)
                                                                                                                          tttCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGc                         >  1:612785/1‑150 (MQ=255)
                                                                                                                          tttCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGc                         >  2:1802470/1‑150 (MQ=255)
                                                                                                                           ttCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCt                        <  2:237993/150‑1 (MQ=255)
                                                                                                                            tCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATg                                 <  2:94911/140‑1 (MQ=255)
                                                                                                                            tCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATg                                 >  1:94911/1‑140 (MQ=255)
                                                                                                                            tCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTg                       <  1:359122/150‑1 (MQ=255)
                                                                                                                            tCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTg                       <  1:74557/150‑1 (MQ=255)
                                                                                                                             ctgctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTgg                      <  1:16313/150‑1 (MQ=255)
                                                                                                                              tgctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGa                     >  2:494916/1‑150 (MQ=255)
                                                                                                                               gctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTg                       <  2:1538788/147‑1 (MQ=255)
                                                                                                                               gctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTg                       >  1:1538788/1‑147 (MQ=255)
                                                                                                                               gctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGAt                    <  2:1040546/150‑1 (MQ=255)
                                                                                                                                ctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGAtt                   >  1:792569/1‑150 (MQ=255)
                                                                                                                                 tgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTg                  <  2:1052828/150‑1 (MQ=255)
                                                                                                                                 tgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTg                  <  1:918147/150‑1 (MQ=255)
                                                                                                                                 tgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTg                  >  1:1537546/1‑150 (MQ=255)
                                                                                                                                 tgGCGCTAAAACGCGGCAGTGACCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTg                  <  2:1529749/150‑1 (MQ=255)
                                                                                                                                  ggCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGTTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGa                 >  2:1240545/1‑150 (MQ=255)
                                                                                                                                  ggCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCGCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGGTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGa                 <  1:869606/150‑1 (MQ=255)
                                                                                                                                   gcgcTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGaa                <  2:285650/150‑1 (MQ=255)
                                                                                                                                    cgcTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAg               <  1:1146216/150‑1 (MQ=255)
                                                                                                                                     gcTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGa              >  1:228792/1‑150 (MQ=255)
                                                                                                                                        aaaaCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGa           >  2:728185/1‑150 (MQ=255)
                                                                                                                                         aagcccgGCAGGGCCACTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGaa          >  2:1608355/6‑150 (MQ=255)
                                                                                                                                         aaaCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTg                                                                               <  1:283105/81‑1 (MQ=255)
                                                                                                                                         aaaCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGaa          <  2:1605718/150‑1 (MQ=255)
                                                                                                                                           gcccgGCAGGGCCACTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAg        >  1:1442511/4‑150 (MQ=255)
                                                                                                                                           gcccgGCAGGGCCACTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAg        >  2:1667038/4‑150 (MQ=255)
                                                                                                                                           aCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAg        <  1:1026003/150‑1 (MQ=255)
                                                                                                                                            cccgGCAGGGCCACTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGt       >  1:1682162/3‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGCAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCA  >  NC_000913/212578‑212872

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