Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,379,266 (C)7→8 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.C94.6% 368.9 / 13.1 112coding (739/1128 nt)zapEdivisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base C (47/59);  total (50/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353
                                                                                                      |                                                                                                  
cGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                     >  2:1094534/1‑100 (MQ=255)
  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cc                                                                                                 >  1:805228/1‑100 (MQ=255)
     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccccc                                                                                            >  2:1915206/1‑100 (MQ=255)
      gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccccc                                                                                           <  2:2112696/100‑1 (MQ=255)
      gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑cccccc                                                                                           >  1:431870/1‑100 (MQ=255)
       cGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccccccc                                                                                         <  1:1829815/100‑1 (MQ=255)
         tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccccccccg                                                                                       <  1:2407886/100‑3 (MQ=255)
           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccccccccgcc                                                                                     >  2:566052/1‑96 (MQ=255)
           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccccccccgcc                                                                                     <  1:1915179/100‑5 (MQ=255)
           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccccccccgcc                                                                                     >  2:2141189/1‑96 (MQ=255)
             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAg                                                                                    >  2:1063709/1‑100 (MQ=255)
             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAg                                                                                    >  1:1849268/1‑100 (MQ=255)
             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAgc                                                                                   >  2:1948879/1‑100 (MQ=255)
                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccccccccgcccgcg                                                                                 >  1:728170/1‑91 (MQ=255)
                  agaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCa                                                                               <  2:1827638/100‑1 (MQ=255)
                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTATTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAcc                                                                            <  2:2036524/100‑1 (MQ=255)
                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAcc                                                                            <  1:900620/98‑1 (MQ=255)
                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCAt                                                                          <  2:32642/100‑1 (MQ=255)
                        ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATa                                                                         <  1:1099113/100‑1 (MQ=255)
                         cAGCGGCTGGGTCTCGACCCCCCATGTTTCCAATTGTCGATGGGTGGATTCCAACCTCCGGTAAATTTTCCGGTTCCCCCCCCCGCCAGCGCCAAc                                                                             <  1:1623062/96‑1 (MQ=255)
                          gtcgagaacttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGt                                                                       >  1:1023365/10‑100 (MQ=255)
                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTATCGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGtt                                                                      <  2:2010726/100‑1 (MQ=255)
                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGtt                                                                      <  1:2192755/100‑1 (MQ=255)
                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGttt                                                                     <  1:1164325/100‑1 (MQ=255)
                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCCTAGTTTAt                                                                   >  1:746827/1‑100 (MQ=255)
                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATcc                                                                 >  1:1278264/1‑100 (MQ=255)
                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACCCCCGGGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATc                                                                  <  1:480883/99‑1 (MQ=255)
                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCa                                                                <  1:1799439/100‑1 (MQ=255)
                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCa                                                                <  2:2539050/100‑1 (MQ=255)
                                  gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATccct                                                               >  2:928547/1‑98 (MQ=255)
                                  gTTCTCGACGCCCATTGTTGCTAATGACCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCAt                                                               >  2:1766930/1‑100 (MQ=255)
                                     ctcGACGCCCATTGTTGCTAATGGCTGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1597044/100‑1 (MQ=255)
                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:2103472/100‑1 (MQ=255)
                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:2163424/100‑1 (MQ=255)
                                      tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                           >  1:920117/1‑100 (MQ=255)
                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTTATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:2206028/99‑1 (MQ=255)
                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  1:2190513/1‑100 (MQ=255)
                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGAGCCAACCATAGTTTATCAATTTgcg                                                          >  1:2458040/1‑100 (MQ=255)
                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCACCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  1:1534996/1‑100 (MQ=255)
                                       cGACGCCCATTGTTGCTAATGGCCGAAGGTTGCTTTCTAACGTCGGTGAATTTTCCCGTTTCCCCCCCCGCCAGCGCCAACCATAGCTTCTCGATTTgc                                                           >  1:2247485/1‑99 (MQ=255)
                                         aCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGcc                                                        <  2:2398484/100‑1 (MQ=255)
                                         aCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGcc                                                        >  2:2286185/1‑100 (MQ=255)
                                           gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGCTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      >  2:259626/1‑100 (MQ=255)
                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  1:2394969/1‑100 (MQ=255)
                                                 tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                >  1:829672/1‑100 (MQ=255)
                                                 tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                <  1:1241698/100‑1 (MQ=255)
                                                 tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                >  2:2052416/1‑100 (MQ=255)
                                                    tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  2:2421331/1‑100 (MQ=255)
                                                    tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  1:1313566/1‑100 (MQ=255)
                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  1:522039/100‑1 (MQ=255)
                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         <  2:155432/100‑1 (MQ=255)
                                                            gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     >  1:2204412/1‑100 (MQ=255)
                                                            gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     <  2:1702703/100‑1 (MQ=255)
                                                            gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTC‑‑CCTCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     <  1:2042825/100‑1 (MQ=255)
                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  2:1023376/100‑1 (MQ=255)
                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   >  1:1687127/1‑100 (MQ=255)
                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  2:1269233/100‑1 (MQ=255)
                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             >  1:1300602/1‑100 (MQ=255)
                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             >  2:2592564/1‑100 (MQ=255)
                                                                      gATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           <  1:2140826/100‑1 (MQ=255)
                                                                        tttCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  1:1064312/100‑1 (MQ=255)
                                                                        tttCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  2:888292/100‑1 (MQ=255)
                                                                         ttCTAACGTCGGTGAATTTTCTCGTGTC‑‑CCCCCCCGCCAGCGCCAACCATAGGTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        <  1:1208022/100‑1 (MQ=255)
                                                                         ttCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        <  2:746832/100‑1 (MQ=255)
                                                                         ttCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        >  1:1902310/1‑100 (MQ=255)
                                                                         ttCTAACGTCGGTGAATATTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        <  2:1936934/100‑1 (MQ=255)
                                                                           cTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      <  2:1858578/100‑1 (MQ=255)
                                                                            tAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     >  2:1910713/1‑100 (MQ=255)
                                                                            tAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    >  2:369794/1‑100 (MQ=255)
                                                                             aaCGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    <  2:1350694/100‑1 (MQ=255)
                                                                              aCGTCGGTGAATTTTCTCGTTCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     <  1:1026639/98‑1 (MQ=255)
                                                                                gTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  2:2190540/100‑1 (MQ=255)
                                                                                   cGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               >  1:2811539/1‑100 (MQ=255)
                                                                                   cGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  2:575902/100‑1 (MQ=255)
                                                                                     gTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             >  2:1650675/1‑100 (MQ=255)
                                                                                     gTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  1:2769711/100‑1 (MQ=255)
                                                                                     gTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             <  1:1043240/100‑1 (MQ=255)
                                                                                     gTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             >  2:43681/1‑100 (MQ=255)
                                                                                     gTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             >  2:2762625/1‑100 (MQ=255)
                                                                                      tGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            <  2:1849294/100‑1 (MQ=255)
                                                                                        aaTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga          <  2:2217894/100‑1 (MQ=255)
                                                                                         aTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:1393739/100‑1 (MQ=255)
                                                                                         aTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  2:862779/1‑100 (MQ=255)
                                                                                         aTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:885237/100‑1 (MQ=255)
                                                                                         aTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:2245888/100‑1 (MQ=255)
                                                                                         aTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  1:1895391/1‑100 (MQ=255)
                                                                                         aTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  1:2273759/1‑100 (MQ=255)
                                                                                         aTTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCCCCCGCGTTTCATCCTGAAGTGGCTAAATCCACAGATGCGCCtgattgag         >  1:2423467/1‑100 (MQ=255)
                                                                                          ttttCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt        <  1:1391459/100‑1 (MQ=255)
                                                                                           tttCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGTTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:728176/98‑1 (MQ=255)
                                                                                           tttCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:829678/100‑1 (MQ=255)
                                                                                           tttCTCGTTT‑‑CCCCCCTCGCCAGCGCCAACCATAGTTTATCCAGTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:2893635/98‑1 (MQ=255)
                                                                                            ttCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  2:1980143/1‑100 (MQ=255)
                                                                                             tctcGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     >  1:566551/1‑100 (MQ=255)
                                                                                             tctcGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     >  1:1398555/1‑100 (MQ=255)
                                                                                             tctcGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  2:1826344/100‑1 (MQ=255)
                                                                                               tcGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  2:1202047/100‑1 (MQ=255)
                                                                                                cGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGACGTGGCGAAAACCACAGATGCGCCTGAGTGAGTGTACg   >  1:542742/1‑99 (MQ=255)
                                                                                                cGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  2:2449292/100‑1 (MQ=255)
                                                                                                cGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  1:1870730/1‑100 (MQ=255)
                                                                                                      |                                                                                                  
CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑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: