Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,379,768:1 +C 100% coding (230/1128 nt) zapE ← divisome ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,379,7641.C55.8% 312.7 / 300.5 163coding (234/1128 nt)zapEdivisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (40/32);  new base C (50/41);  total (90/73)
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 177 total aligned reads displayed.

CCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCT  >  NC_000913/3379665‑3379863
                                                                                                    |                                                                                                     
 cGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCg                                                                                                        <  2:681382/99‑1 (MQ=255)
    agcaTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGcccc                                                                                                  <  1:3623778/100‑1 (MQ=255)
     gcaTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGccccc                                                                                                 >  1:1939632/1‑99 (MQ=255)
      caTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGccccca                                                                                                <  1:1683052/100‑3 (MQ=255)
        tAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGcccccaca                                                                                              <  2:3674963/100‑5 (MQ=255)
        tAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGcccccaca                                                                                              >  1:284741/1‑96 (MQ=255)
         aaaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGcccccacat                                                                                             >  1:714882/1‑95 (MQ=255)
         aaaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACata                                                                                            >  1:474270/1‑100 (MQ=255)
          aaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGcccccaca                                                                                              >  1:1537175/1‑94 (MQ=255)
           aaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACatat                                                                                           <  1:1171906/100‑1 (MQ=255)
             cGGTGAAAGTGCAGGCGCTGTTTCCGCACTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCACCCCACAtatata                                                                                        >  2:4114165/1‑100 (MQ=255)
              ggggAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACAtatata                                                                                        <  2:1928234/97‑1 (MQ=255)
              ggTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACAtatata                                                                                        <  2:3604271/100‑1 (MQ=255)
                tGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAg                                                                                      <  2:1152931/100‑1 (MQ=255)
                 gAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGACAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACCCCGCCCCCACATATATAAGc                                                                                     >  1:432681/1‑100 (MQ=255)
                 gAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGc                                                                                     >  1:1428626/1‑100 (MQ=255)
                 gAAAGTGCAGGCGCCGTGTACGCTAACGCGGCAGTCTTGGATCAACAACGTCCACCAGCCAGGTTTTCCCCCCTACTACACCGCCACcacacacatcca                                                                                      >  2:3391727/1‑91 (MQ=255)
                     gTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACg                                                                                 >  1:2655586/1‑100 (MQ=255)
                     gTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGc                                                                                >  1:1818712/1‑100 (MQ=255)
                     gTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGc                                                                                >  1:3990036/1‑100 (MQ=255)
                      tGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGc                                                                                >  2:1904552/1‑100 (MQ=255)
                      tGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGc                                                                                >  1:2661741/1‑100 (MQ=255)
                       gCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCGATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCa                                                                               >  2:2285781/1‑100 (MQ=255)
                       gCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCAc                                                                              >  1:2037064/1‑100 (MQ=255)
                       gCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCAc                                                                              <  1:2806714/100‑1 (MQ=255)
                          ggCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTg                                                                            >  1:3343098/1‑100 (MQ=255)
                            cgcTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCg                                                                         <  2:2595307/100‑1 (MQ=255)
                              cTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGt                                                                        >  1:2037106/1‑100 (MQ=255)
                               tGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTa                                                                       <  1:2936287/100‑1 (MQ=255)
                               tGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTa                                                                       <  1:909594/100‑1 (MQ=255)
                               tGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTAt                                                                      <  2:581500/100‑1 (MQ=255)
                                gTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATg                                                                     >  1:2745206/1‑100 (MQ=255)
                                   tCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCt                                                                   <  1:587720/100‑1 (MQ=255)
                                     cGCTCTCCCGGCAGGCTTTGATAGTAAAGGTCCATCAGCCGGGTTTCCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtg                                                                <  2:138887/100‑1 (MQ=255)
                                     cGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCttt                                                                 >  2:1790773/1‑100 (MQ=255)
                                      gCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTtg                                                                <  1:2271717/100‑1 (MQ=255)
                                      gCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgt                                                               >  1:2612477/1‑100 (MQ=255)
                                      gCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgt                                                               >  1:2459499/1‑100 (MQ=255)
                                       ctctCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgt                                                               <  1:2688870/100‑1 (MQ=255)
                                        tctcCGGGCAGGCTTTGATAGAAAAGGTCCCTCAGCCAGGTTTTCCCGCGTACTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgtg                                                              <  2:614549/100‑1 (MQ=255)
                                        tctcCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgtg                                                              <  1:641329/100‑1 (MQ=255)
                                        tctcCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTtgtg                                                              <  1:2946459/100‑1 (MQ=255)
                                          tcCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTc                                                            <  1:2277066/100‑1 (MQ=255)
                                          tcCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTc                                                            >  1:1724178/1‑100 (MQ=255)
                                           cccGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCAGATATATAAGCCACGCACTGGCGTATGCTTTGTGTCt                                                           >  1:3941345/1‑100 (MQ=255)
                                             cGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTc                                                         >  1:1021313/1‑100 (MQ=255)
                                             cGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTc                                                         >  2:3732455/1‑100 (MQ=255)
                                                cAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTcgcg                                                      <  1:1306309/100‑1 (MQ=255)
                                                cAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGt                                                     <  1:1230430/100‑1 (MQ=255)
                                                  ggCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGtc                                                    <  2:2847837/100‑2 (MQ=255)
                                                   gCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGttt                                                   >  2:2760852/1‑100 (MQ=255)
                                                    cTTTGATAGAAAAGGTCCATCAGCCAGTTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAc                                                 >  2:653852/1‑100 (MQ=255)
                                                    cTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTa                                                  >  2:1047397/1‑100 (MQ=255)
                                                       tGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcccc                                              <  1:1215707/100‑1 (MQ=255)
                                                        gctAGAAAGGGTCCATCAGCCAGGTTTTCCCGCGTTATACAGCGCACCCAATATATAAGCCACGCACTGGAGTATGCACTGTGTCTTCGCGTTAACCGca                                             >  2:31524/3‑100 (MQ=255)
                                                        gATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcccc                                              >  2:212766/1‑100 (MQ=255)
                                                        gATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAcccc                                              >  2:1842639/1‑100 (MQ=255)
                                                          tAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCcaca                                           <  2:1320499/100‑1 (MQ=255)
                                                          tAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCcaca                                           >  1:1077714/1‑100 (MQ=255)
                                                           aGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCcaca                                           <  1:3832394/100‑1 (MQ=255)
                                                           aGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATCTAAGCCACGCCCTGGCGTATGCTTTGTGTCTTCCCGTTTACCCCACAg                                          >  1:812871/1‑100 (MQ=255)
                                                            gAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAg                                          >  1:220534/1‑100 (MQ=255)
                                                             aaaaGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGc                                         <  2:1342207/100‑1 (MQ=255)
                                                               aaGGTCCATCAGCCAGGTTTTCCCGCTTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGc                                         <  2:3614850/98‑1 (MQ=255)
                                                               aaGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTTTTCGCGTTTACCCCACAGCTTa                                      <  2:2400688/100‑1 (MQ=255)
                                                                 ggTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTCGTGTCTTCGCGTTTACCCCACAGCTTAc                                     <  2:3844832/100‑1 (MQ=255)
                                                                 ggTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTAcc                                    <  2:1614301/100‑1 (MQ=255)
                                                                 ggTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTAcc                                    >  2:3845043/1‑100 (MQ=255)
                                                                  gTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTAcc                                    <  1:543167/100‑1 (MQ=255)
                                                                  gTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCg                                   >  1:3083515/1‑100 (MQ=255)
                                                                    ccATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGa                                  >  1:1610338/1‑100 (MQ=255)
                                                                      aTCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAcc                                <  1:1108406/100‑1 (MQ=255)
                                                                       tCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAccc                               >  2:2697623/1‑100 (MQ=255)
                                                                         aGCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGC                             >  2:2825978‑M1/1‑100 (MQ=255)
                                                                         aGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCG                            <  1:1336949‑M1/100‑1 (MQ=255)
                                                                          gCCAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCG                            >  1:3599807‑M1/1‑100 (MQ=255)
                                                                           ccAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGC                           <  2:3941319‑M1/100‑1 (MQ=255)
                                                                            cAGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCC                          <  2:2061082‑M1/100‑1 (MQ=255)
                                                                            cAGGTTGTCCCGCGTCCGACAGCGCCCCCACCTATATACGCCACGCACTGGCGTATGCGTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCgcgc                           <  2:3518858/99‑1 (MQ=255)
                                                                             aGGTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCA                         <  1:3668168‑M1/100‑1 (MQ=255)
                                                                              ggTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCAt                        >  1:2091760/1‑100 (MQ=255)
                                                                               gtcttCCCGCGTCCTATACCGCTCCTCCTATCTAATCCTCTCACCGTCGTATGCTATGTCCCTTCGCGTTTTCCCCACAGCTTACCTACCCGCGCCAt                        <  2:3201178/95‑1 (MQ=255)
                                                                               ggtttCCCGCGTCCTACACCGCCCCGACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCTCAGCTTACCGACCCGCGCCAtt                       <  2:224485/98‑1 (MQ=255)
                                                                               gTTTTCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCAtt                       >  2:1665931/1‑100 (MQ=255)
                                                                                 tttCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAg                     <  1:953241/100‑1 (MQ=255)
                                                                                  ttCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGcc                   <  2:3286313/100‑1 (MQ=255)
                                                                                  ttCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGcc                   <  1:1496081/100‑1 (MQ=255)
                                                                                   tCCCGCGTCCTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGcc                   >  1:3998868/1‑100 (MQ=255)
                                                                                     ccGCGTCCTACACCGCCCCACATATATAATCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTATCCCACAGCTTACCGACCCGCGCCATTAGCCCAc                <  2:2481934/100‑1 (MQ=255)
                                                                                           ccTACACCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTc           >  2:1528519/1‑100 (MQ=255)
                                                                                              acacCGCCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTg        >  2:555240/1‑100 (MQ=255)
                                                                                              acacCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTgg       >  1:2699348/1‑100 (MQ=255)
                                                                                              acacCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTgg       >  1:3437907/1‑100 (MQ=255)
                                                                                               cacCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTggg      <  1:3151464/100‑1 (MQ=255)
                                                                                                acCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTgggg     <  1:3087579/100‑1 (MQ=255)
                                                                                                 ccGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTggggg    >  2:27900/1‑100 (MQ=255)
                                                                                                 ccGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTggggg    >  2:2234417/1‑100 (MQ=255)
                                                                                                 ccGCCCCACATATATAAGCCACGCACTGGCGTAAGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTggggg    >  2:2485875/1‑100 (MQ=255)
                                                                                                  cGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGc   <  2:3032319/100‑1 (MQ=255)
                                                                                                   gCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCt  <  1:1381777/100‑1 (MQ=255)
                                                                                                    |                                                                                                     
CCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCCACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCT  >  NC_000913/3379665‑3379863

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: