Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,418,629 T→C 16.4% intergenic (+18/‑413) yhdV → / → yhdW putative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,418,6290TC16.4% 567.2 / 105.9 267intergenic (+18/‑413)yhdV/yhdWputative outer membrane protein/pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (121/101);  new base C (28/15);  total (151/116)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-03

Only 100 of 270 total aligned reads displayed.

CTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGT  >  NC_000913/3418530‑3418727
                                                                                                   |                                                                                                  
cTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGt                                                                                                    <  1:4050902/100‑1 (MQ=255)
    gAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGCTCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTccc                                                                                                 >  2:3528165/1‑99 (MQ=255)
     aaaCGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCcgc                                                                                               <  2:3164460/100‑1 (MQ=255)
       aCGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCgcgc                                                                                             <  1:2607943/100‑1 (MQ=255)
       aCGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCgcgc                                                                                             <  2:654727/100‑1 (MQ=255)
          gCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGctct                                                                                          >  1:1234212/1‑100 (MQ=255)
          gCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                          <  1:2798454/100‑1 (MQ=255)
          gCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                          <  1:2579988/100‑1 (MQ=255)
          gCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                          <  1:3204999/100‑1 (MQ=255)
          gCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                          <  1:815352/100‑1 (MQ=255)
          gCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccct                                                                                          <  2:3689746/100‑1 (MQ=255)
           cTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccctg                                                                                         <  1:43093/100‑2 (MQ=255)
           cTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGAGCGCGGAACTCCGCGccctg                                                                                         <  2:1321346/100‑2 (MQ=255)
            tGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGccctgt                                                                                        <  1:3917470/100‑3 (MQ=255)
             ggggCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCtttt                                                                                       <  1:48032/100‑1 (MQ=255)
              gggCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCttttt                                                                                      >  2:627055/1‑100 (MQ=255)
               ggCTATTGCCGGCGTGGCAGCTGCTGTTTCTGGCCTGACGATAGGGATTATCGCTCTATCTAAATAACTGAACAGGGCGCGGAGTTCCGCGCTtttttt                                                                                      >  2:4001747/1‑99 (MQ=255)
               ggCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTa                                                                                     >  1:1900232/1‑100 (MQ=255)
               ggCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTa                                                                                     <  2:877867/100‑1 (MQ=255)
                 cTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAAc                                                                                   >  2:3622562/1‑100 (MQ=255)
                 cTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAAc                                                                                   <  1:1395704/100‑1 (MQ=255)
                 cTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCgcgccctgttcagt                                                                                   <  1:2917621/100‑11 (MQ=38)
                 cTAGTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAACTCCGCGCCCtgttcag                                                                                    <  2:970784/99‑4 (MQ=38)
                     tGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATa                                                                               >  1:2134298/1‑100 (MQ=255)
                      gCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATAc                                                                              >  1:3801205/1‑100 (MQ=255)
                      gCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATAc                                                                              <  2:543702/100‑1 (MQ=255)
                      gCCGGCGGGGCAGCTGCTGTTGCGGGCGTGACGATGGGGATCATCGCTCTATCGAAATAACGGGGCAGGGCGCGGAGTTCCGCGCTCGTTTGAAc                                                                                   <  2:3844718/95‑1 (MQ=255)
                        cggcggGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAg                                                                            <  2:1234211/100‑1 (MQ=255)
                        cggcggGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAg                                                                            >  1:3080791/1‑100 (MQ=255)
                         ngcggGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGt                                                                           >  1:3691554/2‑100 (MQ=255)
                         ggcggGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGt                                                                           <  2:270839/100‑1 (MQ=255)
                         ggcggGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGt                                                                           >  2:2971153/1‑100 (MQ=255)
                             gggCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtct                                                                       <  1:2033539/100‑1 (MQ=255)
                              ggCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtctc                                                                      <  1:3275076/100‑1 (MQ=255)
                              ggCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtctc                                                                      >  1:2600112/1‑100 (MQ=255)
                                cAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTAtctctc                                                                    >  1:183260/1‑100 (MQ=255)
                                 aGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGCAATAACTCAACAGCGCGCGCAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTcc                                                                   <  2:343538/100‑1 (MQ=255)
                                  gctgctGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCa                                                                  >  1:3313875/1‑100 (MQ=255)
                                  gctgctGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCa                                                                  >  2:1702724/1‑100 (MQ=255)
                                    tgctgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATa                                                                >  1:1026108/1‑100 (MQ=255)
                                    tgctgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATa                                                                >  2:1616514/1‑100 (MQ=255)
                                    tgctgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATa                                                                <  2:183260/100‑1 (MQ=255)
                                    tgctgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGATCCGCGCTCTTTTTAACGATACAGTATCTCTCCATa                                                                >  1:2944496/1‑100 (MQ=255)
                                     gctgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATaa                                                               <  1:1129438/100‑1 (MQ=255)
                                     gctgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATaa                                                               >  2:2212156/1‑100 (MQ=255)
                                     gctgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATaa                                                               <  2:1414567/100‑1 (MQ=255)
                                      ctgTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAg                                                              >  2:3635291/1‑100 (MQ=255)
                                        gTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGAt                                                            >  1:2644730/1‑100 (MQ=255)
                                          tGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaa                                                          <  1:2083092/100‑1 (MQ=255)
                                          tGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCAGCGCTCTTTTTAACGATAAAGTATGTTTACATAAGATaa                                                          >  2:1951320/1‑100 (MQ=255)
                                           gCTGGCCTGACGATGGGGATCATCGCTCTATCGATATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATaaa                                                         >  2:1667857/1‑100 (MQ=255)
                                             tGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAat                                                       >  1:2452864/1‑100 (MQ=255)
                                              ggCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAata                                                      >  2:3058525/1‑100 (MQ=255)
                                                ccTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAataga                                                    >  1:647539/1‑98 (MQ=255)
                                                ccTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATAtt                                                    >  2:564398/1‑100 (MQ=255)
                                                ccTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATAtt                                                    >  1:1373680/1‑100 (MQ=255)
                                                  tGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGc                                                  <  1:1554831/100‑1 (MQ=255)
                                                  tGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGc                                                  >  1:1497713/1‑100 (MQ=255)
                                                   gACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCt                                                 >  1:3628197/1‑100 (MQ=255)
                                                    atcttGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCta                                                <  2:4133485/96‑1 (MQ=255)
                                                    aCGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCta                                                <  1:1616518/100‑1 (MQ=255)
                                                     cGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCtat                                               >  2:2151306/1‑100 (MQ=255)
                                                       aTGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTatat                                             >  2:320905/1‑100 (MQ=255)
                                                        tGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATGAAATATTGCtaagat                                            >  1:4133514/1‑96 (MQ=255)
                                                          gggATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTct                                          <  2:431876/100‑1 (MQ=255)
                                                          gggATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTct                                          >  1:3487668/1‑100 (MQ=255)
                                                          gggATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTct                                          >  2:2851804/1‑100 (MQ=255)
                                                              tcatcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTctgctg                                      <  2:3886432/100‑1 (MQ=255)
                                                                atcGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTggg                                    <  1:912040/100‑1 (MQ=255)
                                                                atcGCTCTATCGAAATAACTGAACAAGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTggg                                    <  2:2502127/100‑1 (MQ=255)
                                                                    ctctATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAAt                                >  1:642131/1‑100 (MQ=255)
                                                                         tCGAAATAACTGAACAGTGCGCGGAGTTTCGCGATCTTTTTATCGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCct                           <  2:2346515/100‑1 (MQ=255)
                                                                         tCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCct                           >  2:1920649/1‑100 (MQ=255)
                                                                          cGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctg                          >  2:2762085/1‑100 (MQ=255)
                                                                           gAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgc                         >  1:2309312/1‑100 (MQ=255)
                                                                           gAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgc                         >  1:2765039/1‑100 (MQ=255)
                                                                            aaaTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgct                        <  1:1677492/100‑1 (MQ=255)
                                                                            aaaTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgct                        <  2:144564/100‑1 (MQ=255)
                                                                             aaTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgctg                       >  2:844006/1‑100 (MQ=255)
                                                                             aaTAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCctgctg                       >  2:4164348/1‑100 (MQ=255)
                                                                                aaCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGAtt                    >  1:1962232/1‑100 (MQ=255)
                                                                                 gttaaaaagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGTTATATTGTGCTGGGTAATTCCCTGCTGAttt                   >  2:2251437/11‑100 (MQ=38)
                                                                                 gttaaaaagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGAttt                   >  1:575805/11‑100 (MQ=38)
                                                                                 gttaaaaagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGAttt                   >  1:876175/11‑100 (MQ=38)
                                                                                 aCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGAttt                   <  1:1865334/100‑1 (MQ=255)
                                                                                 aCTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGAttt                   >  2:1171739/1‑100 (MQ=255)
                                                                                  cTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTa                  >  2:1855659/1‑100 (MQ=255)
                                                                                   tGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAg                 <  1:1778762/100‑1 (MQ=255)
                                                                                    gAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATGCAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGc                <  1:3502253/100‑1 (MQ=255)
                                                                                    gAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGAGACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGc                <  2:3919826/100‑1 (MQ=255)
                                                                                    aaaaagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATCTTGCTATATTGTACTGGTTAATTCCCTGCTGAtttctc                >  2:1071226/8‑97 (MQ=255)
                                                                                    aaaaagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGc                >  1:2699061/8‑100 (MQ=255)
                                                                                    aaaaagagcgcGGAACGCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTa                  >  2:3414817/8‑98 (MQ=255)
                                                                                       cAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATg             >  1:3652466/1‑100 (MQ=255)
                                                                                       cAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATg             >  2:1359562/1‑100 (MQ=255)
                                                                                       aagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATg             >  2:167976/5‑100 (MQ=255)
                                                                                       aagagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATg             >  1:510772/5‑100 (MQ=255)
                                                                                        agagcgcGGAACTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGt            >  1:1005165/4‑100 (MQ=255)
                                                                                                  gTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGt  <  2:2657206/100‑1 (MQ=255)
                                                                                                  gTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGt  >  2:2118752/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CTGGGAAACGGCTGGGGCTATTGCCGGCGGGGCAGCTGCTGTTGCTGGCCTGACGATGGGGATCATCGCTCTATCGAAATAACTGAACAGGGCGCGGAGTTCCGCGCTCTTTTTAACGATACAGTATCTCTCCATAAGATAAAATATTGCTATATTCTGCTGGGTAATTCCCTGCTGATTTAGCATGTTACATCGCGT  >  NC_000913/3418530‑3418727

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