Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,543 C→A P1100Q (CCG→CAG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5430CA87.5% 400.9 / 52.0 152P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (9/10);  new base A (73/60);  total (82/70)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

Only 100 of 153 total aligned reads displayed.

GGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAA  >  NC_000913/4184446‑4184642
                                                                                                 |                                                                                                   
ggTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGc                                                                                                   >  1:3028694/1‑100 (MQ=255)
  tcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTg                                                                                                 <  1:3809516/100‑1 (MQ=255)
   cgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                >  2:3881393/1‑100 (MQ=255)
   cgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                >  2:4312633/1‑100 (MQ=255)
    gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTggg                                                                                               <  1:59402/100‑1 (MQ=255)
     tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                              <  1:4067233/100‑1 (MQ=255)
     tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                              <  1:125045/100‑1 (MQ=255)
      cACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCg                                                                                             >  2:2415498/1‑100 (MQ=255)
           tAACAAGGGTGTAATTTCTTAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCg                                                                                        <  2:3867817/100‑1 (MQ=255)
           naaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCg                                                                                        >  1:3525537/2‑100 (MQ=255)
             aCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtc                                                                                      >  1:3000921/1‑100 (MQ=255)
             aCAAGGGTGTAATCTCTAAGATCAACCCCATCGAAGAAATGCCTCACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtc                                                                                      <  2:379987/100‑1 (MQ=255)
             aCAAGGGAGTAATATCTAAGATCAACCCGATCGAAGCTATGCCTTACGATGAAAACGTTACGCAGGTAGACCTCGTACTGACCCAGCTGTGCGTCCAGtc                                                                                      >  2:2790793/1‑100 (MQ=255)
              cAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtct                                                                                     >  1:4412376/1‑100 (MQ=255)
              cAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtct                                                                                     >  2:2250304/1‑100 (MQ=255)
                aGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCg                                                                                   >  2:4320760/1‑100 (MQ=255)
                aGGGTGTAATTTCTAAGATCAACCCGAGCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCAGCTGGGCGTACCGTCCCg                                                                                   >  2:1493198/1‑100 (MQ=255)
                 gggTGAAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGt                                                                                  >  2:3811689/1‑100 (MQ=255)
                  ggTGTAATTTATAAGACCAACCCGATCGAAGATATCCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGACGTACCGTCTCGTa                                                                                 <  2:2137257/100‑1 (MQ=255)
                   gtgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCTGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                <  2:1296325/100‑1 (MQ=255)
                   gtgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                >  1:1463419/1‑100 (MQ=255)
                    tgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATg                                                                               >  1:4366084/1‑100 (MQ=255)
                    tgtAATTGGTAAGATAAACCCGCTCGATGAGATGCCTAACGACAAAAACGGTACGCCGGTAGACGTCGAATCGAGCCAGCTGGACGTACCGTCCCGGATg                                                                               <  2:1779223/100‑1 (MQ=255)
                     gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGCAAACGCTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGa                                                                              >  2:2681869/1‑100 (MQ=255)
                     gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGa                                                                              <  2:1506654/100‑1 (MQ=255)
                        aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                           >  2:2244678/1‑100 (MQ=255)
                        aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                           <  2:2866725/100‑1 (MQ=255)
                        aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                           >  1:1532944/1‑100 (MQ=255)
                        aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                           >  1:666387/1‑100 (MQ=255)
                         tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAt                                                                          >  1:597147/1‑100 (MQ=255)
                           tCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCg                                                                        >  1:1928083/1‑100 (MQ=255)
                            cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCgg                                                                       >  2:3938166/1‑100 (MQ=255)
                            cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCgg                                                                       >  2:1667568/1‑100 (MQ=255)
                              aaGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                     <  1:301626/100‑1 (MQ=255)
                              aaGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                     >  2:2263282/1‑100 (MQ=255)
                               aGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                    <  2:1899811/100‑1 (MQ=255)
                               aGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                    <  2:439223/100‑1 (MQ=255)
                                gATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAg                                                                   >  1:1129557/1‑100 (MQ=255)
                                gATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAg                                                                   >  2:3552802/1‑100 (MQ=255)
                                 aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGa                                                                  <  1:4253021/100‑1 (MQ=255)
                                  tCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                 >  2:3067736/1‑100 (MQ=255)
                                   cAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                <  1:48278/100‑1 (MQ=255)
                                   cAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                >  2:3140778/1‑100 (MQ=255)
                                   cAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                >  2:4192362/1‑100 (MQ=255)
                                    aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATcc                                                               >  1:544530/1‑100 (MQ=255)
                                    aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATcc                                                               <  2:3486493/100‑1 (MQ=255)
                                    aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATcc                                                               >  1:2758237/1‑100 (MQ=255)
                                    aaCCCGATCGAAGATATGCCTGACGGTGAAAACGGTGCGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATcc                                                               <  2:4095693/100‑1 (MQ=255)
                                       ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCg                                                            <  1:1361184/100‑1 (MQ=255)
                                        cGATCGAAGATATGCCTTACGATGAAATCGGTACGCCGGTATACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGa                                                           <  2:334716/100‑1 (MQ=255)
                                         gATCGAAGACATGCCTTACGATGACAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaa                                                          <  2:958532/100‑1 (MQ=255)
                                           tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                        <  2:1726500/100‑1 (MQ=255)
                                           tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                        >  1:3163252/1‑100 (MQ=255)
                                           tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                        <  2:617688/100‑1 (MQ=255)
                                            cGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                       >  2:367350/1‑100 (MQ=255)
                                             gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAccc                                                      <  1:3374348/100‑1 (MQ=255)
                                             gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATAAACATCGGTCAGATCCTCGAAAccc                                                      >  2:2571381/1‑100 (MQ=255)
                                              aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                     >  1:981194/1‑100 (MQ=255)
                                                gATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                   >  2:822361/1‑100 (MQ=255)
                                                gATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                   >  2:821299/1‑100 (MQ=255)
                                                  tatGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTg                                                 >  2:2920391/1‑100 (MQ=255)
                                                   atGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                <  1:3609054/100‑1 (MQ=255)
                                                    tGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                               <  1:4051108/100‑1 (MQ=255)
                                                    tGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                               <  2:2110333/100‑1 (MQ=255)
                                                    tGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                               <  1:3296600/100‑1 (MQ=255)
                                                      ccTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTa                                             <  1:2381601/100‑1 (MQ=255)
                                                       cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAt                                            >  2:440200/1‑100 (MQ=255)
                                                       cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAt                                            <  1:3842764/100‑1 (MQ=255)
                                                         tACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGACTCGTATGTACATAGGTAAGATCCTCGAAACAAACCTGGGTATgg                                          >  2:719719/1‑100 (MQ=255)
                                                          aCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGc                                         <  1:2295724/100‑1 (MQ=255)
                                                           cGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCt                                        <  1:204999/100‑1 (MQ=255)
                                                           cGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCt                                        >  1:1319399/1‑100 (MQ=255)
                                                            gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTg                                       >  2:2462827/1‑100 (MQ=255)
                                                             aTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGc                                      >  1:2220052/1‑100 (MQ=255)
                                                                 aaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                  >  2:273919/1‑100 (MQ=255)
                                                                   aCGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAgg                                >  2:2038241/1‑100 (MQ=255)
                                                                   aCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAgg                                <  1:2052241/100‑1 (MQ=255)
                                                                   aCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAgg                                <  2:1891943/100‑1 (MQ=255)
                                                                     ggTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTa                              >  1:3136995/1‑100 (MQ=255)
                                                                       tACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATc                            >  2:1960734/1‑100 (MQ=255)
                                                                        aCGCCGGTAGACAGCGTACTGAACCAGCTGTGCGTACCGTGGCGTATGAACATCGGTCAGATCCTCTAAACCCACCTGGGTATGGCTGCGAAATGTATCg                           <  2:1923235/100‑1 (MQ=255)
                                                                         cGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                          >  2:4076181/1‑100 (MQ=255)
                                                                          gCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGc                         <  1:1357885/100‑1 (MQ=255)
                                                                            cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                       >  1:1095211/1‑100 (MQ=255)
                                                                             ggTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                      <  2:2052365/100‑1 (MQ=255)
                                                                                   cATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGATCATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATTGGCGACAAGATc                <  2:4174443/100‑1 (MQ=255)
                                                                                   cATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATc                >  1:2537555/1‑100 (MQ=255)
                                                                                   cATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATc                >  2:1786820/1‑100 (MQ=255)
                                                                                    aTCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCa               >  2:3347121/1‑100 (MQ=255)
                                                                                    aTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCa               >  1:2762869/1‑100 (MQ=255)
                                                                                     tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCaa              <  2:3463866/100‑1 (MQ=255)
                                                                                      cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc             >  1:2372981/1‑100 (MQ=255)
                                                                                        tACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACCAGATCAACGc           >  1:1988714/1‑100 (MQ=255)
                                                                                         aCTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc          <  1:837182/100‑1 (MQ=255)
                                                                                          cTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCa         >  2:1174633/1‑100 (MQ=255)
                                                                                           tGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCAt        <  2:1587798/100‑1 (MQ=255)
                                                                                            gAACCAGCCGGGCGTATCGTCTCGTATGAACATCTTTCAGATCCACGAAACCCACCTGTGGATGGCTCCGAAAGGAATCGGCGACAAGATCAACGCCATg       >  2:1070414/1‑100 (MQ=255)
                                                                                             aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGc      <  1:2994810/100‑1 (MQ=255)
                                                                                              aCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCt     >  2:604370/1‑100 (MQ=255)
                                                                                                 cGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTTGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaa  >  2:678146/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAA  >  NC_000913/4184446‑4184642

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