Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 212,728 G→C 5.5% G133A (GGC→GCC)  tilS → tRNA(Ile)‑lysidine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913212,7280GC 5.5% 577.7 / 15.5 172G133A (GGC→GCC) tilStRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base G (84/77);  new base C (3/7);  total (88/84)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

GCAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAGGACG  >  NC_000913/212579‑212877
                                                                                                                                                     |                                                                                                                                                     
gcAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTgg                                                                                                                                                       >  2:852396/1‑150 (MQ=255)
gcAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTgg                                                                                                                                                       >  2:561826/1‑150 (MQ=255)
 cAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGc                                                                                                                                                      >  1:1145819/1‑150 (MQ=255)
 cAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGc                                                                                                                                                      >  1:13556/1‑150 (MQ=255)
  aagaagGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGcc                                                                                                                                                     <  2:1624348/150‑1 (MQ=255)
      agGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGc                                                                                                                                                 >  2:746459/1‑150 (MQ=255)
       ggACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGcc                                                                                                                                                >  1:33237/1‑150 (MQ=255)
         aCTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCgg                                                                                                                                              <  1:1979433/150‑1 (MQ=255)
          cTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGGGCCACTGCCgcg                                                                                                                                             <  1:281563/150‑3 (MQ=255)
            gggCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCt                                                                                                                                           >  2:50536/1‑150 (MQ=255)
              gCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTGGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCttt                                                                                                                                         <  1:1693897/150‑1 (MQ=255)
               cATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTc                                                                                                                                        >  1:520504/1‑150 (MQ=255)
               cATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTc                                                                                                                                        >  2:1885182/1‑150 (MQ=255)
                aTTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCg                                                                                                                                       >  2:467392/1‑150 (MQ=255)
                aTTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCg                                                                                                                                       <  1:681376/150‑1 (MQ=255)
                aTTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCGCGACGATCAATATGAAACAATTATGCTGGCGCTAAAACGCGGCAGAGCCCCTGCCGGGCtggcg                                                                                                                                       >  2:1203163/1‑146 (MQ=255)
                 ttaaGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGGGCCACTGCCGCGttttagc                                                                                                                                      <  2:1541851/147‑4 (MQ=255)
                 ttGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCgg                                                                                                                                      >  1:1434193/1‑150 (MQ=255)
                 ttGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCgg                                                                                                                                      <  2:369615/150‑1 (MQ=255)
                 ttGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCGGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCgg                                                                                                                                      <  2:1271556/150‑1 (MQ=255)
                     ggCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAAAATCTCGACGATAAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGTAATTCGGCGAt                                                                                                                                  >  1:1052345/1‑150 (MQ=255)
                       cccAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCCCCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATgg                                                                                                                                >  2:696370/1‑150 (MQ=255)
                       cccAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTCCAGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGACGCGCTTTCGGCTATgg                                                                                                                                >  2:1551412/1‑150 (MQ=255)
                        ccAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGCAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGACCTGCCGGGATTTCGTCTATGGc                                                                                                                               >  1:1814382/1‑150 (MQ=255)
                        ccAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGc                                                                                                                               <  2:32101/150‑1 (MQ=255)
                         cAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCg                                                                                                                              >  1:2008056/1‑150 (MQ=255)
                         cAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCg                                                                                                                              >  2:1790460/1‑150 (MQ=255)
                          aGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGGGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCgg                                                                                                                             <  2:1463649/150‑1 (MQ=255)
                             cgcgGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAg                                                                                                                          <  2:1701854/150‑1 (MQ=255)
                                    ggcaCGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGa                                                                                                                   <  1:1687240/150‑1 (MQ=255)
                                     gcaCGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAg                                                                                                                  <  1:307766/150‑1 (MQ=255)
                                            aTCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGGTTTGGCGGAAGTCGCGGAGTTTGCCg                                                                                                           >  1:1870904/1‑150 (MQ=255)
                                              cAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGa                                                                                                         >  1:1489460/1‑150 (MQ=255)
                                              cAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGa                                                                                                         >  1:127703/1‑150 (MQ=255)
                                               aGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCGGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGaa                                                                                                        <  1:286090/150‑1 (MQ=255)
                                                   aTTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAAcgcg                                                                                                    >  1:603440/1‑150 (MQ=255)
                                                    tttGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCgg                                                                                                   >  2:432010/1‑150 (MQ=255)
                                                    tttGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGGCACCGCGCAACATCTCGACGATCAATGAGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCgg                                                                                                   <  2:1564394/150‑1 (MQ=255)
                                                       gCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCtt                                                                                                <  2:1228033/150‑1 (MQ=255)
                                                        cccGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCtga                                                                                               >  1:1929942/1‑148 (MQ=255)
                                                        cccGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTa                                                                                               >  1:887626/1‑150 (MQ=255)
                                                          cGCACCTTGTTGCCCGGTGAAGTGCTGGGCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTAtt                                                                                             <  1:935726/150‑1 (MQ=255)
                                                            cACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCg                                                                                           <  1:1273340/150‑1 (MQ=255)
                                                              ccTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGcc                                                                                         <  1:1193361/150‑1 (MQ=255)
                                                                  gttgCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGtt                                                                                     <  2:265637/150‑1 (MQ=255)
                                                                  gttgCCCGGTGAAGTGCTGGGCACCGCGCAACATCTCGACGCTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGtt                                                                                     <  2:2018347/150‑1 (MQ=255)
                                                                    tgCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAACCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGATCGCGGCTTATTCGCCCGTTGc                                                                                   >  1:636385/1‑150 (MQ=255)
                                                                    tgCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGc                                                                                   <  1:1885182/150‑1 (MQ=255)
                                                                     gCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCt                                                                                  >  1:1799086/1‑150 (MQ=255)
                                                                       ccGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCg                                                                                <  1:1231385/150‑1 (MQ=255)
                                                                       ccGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCg                                                                                <  1:1734067/150‑1 (MQ=255)
                                                                         ggTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGcc                                                                              <  2:1490695/150‑1 (MQ=255)
                                                                           tGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCg                                                                            <  2:423620/150‑1 (MQ=255)
                                                                           tGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCg                                                                            >  2:649224/1‑150 (MQ=255)
                                                                              aGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGTAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCAc                                                                         <  2:1076141/150‑1 (MQ=255)
                                                                                 gCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCAcgcg                                                                      >  1:1413336/1‑150 (MQ=255)
                                                                                 gCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCAcgcg                                                                      >  2:636774/1‑150 (MQ=255)
                                                                                 gCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGATAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCAcgcg                                                                      >  2:1955074/1‑150 (MQ=255)
                                                                                    ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggg                                                                   >  2:102483/1‑150 (MQ=255)
                                                                                     gTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCggggg                                                                  >  1:493513/1‑150 (MQ=255)
                                                                                      tCACCGCACAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggggg                                                                 >  1:1796244/1‑150 (MQ=255)
                                                                                      gcACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTAGGGCGGAAGGCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggggg                                                                 <  1:99778/149‑1 (MQ=255)
                                                                                       cACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGa                                                                >  2:282606/1‑150 (MQ=255)
                                                                                       cACCGCGCAACATCTCGACGATCAATGGGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGa                                                                <  1:868596/150‑1 (MQ=255)
                                                                                        aCCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCGGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCCCGCGGGGGGaa                                                               >  2:1201502/1‑150 (MQ=255)
                                                                                          cgcgCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACt                                                             >  1:390643/1‑150 (MQ=255)
                                                                                            cgcAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTgg                                                           <  2:1297856/150‑1 (MQ=255)
                                                                                            cgcAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTgg                                                           <  2:33237/150‑1 (MQ=255)
                                                                                            cgcAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCGTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTgg                                                           <  2:1434193/150‑1 (MQ=255)
                                                                                               aaCATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGc                                                        <  2:726849/150‑1 (MQ=255)
                                                                                                         cgaTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAGGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGcgcg                                              >  1:454020/1‑150 (MQ=255)
                                                                                                          gaTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGcgcgg                                             >  1:618514/1‑149 (MQ=255)
                                                                                                           aTCAATGTGAAACCTTGCTGCTGGCGCTACAACGCGGCAGGGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTc                                            <  1:588405/150‑1 (MQ=255)
                                                                                                            tCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCa                                           >  2:938524/1‑150 (MQ=255)
                                                                                                            tCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCa                                           >  1:424588/1‑150 (MQ=255)
                                                                                                             cAATGTGAACCCTTTCTGCTGGCGCTAACACGCGGCAGTGGCCCTGCCGGGCTTTCGGCAATGGCGCAAGGCTCGGAGATTGCCGGACCGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAg                                          <  1:55777/150‑1 (MQ=255)
                                                                                                               aTGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTa                                        <  1:942982/150‑1 (MQ=255)
                                                                                                                tgtgAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGTAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTAt                                       <  1:800875/150‑1 (MQ=255)
                                                                                                                tgtgAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTAt                                       <  1:917246/150‑1 (MQ=255)
                                                                                                                    aacccTATCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGTGAACTGGTGCAGTGGGCGCGTCAGTATGAtt                                   <  1:1055235/147‑1 (MQ=255)
                                                                                                                            ctgctgGCGCTAACACTCGGCAGTGGCCCTGCCGGGCTTTCGGCTAAGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGATGCTCGCCCGCACGCGGGTGTAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTgg                           <  1:1768844/150‑1 (MQ=255)
                                                                                                                              gctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAGGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGAt                         <  2:454020/150‑1 (MQ=255)
                                                                                                                               ctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGAtt                        <  2:493513/150‑1 (MQ=255)
                                                                                                                               ctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGAtt                        >  1:1936165/1‑150 (MQ=255)
                                                                                                                               cgggCGCAACACCGCGGCCGGGGCCCTGCCGGGCTTTCGGCTATGGCGGACGGCTCTGAGTGTGGCGGAACGCTCCTTATTCGCCCGTTGCTCTCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGAtt                        <  2:314140/148‑1 (MQ=255)
                                                                                                                                 ggCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGa                      >  2:737351/1‑150 (MQ=255)
                                                                                                                                  gcgcTAAACCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGaa                     <  1:1261963/150‑1 (MQ=255)
                                                                                                                                   cgcTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAg                    >  2:1594680/1‑150 (MQ=255)
                                                                                                                                   cgcTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCCGGATTGAAg                    <  2:2008056/150‑1 (MQ=255)
                                                                                                                                     cTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGTCTATGGCGGAAGTCTCTGAGTTTGCCGGTACGCGGCTTATTCGCCCGTTGCTCGCGCGCACGCGGGGGGAACTGGTGCAGTTGGCGCGTCAGTATGATTTACGCTGGATTGATGAc                  <  2:1052345/150‑1 (MQ=255)
                                                                                                                                     cTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGAc                  >  2:1296549/1‑150 (MQ=255)
                                                                                                                                      tAAACCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGTGGGAACTGGTGCAGTGGGCGCGTCATTATGATTTACGCTGGCTTGAAGACg                 >  1:1240382/1‑150 (MQ=255)
                                                                                                                                         aaCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGaaa              >  1:1511415/1‑150 (MQ=255)
                                                                                                                                          aCGCGGCAGTGGCCCTGCCGGGCGTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAg             <  2:1888298/150‑1 (MQ=255)
                                                                                                                                             cgGCAGTGGCCCTGCCGGGCTTTCGGCGATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGAGTGAAGACGAAAGTaa          <  1:1401245/150‑1 (MQ=255)
                                                                                                                                              ggCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAAt         >  1:1089554/1‑150 (MQ=255)
                                                                                                                                              ggCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAAt         <  2:1844203/150‑1 (MQ=255)
                                                                                                                                               gCAGTGGCCCGGCCGGGCTTTCGGCTATGGCGGAAGGCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACTCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTCCGCTGGATTGAAGACGAAAGTAATc        <  2:977278/150‑1 (MQ=255)
                                                                                                                                                    ggCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAGgac   >  2:490098/1‑150 (MQ=255)
                                                                                                                                                     gCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAGgact  >  2:241309/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAGGACG  >  NC_000913/212579‑212877

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