Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,189,448 C→T Q1367* (CAG→TAG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,189,4480CT99.3% 485.8 / ‑4.6 140Q1367* (CAG→TAG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (67/71);  total (68/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 148 total aligned reads displayed.

AGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCA  >  NC_000913/4189351‑4189546
                                                                                                 |                                                                                                  
 gcagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAgg                                                                                                 >  2:201373/1‑100 (MQ=255)
  cagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGa                                                                                                >  2:1100215/1‑100 (MQ=255)
    gcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATc                                                                                              >  1:938759/1‑100 (MQ=255)
    gcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATc                                                                                              >  1:2210637/1‑100 (MQ=255)
       gTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTa                                                                                           >  2:2375737/1‑100 (MQ=255)
       gTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTa                                                                                           >  2:1664848/1‑100 (MQ=255)
         tGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATg                                                                                         <  1:768233/100‑1 (MQ=255)
             ggCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTc                                                                                     >  1:721566/1‑100 (MQ=255)
               cAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGc                                                                                   <  2:972708/100‑1 (MQ=255)
               cAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGc                                                                                   >  1:67774/1‑100 (MQ=255)
                 aaCGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCg                                                                                 >  2:580283/1‑100 (MQ=255)
                 aaCGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCg                                                                                 <  2:1822880/100‑1 (MQ=255)
                 aaCGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCg                                                                                 <  1:353423/100‑1 (MQ=255)
                   cgcgACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGtg                                                                               <  2:2370426/100‑1 (MQ=255)
                     cgacgaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGtgct                                                                             >  2:2775159/1‑100 (MQ=255)
                          aaCTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCggg                                                                        >  2:476227/1‑100 (MQ=255)
                           aCTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGt                                                                       <  2:1604624/100‑1 (MQ=255)
                            cTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTg                                                                      <  2:1321074/100‑1 (MQ=255)
                             tGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCTGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGa                                                                     <  1:2105137/100‑1 (MQ=255)
                             tGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGa                                                                     <  2:1549752/100‑1 (MQ=255)
                              gcgcGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGaa                                                                    <  1:1587279/100‑1 (MQ=255)
                              gcgcGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGaa                                                                    <  2:721568/100‑1 (MQ=255)
                              gcgcGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGaa                                                                    <  2:402375/100‑1 (MQ=255)
                                gcgGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGc                                                                  <  1:2510494/100‑1 (MQ=255)
                                 cgGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCt                                                                 <  2:2080395/100‑1 (MQ=255)
                                 cgGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCt                                                                 >  1:1071643/1‑100 (MQ=255)
                                 cgGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGAGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCt                                                                 <  2:1658153/100‑1 (MQ=255)
                                  ggCCTGAAAGAGAACGTGATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGGGCTGCGGGTGAAGCTc                                                                <  1:273962/100‑1 (MQ=255)
                                   gCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTcc                                                               <  2:326153/100‑1 (MQ=255)
                                     cTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACGAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTcccg                                                             >  2:631543/1‑98 (MQ=255)
                                      tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGTTCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGc                                                            >  2:438816/1‑100 (MQ=255)
                                      tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGc                                                            <  1:201372/100‑1 (MQ=255)
                                      tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGc                                                            <  2:909689/100‑1 (MQ=255)
                                      tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGc                                                            <  1:2710584/100‑1 (MQ=255)
                                       gAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCt                                                           <  1:2151572/100‑1 (MQ=255)
                                         aaGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGc                                                         >  1:1269880/1‑100 (MQ=255)
                                           gagaACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCAc                                                       <  1:2640214/100‑1 (MQ=255)
                                           gagaACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCAc                                                       >  2:741848/1‑100 (MQ=255)
                                           gagaACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCAc                                                       <  1:1070188/100‑1 (MQ=255)
                                             gaACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCg                                                     >  1:1702516/1‑100 (MQ=255)
                                             gaACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCg                                                     >  2:1989359/1‑100 (MQ=255)
                                               aCGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCa                                                   <  1:1314050/100‑1 (MQ=255)
                                               aCGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGGGCTGCGGGTGAAGCTCCGGCTGCACCGCa                                                   >  1:1516534/1‑100 (MQ=255)
                                                cGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAg                                                  <  2:2589538/100‑1 (MQ=255)
                                                cGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAg                                                  >  1:820668/1‑100 (MQ=255)
                                                 gTTGTCGTGGGTCGTCTGATCCCGGGAGGTACCGGTTGCGCGTGCGACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTGCGGCTGCACCGCa                                                   <  1:2504771/98‑1 (MQ=255)
                                                 gTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAgg                                                 <  1:2880808/100‑1 (MQ=255)
                                                    aTCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGCTCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGa                                              >  2:2741930/1‑100 (MQ=255)
                                                    aTCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGa                                              >  1:949195/1‑100 (MQ=255)
                                                     tCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGAc                                             <  1:1485706/100‑1 (MQ=255)
                                                        tGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGc                                          <  2:2231564/100‑1 (MQ=255)
                                                        tGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGc                                          <  1:424205/100‑1 (MQ=255)
                                                         gggTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCa                                         >  2:476819/1‑100 (MQ=255)
                                                          ggTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCag                                        >  1:516842/1‑100 (MQ=255)
                                                           gtcgtcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCaga                                       <  2:32623/100‑1 (MQ=255)
                                                              gtcTGATCCCGGCGGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTTCGGCTGCACCGCAGGTGACTGCagaaga                                    >  2:1390276/1‑100 (MQ=255)
                                                              gtcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATGGTATGCGTCGCCGTGCTGCAGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCagaaga                                    >  1:1636668/1‑100 (MQ=255)
                                                              gtcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCagaaga                                    >  1:1208596/1‑100 (MQ=255)
                                                              gtcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCagaaga                                    >  1:2069444/1‑100 (MQ=255)
                                                               tcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGAc                                   >  1:1250033/1‑100 (MQ=255)
                                                               tcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGAc                                   >  1:429319/1‑100 (MQ=255)
                                                               tcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGAc                                   >  2:1659877/1‑100 (MQ=255)
                                                                 tGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGc                                 >  1:598383/1‑100 (MQ=255)
                                                                  gATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCa                                <  1:1685003/100‑1 (MQ=255)
                                                                  gATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCa                                >  2:2237349/1‑100 (MQ=255)
                                                                  gATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCa                                >  1:2221174/1‑100 (MQ=255)
                                                                   aTCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCAt                               >  2:135215/1‑100 (MQ=255)
                                                                      ccGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCtn                            <  2:778717/100‑2 (MQ=255)
                                                                        ggCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGcc                          <  2:157337/100‑1 (MQ=255)
                                                                            ggTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGcc                      >  2:785840/1‑100 (MQ=255)
                                                                            ggTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGcc                      >  1:697141/1‑100 (MQ=255)
                                                                            ggTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGcc                      >  2:1093410/1‑100 (MQ=255)
                                                                            ggTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGcc                      >  1:1957905/1‑100 (MQ=255)
                                                                             gTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGGCt                     <  2:236754/100‑1 (MQ=255)
                                                                             gTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCt                     >  2:1696992/1‑100 (MQ=255)
                                                                               aCCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACCGCAGAAGACGCATCTGCCAGCCTgg                   <  1:476817/100‑1 (MQ=255)
                                                                                ccGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGc                  <  1:2340476/100‑1 (MQ=255)
                                                                                ccGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGc                  <  1:2367093/100‑1 (MQ=255)
                                                                                ccGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGc                  >  1:2917055/1‑100 (MQ=255)
                                                                                 cGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGTTGAAGTTGCGGCTGCACCGCAGGTGACTGTAGAAGACGCATCTGCCAGCCTg                    <  1:1874247/97‑1 (MQ=255)
                                                                                 cGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCa                 >  2:1611338/1‑100 (MQ=255)
                                                                                  ggTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAg                >  2:1053744/1‑100 (MQ=255)
                                                                                  ggTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAg                <  1:416825/100‑1 (MQ=255)
                                                                                  ggTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAg                >  2:1806484/1‑100 (MQ=255)
                                                                                  ggTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAg                >  2:2139814/1‑100 (MQ=255)
                                                                                     tACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGGCTGGCAGatg             <  2:287052/100‑3 (MQ=255)
                                                                                     tACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAAc             <  1:1664256/100‑1 (MQ=255)
                                                                                     tACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAAc             <  2:2419190/100‑1 (MQ=255)
                                                                                      aCGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAAct            <  2:938763/100‑1 (MQ=255)
                                                                                      aCGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAAct            >  1:1735720/1‑100 (MQ=255)
                                                                                       cgcgTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAActg           <  1:1100210/100‑1 (MQ=255)
                                                                                        gcgTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAActgc          >  2:2462357/1‑100 (MQ=255)
                                                                                        gcgTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAActgc          <  1:1605094/100‑1 (MQ=255)
                                                                                         cgTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGACctgcg         >  2:891853/1‑99 (MQ=255)
                                                                                          gTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAActgctn        <  1:1416649/100‑2 (MQ=255)
                                                                                           tACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGa       <  1:476225/100‑1 (MQ=255)
                                                                                            accactAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGaa      <  1:51836/100‑1 (MQ=255)
                                                                                            accactAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGaa      >  2:1371289/1‑100 (MQ=255)
                                                                                             ccactaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGatcggaagagc     >  1:2643674/1‑90 (MQ=255)
                                                                                             ccactaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAAc     >  1:1200610/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
AGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCA  >  NC_000913/4189351‑4189546

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