Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,379,266:1 (C)7→8 100% coding (732/1128 nt) zapE ← divisome ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,379,2591.C92.2% 423.4 / 32.8 141coding (739/1128 nt)zapEdivisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (4/7);  new base C (65/64);  total (69/72)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 161 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353
                                                                                                       |                                                                                                   
 gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑c                                                                                                  >  1:2134326/1‑100 (MQ=255)
 gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGACCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑Gttt                                                                                                      >  2:572176/1‑99 (MQ=255)
  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑cc                                                                                                 >  2:136861/1‑100 (MQ=255)
  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑cc                                                                                                 >  2:3752667/1‑100 (MQ=255)
  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGACGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑cc                                                                                                 >  2:3405704/1‑100 (MQ=255)
  cacaGAGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑cc                                                                                                 <  2:3039165/100‑1 (MQ=255)
     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑ccccc                                                                                            <  2:1622804/100‑1 (MQ=255)
       cGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑ccccccc                                                                                         <  1:3642259/100‑1 (MQ=255)
        gtagtaAAAGAGACCGCCAGCGTCTGGTTCTTGACGTTCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑cccccccc                                                                                        <  1:3529433/100‑2 (MQ=255)
        gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGATTTTTCTC‑‑GTTT‑‑‑cccccccc                                                                                        <  2:1067815/100‑2 (MQ=255)
           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑ccccccccgcc                                                                                     >  2:2727908/1‑96 (MQ=255)
           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑ccccccccgcc                                                                                     >  2:1814867/1‑96 (MQ=255)
             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGATAAATTTTCTC‑‑ATTT‑‑‑CCCACCCCACCAg                                                                                    >  1:2573873/1‑100 (MQ=255)
             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGACGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAg                                                                                    >  1:3317355/1‑100 (MQ=255)
                agagACCGCCAGCGTCTGGTTCTCTACGCCCATTGTTGGTAATGTCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAgcgc                                                                                 <  2:2594393/100‑1 (MQ=255)
                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAgcgc                                                                                 <  2:2259325/100‑1 (MQ=255)
                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGcc                                                                                <  2:2824981/100‑1 (MQ=255)
                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGcc                                                                                >  2:3872875/1‑100 (MQ=255)
                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCaa                                                                              >  1:2209853/1‑100 (MQ=255)
                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGCTGAATTTTCTC‑‑GTTTCCCCCCCCCGCCAGCGCCa                                                                               >  2:2261154/1‑100 (MQ=255)
                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTTAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAAcc                                                                            >  2:2336460/1‑100 (MQ=255)
                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGACAGCGCCAAc                                                                             >  2:767482/1‑99 (MQ=255)
                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCAt                                                                          <  1:3581374/100‑1 (MQ=255)
                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCAt                                                                          >  2:2770058/1‑100 (MQ=255)
                           gCGTCTGGTTCTCGACGCCCATTGTTGTTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGtt                                                                      <  1:2799301/100‑1 (MQ=255)
                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTAt                                                                   <  1:1863538/100‑1 (MQ=255)
                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATc                                                                  >  1:3484012/1‑100 (MQ=255)
                                tGGTTCTCGACGCCCATTGTTGTTAATGGCTGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATcc                                                                 <  1:654290/100‑1 (MQ=255)
                                tGGTTCTCGACGCCCATTGTTGCTAATGGTCGATGGTTGATTTATAACGTCGGTGAATGTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATcc                                                                 <  1:425879/100‑1 (MQ=255)
                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATcc                                                                 <  2:1357082/100‑1 (MQ=255)
                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATcc                                                                 <  2:1084156/100‑1 (MQ=255)
                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCa                                                                >  1:990871/1‑100 (MQ=255)
                                   ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCAtt                                                              >  2:3757436/1‑100 (MQ=255)
                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑CCACCCCCGCCAGCGCCAACAATAGTTTATCCAttt                                                             >  1:927446/1‑99 (MQ=255)
                                      tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                           >  2:1006509/1‑100 (MQ=255)
                                      tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                           >  1:1258442/1‑100 (MQ=255)
                                      tcGACGCCCATTGTTGCTAATGGCCGATGGTTAATTCCTAACGTCGGTAAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:422496/99‑1 (MQ=255)
                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  1:1531484/1‑100 (MQ=255)
                                          cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGcc                                                        <  2:498031/99‑1 (MQ=255)
                                           gCCCATTGTTGCTAATGTCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      <  2:463213/100‑1 (MQ=255)
                                             ccgttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTCTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACAATAGTTTATCCATTTACGACCggg                                                    >  2:103285/4‑100 (MQ=255)
                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  1:1428887/100‑1 (MQ=255)
                                               aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                  <  2:3811945/100‑1 (MQ=255)
                                                ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                  <  1:2838111/99‑1 (MQ=255)
                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCa                                               >  1:933341/1‑97 (MQ=255)
                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  1:2826540/100‑1 (MQ=255)
                                                      cTAATGGCCGCTGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACTATTGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           <  2:2580634/100‑1 (MQ=255)
                                                      cTAATGGCCGATGGTTGCTTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCCACCATAGTTTATCAATTTGCGCCCGCGTTTCATCGt                                           >  2:3168594/1‑100 (MQ=255)
                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           >  1:2381664/1‑100 (MQ=255)
                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTACCCATTTGCGCCCGGGTTTCATc                                             >  1:3846036/1‑98 (MQ=255)
                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  1:1692970/1‑100 (MQ=255)
                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          <  1:2581746/100‑1 (MQ=255)
                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         <  2:3600446/100‑1 (MQ=255)
                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCCCCCGCGCCAACCATAGTTTACCCATTTGCGCCCGGGTTTCATCGTGCAg                                       >  2:899679/1‑100 (MQ=255)
                                                           ggCCTAAGGGTGATTTCTTAAGTCGGTTAATTTTTTC‑‑CTTTCCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                        <  2:1476259/100‑1 (MQ=255)
                                                            gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     >  1:2959944/1‑100 (MQ=255)
                                                            gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     >  1:447464/1‑100 (MQ=255)
                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   >  1:2029934/1‑100 (MQ=255)
                                                                aTGGTTGATTTCTAACGTCGGTTAATTTTTTTTTTTTT‑‑‑TTCTTTCTCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:1054134/100‑1 (MQ=255)
                                                                aTGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGa                                 <  2:239699/100‑1 (MQ=255)
                                                                 tGGTTGATTTCTAACGTCGGTGAAATTTCTC‑‑GTTTC‑‑CCCCCCCTCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                                >  2:2650899/1‑100 (MQ=255)
                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTC‑‑TTTTC‑‑CCCCCCCGCCAGCGCCAAACATAGTTTATCCATGTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               <  1:2472368/100‑1 (MQ=255)
                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGAGCCCGGGTTTCATCGTGAAGTGGCGaaa                               >  1:394159/1‑100 (MQ=255)
                                                                     tGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCGATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            >  2:2640470/1‑100 (MQ=255)
                                                                      gATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑CCACCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  1:3723740/1‑100 (MQ=255)
                                                                       aTTTCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          <  2:3483973/100‑1 (MQ=255)
                                                                        tttCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  1:1006520/100‑1 (MQ=255)
                                                                        tttCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGAGCCAACCATAGCTTATCCATTTGCGCCCGGGTTTCCTCCTGAAGTGACGAAAGCcaca                         >  2:2517385/1‑100 (MQ=255)
                                                                        tttCTAACGTCGGTGAATTTGCTC‑‑GTTA‑‑‑CCTCCCCGCCCGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          <  2:1976300/98‑1 (MQ=255)
                                                                         ttttAACGTCTGTAATTTTTTTC‑‑TTTT‑‑‑CCCTCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        <  1:1079478/97‑1 (MQ=255)
                                                                         ttCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        >  1:3798220/1‑100 (MQ=255)
                                                                         ttCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCATAAAGTGGCGAAAGCcaca                         >  1:921902/1‑99 (MQ=255)
                                                                          tCTAACGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       <  2:2454709/100‑1 (MQ=255)
                                                                              aCGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  2:415888/100‑1 (MQ=255)
                                                                               cGTCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCACCCGAGTTTAACCGTAACGTGGCGAAACACACCGATacgc                  >  1:3672210/1‑100 (MQ=255)
                                                                                 tCGGTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:3752706/99‑1 (MQ=255)
                                                                                  cGGTGAATTTTCTC‑‑TTTTC‑‑CCCCCCCGCCAGCGCCAATCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  1:1536832/100‑1 (MQ=255)
                                                                                   ggTGAATTTTCTC‑‑GTTTCCCCCCCCCGCCAGCTCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  2:2746593/100‑1 (MQ=255)
                                                                                    gTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCATCCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  1:1581481/100‑1 (MQ=255)
                                                                                    gTGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  1:3191987/100‑1 (MQ=255)
                                                                                     tGAATTTTCTT‑‑GTTTC‑‑CCCCCCCGCCATCGCTAACCATTGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  1:384249/99‑1 (MQ=255)
                                                                                     tGAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            >  2:2488767/1‑100 (MQ=255)
                                                                                      gAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgaggg           >  2:461213/1‑98 (MQ=255)
                                                                                      gAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAATAGCGAAAGCCACAGATGCGCCtgagtg           >  1:3363066/1‑100 (MQ=255)
                                                                                      gAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTCAAGTGGCGAAAGCCACACATGCGCCtgagtg           <  1:2640494/100‑1 (MQ=255)
                                                                                      gAATTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCACCAGAGTTTCAACATGAAATAGCGAAAGCCACAGATGCGCCtgagt            >  1:187077/1‑99 (MQ=255)
                                                                                        aTTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAATCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:1164260/100‑1 (MQ=255)
                                                                                        aTTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  1:715508/1‑100 (MQ=255)
                                                                                        aTTTTCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  2:1044978/1‑100 (MQ=255)
                                                                                         ttttCTC‑‑GTTTT‑‑‑CCCCCCCGCCTGCGTCATCCATAGTTTATCCATTTGCGGCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:2737060/99‑1 (MQ=255)
                                                                                         ttttCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  2:3803882/1‑99 (MQ=255)
                                                                                         ttttCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        <  1:956634/100‑1 (MQ=255)
                                                                                         ttttCTC‑‑GTTT‑‑‑CCCCCTCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:2736616/99‑1 (MQ=255)
                                                                                           ttCTC‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  2:3102077/1‑100 (MQ=255)
                                                                                           ttCTC‑‑GTTTC‑‑CCCCCCCGCAAGCGCCAACCATAGTTTATCAATTTGCGCACTGGTTTCATAGTGAAGTGGCAAAACCCACAGAAGCGCCtcagtcaaggt      >  2:1840128/1‑94 (MQ=255)
                                                                                            tatc‑‑GTTTG‑‑CCCCCCCGCCAGCGCCAAGCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAGGTGGCGAAAGCTACAGATGCGCCTGagtgagtg       <  2:2345849/96‑1 (MQ=255)
                                                                                             ctc‑‑GTTTC‑‑CCCCCCCTCCAGCGCTAACCATAGTTTATCCATTTGCGCCCGAGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  1:289090/99‑1 (MQ=255)
                                                                                             ctc‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    <  2:1231799/100‑1 (MQ=255)
                                                                                               c‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  2:1668996/100‑1 (MQ=255)
                                                                                               c‑‑GTTTC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  2:2802192/1‑100 (MQ=255)
                                                                                                       |                                                                                                   
CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTC‑‑GTTT‑‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353

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