Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,182,994 G→A A23T (GCT→ACT)  rpoC → RNA polymerase beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,182,9940GA99.2% 885.9 / NA 244A23T (GCT→ACT) rpoCRNA polymerase beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (131/110);  minor base C (0/2);  total (131/113)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

GCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAA  >  NC_000913/4182904‑4183083
                                                                                          |                                                                                         
gCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTAct                                                                                         <  1:828694/93‑1 (MQ=255)
 cTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTActc                                                                                        <  2:5702691/93‑1 (MQ=255)
  tAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGCTCAAAATTActct                                                                                       >  2:6286390/1‑93 (MQ=255)
  tAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTActct                                                                                       >  2:5810544/1‑93 (MQ=255)
  tAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTActct                                                                                       <  2:4332922/93‑1 (MQ=255)
  tAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTActct                                                                                       >  2:3597448/1‑93 (MQ=255)
  tAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTActct                                                                                       <  2:6248146/93‑1 (MQ=255)
   aaCTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTg                                                                                      >  2:2779961/1‑93 (MQ=255)
   aaCTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTg                                                                                      >  2:2263382/1‑93 (MQ=255)
   aaCTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTg                                                                                      >  2:1243303/1‑93 (MQ=255)
     cTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGc                                                                                    >  2:2952353/1‑93 (MQ=255)
      tCCGACGGGAGCAAATCCGTGACAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCt                                                                                   <  2:2511657/93‑1 (MQ=255)
      tCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCt                                                                                   >  1:4778674/1‑93 (MQ=255)
       ccGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCtt                                                                                  <  2:758288/93‑1 (MQ=255)
         gACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCg                                                                                >  1:408184/1‑93 (MQ=255)
         gACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCg                                                                                <  1:4668456/93‑1 (MQ=255)
          aCGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGc                                                                               >  2:791273/1‑93 (MQ=255)
          aCGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGc                                                                               <  1:3003702/93‑1 (MQ=255)
           cGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGcc                                                                              <  2:5475951/93‑1 (MQ=255)
            gggAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCa                                                                             >  1:6136305/1‑93 (MQ=255)
             ggAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAg                                                                            >  1:3678862/1‑93 (MQ=255)
             ggAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAg                                                                            >  1:2626624/1‑93 (MQ=255)
              gAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGa                                                                           >  2:3363908/1‑93 (MQ=255)
               aGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGAc                                                                          <  1:5170568/93‑1 (MQ=255)
                 cAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACAt                                                                        >  2:2373856/1‑93 (MQ=255)
                   aaTCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGa                                                                      <  2:2716263/93‑1 (MQ=255)
                   aaTCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGa                                                                      >  2:6431681/1‑93 (MQ=255)
                     tCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATc                                                                    >  2:613289/1‑93 (MQ=255)
                     tCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATc                                                                    <  1:4420045/93‑1 (MQ=255)
                       cGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCg                                                                  <  2:1910476/93‑1 (MQ=255)
                       cGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCg                                                                  >  1:6393166/1‑93 (MQ=255)
                        gTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGATTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATccct                                                                 >  2:1165036/1‑91 (MQ=255)
                         tGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGtt                                                                >  1:3156474/1‑93 (MQ=255)
                          gAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTc                                                               >  2:531887/1‑93 (MQ=255)
                          gAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTc                                                               >  1:72818/1‑93 (MQ=255)
                            aaGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCAt                                                             <  2:560870/93‑1 (MQ=255)
                            aaGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCAt                                                             <  2:593793/93‑1 (MQ=255)
                             aGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATg                                                            <  1:5691496/93‑1 (MQ=255)
                              gATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATgg                                                           >  1:3745099/1‑93 (MQ=255)
                               aTTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGt                                                          <  2:3040214/93‑1 (MQ=255)
                                  tattaAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCtt                                                       >  1:6201663/1‑93 (MQ=255)
                                  tattaAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCtt                                                       <  1:6374827/93‑1 (MQ=255)
                                   attaAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCttt                                                      <  2:4291551/93‑1 (MQ=255)
                                   attaAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCttt                                                      <  2:3141167/93‑1 (MQ=255)
                                     taAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCg                                                    <  1:4619047/93‑1 (MQ=255)
                                     taAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCg                                                    <  1:134363/93‑1 (MQ=255)
                                         gTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGa                                                >  1:1785694/1‑93 (MQ=255)
                                         gTTTCTGAAAGAGCAGACTAAAACCGAAGAGGTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGa                                                >  2:4046956/1‑93 (MQ=255)
                                          tttCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGaa                                               <  2:4764817/93‑1 (MQ=255)
                                          tttCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGaa                                               <  2:2224380/93‑1 (MQ=255)
                                            tCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGt                                             >  1:4000330/1‑93 (MQ=255)
                                             cTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGtt                                            <  2:3509675/93‑1 (MQ=255)
                                               gAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTaa                                          >  1:3515145/1‑93 (MQ=255)
                                               gAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTaa                                          >  1:2233986/1‑93 (MQ=255)
                                               gAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTaa                                          >  1:2021654/1‑93 (MQ=255)
                                                aaaGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTaaa                                         >  2:5372477/1‑93 (MQ=255)
                                                 aaGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTaaaa                                        <  2:4462482/93‑1 (MQ=255)
                                                    cgcAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGc                                     <  1:3950592/93‑1 (MQ=255)
                                                     gcAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGcc                                    <  1:5279450/93‑1 (MQ=255)
                                                     gcAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGcc                                    >  2:20299/1‑93 (MQ=255)
                                                      cAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCg                                   >  1:336908/1‑93 (MQ=255)
                                                       aGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCgg                                  >  2:4648684/1‑93 (MQ=255)
                                                       aGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCgg                                  >  2:296734/1‑93 (MQ=255)
                                                       aGACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCgg                                  >  1:2339552/1‑93 (MQ=255)
                                                        gACTAAAGCCGAAGAGTTTGATGCTATCAAACTTACTCTGGCTGCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGa                                 <  1:2564094/93‑1 (MQ=255)
                                                        gACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGa                                 >  1:2429007/1‑93 (MQ=255)
                                                        gACTAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGa                                 <  1:2775360/93‑1 (MQ=255)
                                                           tAAAACCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAAc                              <  1:5956254/93‑1 (MQ=255)
                                                           caaaaCCCGAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAAc                              <  2:1964159/92‑1 (MQ=255)
                                                             aaaCCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCa                            <  2:2322667/93‑1 (MQ=255)
                                                              aaCCGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCAt                           >  1:5549001/1‑93 (MQ=255)
                                                                ccGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCa                         >  2:5982494/1‑93 (MQ=255)
                                                                 cGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCaa                        >  2:3788577/1‑93 (MQ=255)
                                                                 cGAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCaa                        >  1:4223910/1‑93 (MQ=255)
                                                                  gAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAAc                       >  2:2292587/1‑93 (MQ=255)
                                                                  gAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAAc                       >  1:3308534/1‑93 (MQ=255)
                                                                  gAAGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAAc                       <  1:4335572/93‑1 (MQ=255)
                                                                   aaGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACt                      >  2:5433471/1‑93 (MQ=255)
                                                                   aaGAGTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACt                      >  2:4595742/1‑93 (MQ=255)
                                                                     gagTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTAc                    >  1:539276/1‑93 (MQ=255)
                                                                      agTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTAcc                   <  2:4905239/93‑1 (MQ=255)
                                                                       gTTTGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCg                  >  2:3909029/1‑93 (MQ=255)
                                                                        tttGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGt                 >  1:2742005/1‑93 (MQ=255)
                                                                        tttGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGt                 <  2:2211218/93‑1 (MQ=255)
                                                                         ttGATGCGATCAAAATTTCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTa                <  2:1944344/93‑1 (MQ=255)
                                                                         ttGATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTa                <  1:5870692/93‑1 (MQ=255)
                                                                         ttGATGCGATAAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTa                <  1:4906693/93‑1 (MQ=255)
                                                                           gATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACg              >  2:3820055/1‑93 (MQ=255)
                                                                           gATGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACg              >  2:3119021/1‑93 (MQ=255)
                                                                            aTGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGt             <  2:151116/93‑1 (MQ=255)
                                                                            aTGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGt             >  1:5370121/1‑93 (MQ=255)
                                                                             tGCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGtt            >  2:210171/1‑93 (MQ=255)
                                                                              gCGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTc           >  1:1083337/1‑93 (MQ=255)
                                                                              gCGATCAAAATTACTCGGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTc           <  2:5075667/93‑1 (MQ=255)
                                                                               cGATCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCa          >  2:1326860/1‑93 (MQ=255)
                                                                                 aTCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCaaa        >  2:2595043/1‑93 (MQ=255)
                                                                                 aTCAAAATTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCaaa        >  1:2126395/1‑93 (MQ=255)
                                                                                    aaaaTTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCa     >  1:5673993/1‑93 (MQ=255)
                                                                                      aaTTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGa   <  1:1040392/93‑1 (MQ=255)
                                                                                       aTTACTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGaa  <  2:5857874/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCTAACTCCGACGGGAGCAAATCCGTGAAAGATTTATTAAAGTTTCTGAAAGCGCAGACTAAAACCGAAGAGTTTGATGCGATCAAAATTGCTCTGGCTTCGCCAGACATGATCCGTTCATGGTCTTTCGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAA  >  NC_000913/4182904‑4183083

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