Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,186,296 T→G 19.2% I316S (ATC→AGC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,186,2960TG19.2% 406.5 / 103.9 187I316S (ATC→AGC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (69/82);  new base G (17/19);  total (86/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

Only 100 of 190 total aligned reads displayed.

TGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTGACTA  >  NC_000913/4186197‑4186395
                                                                                                   |                                                                                                   
 ggATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATc                                                                                                    <  1:1326976/100‑1 (MQ=255)
 ggATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATc                                                                                                    >  2:2262347/1‑100 (MQ=255)
 ggATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATc                                                                                                    >  2:3570551/1‑100 (MQ=255)
    tCTGGCTGCGCCGGACATTATCGTACGTAACGAAAAACGTATGCGGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCAcc                                                                                                 <  2:1632478/100‑1 (MQ=255)
    tCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCAcc                                                                                                 <  1:3853216/100‑1 (MQ=255)
    tCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCAcc                                                                                                 <  2:1266143/100‑1 (MQ=255)
       ggCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGt                                                                                              >  2:1041439/1‑100 (MQ=255)
       ggCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGt                                                                                              >  2:20309/1‑100 (MQ=255)
       ggCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGt                                                                                              >  2:337415/1‑100 (MQ=255)
         cTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTc                                                                                            >  1:1596773/1‑100 (MQ=255)
           gcgcCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTa                                                                                          >  2:2684905/1‑100 (MQ=255)
            cgcCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGATCTaa                                                                                         >  1:877902/1‑100 (MQ=255)
            cgcCGGACAGCAGCGTACTGAACGAAAGGCGGATGCGGCAGGAGGCCGTATACGCCCTGCTGGATAACGGTCTTCGCGGTCGTGCGATCACCGGTTCTa                                                                                          <  2:2487359/99‑1 (MQ=255)
             gcCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAAc                                                                                        <  2:445410/100‑1 (MQ=255)
              ccGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACa                                                                                       <  1:64901/100‑1 (MQ=255)
                 atcttcatcGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAgn                                                                                    <  1:3930038/96‑2 (MQ=255)
                  aCATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCg                                                                                   <  2:3807087/100‑1 (MQ=255)
                   catcatCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGt                                                                                  >  1:3765754/1‑100 (MQ=255)
                   catcatCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCAGTCGTGCGATCACCGGTTCTAACAAGCg                                                                                   >  2:2104154/1‑99 (MQ=255)
                    atcatcGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTc                                                                                 <  2:842028/100‑1 (MQ=255)
                    atcatcGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTc                                                                                 <  2:843855/100‑1 (MQ=255)
                       atcGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGGCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCctc                                                                              <  1:2741596/100‑1 (MQ=255)
                       atcGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCctc                                                                              >  1:2344724/1‑100 (MQ=255)
                        tcGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCctct                                                                             <  2:701872/100‑1 (MQ=255)
                         cgtacgtaACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGGTCTAACAAGCGTCCTCTg                                                                            >  1:3740059/1‑100 (MQ=255)
                         cgtacctaACGAAAAACGTATACTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGGACCGGATCTAAAAAGCGGCctct                                                                             >  1:1731786/1‑99 (MQ=255)
                           tacgtaACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGaa                                                                          <  1:1224085/100‑1 (MQ=255)
                            acgtaACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGaaa                                                                         <  2:465763/100‑1 (MQ=255)
                            acgtaACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGaaa                                                                         >  1:562721/1‑100 (MQ=255)
                              gtaACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATc                                                                       <  2:2726797/100‑1 (MQ=255)
                              gtaACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATc                                                                       >  2:2323421/1‑100 (MQ=255)
                              gtaACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATc                                                                       >  2:806137/1‑100 (MQ=255)
                                aaCGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCtt                                                                     >  2:1990368/1‑100 (MQ=255)
                                aaCGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGAAATCtt                                                                     <  1:245064/100‑1 (MQ=255)
                                 aCGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGCCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCttt                                                                    <  1:2307706/100‑1 (MQ=255)
                                  ctcaaaaCGTATGCTGCAGGAAGTAGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCTATCACCGGTTCTAACAAGCGTCCTCTGAAATc                                                                       <  1:2609701/93‑1 (MQ=255)
                                   gAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTgg                                                                  <  2:3766320/100‑1 (MQ=255)
                                   gAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTgg                                                                  >  1:106361/1‑100 (MQ=255)
                                        aCGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGan                                                             <  2:1536982/100‑2 (MQ=255)
                                        aCGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGAc                                                             <  1:431925/100‑1 (MQ=255)
                                         cGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCATCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACa                                                            <  1:2627355/100‑1 (MQ=255)
                                         cGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACa                                                            >  2:1024374/1‑100 (MQ=255)
                                          gTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACAt                                                           <  1:375148/100‑1 (MQ=255)
                                             tgctgcAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGAt                                                        <  2:2510291/100‑1 (MQ=255)
                                              gctgcAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATc                                                       >  2:481695/1‑100 (MQ=255)
                                              gctgcAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCAAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATc                                                       >  2:3006941/1‑100 (MQ=255)
                                               ctgcAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCa                                                      >  2:3858493/1‑100 (MQ=255)
                                                tgcAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCaa                                                     >  1:2204008/1‑100 (MQ=255)
                                                 gcAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCaaa                                                    >  2:3547119/1‑100 (MQ=255)
                                                 gcAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCaaa                                                    <  1:2298748/100‑1 (MQ=255)
                                                  cAGGAAGCGGTAGGCGGCCTGCTGGATAACGGTCGTCGAGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCaaa                                                    <  1:3401498/99‑1 (MQ=255)
                                                   aGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAgg                                                  <  1:742605/100‑1 (MQ=255)
                                                    ggAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGt                                                 <  1:1963206/100‑1 (MQ=255)
                                                     gAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTa                                                <  2:1175717/100‑1 (MQ=255)
                                                       aGCGGTAGGCGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTa                                                <  1:3651589/98‑1 (MQ=255)
                                                         cGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACa                                            >  2:1839706/1‑100 (MQ=255)
                                                         cGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACa                                            <  1:371750/100‑1 (MQ=255)
                                                          ggTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAg                                           >  1:1180655/1‑100 (MQ=255)
                                                          ggTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAg                                           >  2:2707589/1‑100 (MQ=255)
                                                              gACGCCCTGCTGGATATCGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTc                                       <  1:2500662/100‑1 (MQ=255)
                                                                 gCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGtt                                    <  1:18672/100‑1 (MQ=255)
                                                                  cccTGCTGGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGttt                                   <  2:1666686/100‑1 (MQ=255)
                                                                   ccTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTc                                  <  2:3165443/100‑1 (MQ=255)
                                                                   ccTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTc                                  <  1:3107575/100‑1 (MQ=255)
                                                                    ctgctgGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTcc                                 <  1:1913502/100‑1 (MQ=255)
                                                                    ctgctgGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTcc                                 <  1:696421/100‑1 (MQ=255)
                                                                     tgctgGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCg                                <  1:1035737/100‑1 (MQ=255)
                                                                     ngctgGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCg                                >  1:1640203/2‑100 (MQ=255)
                                                                      gctgGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGt                               <  1:951853/100‑1 (MQ=255)
                                                                       ctgGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTc                              <  2:2795500/100‑1 (MQ=255)
                                                                       ctgGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTc                              <  2:3898480/100‑1 (MQ=255)
                                                                       ctgGATAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTc                              <  1:4161166/100‑1 (MQ=255)
                                                                        tgGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTGTAACAAGCGTCCTCTGAAATTTTTGGCCGCCATGGTCAACGGTAAACAGGGTCGGTTCCGTCa                             >  1:330674/1‑100 (MQ=255)
                                                                        tgGAGAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGGCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCa                             <  1:2780165/100‑1 (MQ=255)
                                                                         ggATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAg                            >  2:418356/1‑100 (MQ=255)
                                                                          gATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGa                           <  2:3524786/100‑1 (MQ=255)
                                                                            tAACGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAAc                         <  2:4016107/100‑1 (MQ=255)
                                                                             aaCGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAAcc                        <  2:181838/100‑1 (MQ=255)
                                                                              aaggTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCt                       >  2:3957379/3‑100 (MQ=255)
                                                                              aCGGTCGTCGCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCt                       >  1:3482162/1‑100 (MQ=255)
                                                                               cGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTg                      <  1:641693/100‑1 (MQ=255)
                                                                                ggTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGc                     >  2:2845883/1‑100 (MQ=255)
                                                                                ggTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGc                     <  1:4143645/100‑1 (MQ=255)
                                                                                 ttcgtcgCGGTCGTGCGAGCACCGGTTCTAACAAGCGTCCTCTGGAATCTTTGGCCGACATGCTCAAAGGTAAACCGGGTCGTTTCCGTCAGAACCTGc                     >  1:3528014/2‑99 (MQ=255)
                                                                                    gtcgCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCgg                 <  2:4239601/100‑1 (MQ=255)
                                                                                    gtcgCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCgg                 >  1:2078552/1‑100 (MQ=255)
                                                                                     tcgCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGt                >  2:3620331/1‑100 (MQ=255)
                                                                                        cgGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAg             <  2:3693879/100‑1 (MQ=255)
                                                                                        cgGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAg             <  2:4160685/100‑1 (MQ=255)
                                                                                          gTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTaa              <  1:4048556/97‑1 (MQ=255)
                                                                                          gTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCg           <  2:1686673/100‑1 (MQ=255)
                                                                                            cGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCgtg         >  1:423142/1‑100 (MQ=255)
                                                                                            cGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCgtg         >  1:3701151/1‑100 (MQ=255)
                                                                                               gCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTg      <  2:29526/100‑1 (MQ=255)
                                                                                               gCGATCACCGGTTCTAACAAGCGTCCTCCGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTg      >  2:1356362/1‑100 (MQ=255)
                                                                                                cGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTg      >  1:1317718/1‑99 (MQ=255)
                                                                                                cGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTGa     >  2:3038091/1‑100 (MQ=255)
                                                                                                 gATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTGac    <  1:1500941/100‑1 (MQ=255)
                                                                                                 gATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTGac    <  1:3991307/100‑1 (MQ=255)
                                                                                                  aTCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGAGTGTTGact   >  1:3692348/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTGACTA  >  NC_000913/4186197‑4186395

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: