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.C95.8% 331.8 / 5.6 96coding (739/1128 nt)zapEdivisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base C (41/51);  total (42/54)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379161‑3379353
                                                                                                   |                                                                                                 
gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c                                                                                                 <  2:1926059/100‑1 (MQ=255)
 cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc                                                                                               >  2:2201313/1‑100 (MQ=255)
        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATTGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccccccctg                                                                                       <  1:1613713/100‑3 (MQ=255)
        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccccccccg                                                                                       <  1:2094356/100‑3 (MQ=255)
            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAg                                                                                    >  2:611631/1‑100 (MQ=255)
            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccccccccgcc                                                                                     >  1:1491075/1‑94 (MQ=255)
              aaGAGACCGCCAGCGTCTGGTTCTCGACGTCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAgcg                                                                                  <  2:2193143/100‑1 (MQ=255)
              aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAgcg                                                                                  <  1:1742857/100‑1 (MQ=255)
               agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAgcgc                                                                                 <  1:998768/100‑1 (MQ=255)
                  gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCCGCGCCaa                                                                              >  1:199817/1‑100 (MQ=255)
                  gaCCGCCAGCGTATGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCaa                                                                              >  1:1565792/1‑100 (MQ=255)
                   aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAc                                                                             <  2:892704/100‑1 (MQ=255)
                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTCACGTCGGTGAATTTTCTCGTTT‑‑cccccccccccagcggc                                                                               >  2:1024832/1‑89 (MQ=255)
                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAcc                                                                            >  1:253541/1‑100 (MQ=255)
                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAcc                                                                            >  1:709087/1‑100 (MQ=255)
                      gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCAt                                                                          >  2:1153899/1‑100 (MQ=255)
                           cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGttt                                                                     >  2:1935407/1‑100 (MQ=255)
                            gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCTCCAGCGCCAACCATAGtt                                                                      >  2:855097/1‑98 (MQ=255)
                            gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTa                                                                    >  1:2174641/1‑100 (MQ=255)
                              cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATc                                                                  <  1:2165293/100‑1 (MQ=255)
                                ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCCGCCAGCGCCAACCATAGTTTATcc                                                                 <  1:2305589/100‑1 (MQ=255)
                                 gTTCTCAACGCCCATTGTTTCTAATGGCCGATGGTTGATTTCCACCGTCGCTGAATTTTCTCGTGC‑‑CACCCCCCGCAAGCGCCCACCATAGCTTATc                                                                  <  2:1153211/97‑1 (MQ=255)
                                    ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1843268/100‑1 (MQ=255)
                                    ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1362762/100‑1 (MQ=255)
                                      cGACGCTCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          <  1:511300/100‑1 (MQ=255)
                                      cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  1:1773955/1‑100 (MQ=255)
                                       gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                         >  1:1627755/1‑100 (MQ=255)
                                         cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                         >  2:2064421/1‑99 (MQ=255)
                                          gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      >  2:334870/1‑100 (MQ=255)
                                           cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  2:74621/1‑100 (MQ=255)
                                            ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACAATAGTTTATCCATTTGCGCCCgg                                                     >  2:411496/1‑99 (MQ=255)
                                                   tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  2:1010524/1‑100 (MQ=255)
                                                     cTAATCGCCCATTGTTGATCTCTCACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             <  2:2099994/98‑1 (MQ=255)
                                                       aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         >  1:1691488/1‑100 (MQ=255)
                                                       aaTGGCCGATGGTTGATTTCAAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         >  2:1455416/1‑100 (MQ=255)
                                                         tttccGTTGGTTGATTTCTAACTTCGTTGAATCTTCTCGTTTCCCCCCCCTCCAGCTCCAACCTTAGTTTATCCTTTTTCGCCCGGGTTTCATCGTg                                          <  2:609170/94‑1 (MQ=255)
                                                         tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       <  2:329566/100‑1 (MQ=255)
                                                          ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      >  1:475460/1‑100 (MQ=255)
                                                          ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      >  1:475707/1‑100 (MQ=255)
                                                           gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                      <  2:1546802/99‑1 (MQ=255)
                                                            ccGATGGTTGATTTCTAACGTCGGTGTATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCACCCGGGTTTCATCGTGAAgggg                                    <  2:135937/100‑1 (MQ=255)
                                                            ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    >  1:1596696/1‑100 (MQ=255)
                                                            ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    <  1:854616/100‑1 (MQ=255)
                                                            ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    <  2:1596684/100‑1 (MQ=255)
                                                             cgatagctcaTTTCTACCGTCCGTAAATTTTCTCTTTTCCCCCCCCGCCAGCGCCAGCCATAGTTTATCCGTTTGCGCCCGGGTTTCATCGCGAAGTAGc                                   <  2:1662983/91‑1 (MQ=255)
                                                             cGATGGTTGATTTCCAACGTCGGTGAAATTTCTCGTTTCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    <  1:247929/100‑1 (MQ=255)
                                                              gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  2:2356391/100‑1 (MQ=255)
                                                              gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  >  2:728118/1‑100 (MQ=255)
                                                              gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  2:509088/100‑1 (MQ=255)
                                                              gATGGTTGATTTCTAACGTCGGGGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:269914/100‑1 (MQ=255)
                                                               aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:2313673/99‑1 (MQ=255)
                                                                tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCgca                                >  1:2314339/1‑98 (MQ=255)
                                                                tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                                >  1:529417/1‑100 (MQ=255)
                                                                 ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               <  1:2185891/100‑1 (MQ=255)
                                                                 ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTAAAGTGGCGaaa                               >  1:43526/1‑100 (MQ=255)
                                                                  gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  1:9102/1‑100 (MQ=255)
                                                                  gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  2:1535424/100‑1 (MQ=255)
                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             <  1:74622/100‑1 (MQ=255)
                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             <  1:576830/100‑1 (MQ=255)
                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCCTGAAGTGGCCaaa                               >  2:2193224/1‑98 (MQ=255)
                                                                    tGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            >  1:2197474/1‑100 (MQ=255)
                                                                    tGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTCCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            <  1:2267486/100‑1 (MQ=255)
                                                                     gATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  1:406521/1‑100 (MQ=255)
                                                                      aTTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGTCGAAAGCcac                          <  1:113748/100‑1 (MQ=255)
                                                                       tttCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  1:543828/100‑1 (MQ=255)
                                                                        ttCTAAGGTAGGTGAATTTTCTCGTTG‑‑CCCCCGCCGCCAGGGCCAACCATAGTTTATCCATGTGCGCGCGGGTTTCAGCGTGAAGTGGCGAAAGCCACAg                        <  2:826218/100‑1 (MQ=255)
                                                                        ttCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTCTCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        >  2:393199/1‑100 (MQ=255)
                                                                        gatctaCGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        <  1:2137772/95‑1 (MQ=255)
                                                                          cTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      <  1:1948669/100‑1 (MQ=255)
                                                                          cTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      >  2:1413911/1‑100 (MQ=255)
                                                                           tAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCACTTGCGCCCGGCTTTCATCGTGAAGTGGCTAAAGCCACAGacg                     <  1:1611157/100‑3 (MQ=255)
                                                                            aaCGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    <  1:2096836/100‑1 (MQ=255)
                                                                            aaCGTCGGTGAATTTTCTCGTTCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    <  1:1637303/100‑1 (MQ=255)
                                                                             aCGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGagatc                   >  2:2137755/1‑96 (MQ=255)
                                                                             aCGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    >  2:497854/1‑99 (MQ=255)
                                                                             aCGTCGGTCAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  2:181575/100‑1 (MQ=255)
                                                                              cGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGGGGCGAAAGCCACAGATgcg                   >  2:10020/1‑99 (MQ=255)
                                                                               gTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:713497/100‑1 (MQ=255)
                                                                               gTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:244944/100‑1 (MQ=255)
                                                                                 cGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  1:1153905/100‑1 (MQ=255)
                                                                                 cGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  2:1913068/100‑1 (MQ=255)
                                                                                  ggTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga              >  1:469431/1‑100 (MQ=255)
                                                                                   gTGACTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGTGTTTCTTCGTGAAGTGGCGAAAGCCACAGATGCGCCtgaa             <  2:2073056/100‑2 (MQ=255)
                                                                                   gTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  1:928862/100‑1 (MQ=255)
                                                                                    tGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            <  2:2150802/100‑1 (MQ=255)
                                                                                       aTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:866828/100‑1 (MQ=255)
                                                                                       aTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:1688630/100‑1 (MQ=255)
                                                                                       aTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:897046/100‑1 (MQ=255)
                                                                                        ttttCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        >  2:47037/1‑100 (MQ=255)
                                                                                         tttCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  1:995919/100‑1 (MQ=255)
                                                                                         tttCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:475458/100‑1 (MQ=255)
                                                                                         tttCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:475705/100‑1 (MQ=255)
                                                                                          ttCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  1:1041707/1‑100 (MQ=255)
                                                                                          ttCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  1:133973/1‑100 (MQ=255)
                                                                                           tctcGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  1:938821/100‑1 (MQ=255)
                                                                                             tcGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  1:533482/100‑1 (MQ=255)
                                                                                              cGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  1:2316262/100‑1 (MQ=255)
                                                                                              cGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  1:625226/1‑100 (MQ=255)
                                                                                              cGTTTCCCCCCCCGCCAGCGACAACCATAGTTTATCCATTTGCGCCCGGGTTTCCTCGTGAAGTGGCGAAAGCCAAAGATGCGCCTGAGTGAGTGTAc    >  2:633528/1‑98 (MQ=255)
                                                                                              cGTTT‑‑CCCCCTCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGTGATAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  2:1565781/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379161‑3379353

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: