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.C96.8% 218.9 / ‑0.4 63coding (739/1128 nt)zapEdivisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (28/33);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GTGAGAGCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCA  >  NC_000913/3379121‑3379392
                                                                                                                                           |                                                                                                                                      
gTGAGAGCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         <  2:786673/139‑1 (MQ=255)
 tGAGAGCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c                                                                                                                                      >  1:829415/1‑139 (MQ=255)
   agagCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑ccc                                                                                                                                    <  2:398868/139‑1 (MQ=255)
               aGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAgcg                                                                                                                         >  2:375999/1‑139 (MQ=255)
                gTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAgcgc                                                                                                                        <  1:1261738/139‑1 (MQ=255)
                    tGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAgcgc                                                                                                                        <  2:389357/135‑1 (MQ=255)
                    tGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAgcgc                                                                                                                        >  1:389357/1‑135 (MQ=255)
                    tGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAc                                                                                                                    >  1:224143/1‑139 (MQ=255)
                    tGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAc                                                                                                                    >  1:2782/1‑139 (MQ=255)
                     gCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAcc                                                                                                                   <  2:1190991/139‑1 (MQ=255)
                      cTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCa                                                                                                                  <  1:375999/139‑1 (MQ=255)
                        gACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAGCGCCAACCATAg                                                                                                               >  1:1134648/1‑139 (MQ=255)
                           tgagggCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCCACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAcc                                                                                                                   >  2:749926/4‑133 (MQ=255)
                           tGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAcc                                                                                                                   <  1:749926/133‑1 (MQ=255)
                           tGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGtt                                                                                                             >  2:835578/1‑139 (MQ=255)
                           tGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGtt                                                                                                             >  1:538072/1‑139 (MQ=255)
                                   tCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCAtt                                                                                                     >  2:565587/1‑139 (MQ=255)
                                    cGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCGGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                                                                   <  2:1134648/139‑1 (MQ=255)
                                     gACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                                                                   >  2:720404/1‑139 (MQ=255)
                                              cGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTGTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCGCGGGt                                                                                          >  1:804385/1‑139 (MQ=255)
                                                tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                                                        <  1:835578/139‑1 (MQ=255)
                                                tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                                                        <  2:829415/139‑1 (MQ=255)
                                                  gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         <  1:508098/89‑1 (MQ=255)
                                                  gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         >  2:508098/1‑89 (MQ=255)
                                                        gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                                                                >  1:4573/1‑139 (MQ=255)
                                                        gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                                                                >  1:610778/1‑139 (MQ=255)
                                                          gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTTATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         <  1:621295/81‑1 (MQ=255)
                                                          gaCAGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         >  2:621295/1‑81 (MQ=255)
                                                               ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                                                                  >  1:604948/1‑114 (MQ=255)
                                                               ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                                                                  <  2:604948/114‑1 (MQ=255)
                                                                cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCACCGTGAAGTGGCGa                                                                        <  1:565587/139‑1 (MQ=255)
                                                                           tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                                                             <  2:538072/139‑1 (MQ=255)
                                                                             tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTAATCGTGAAGTGGCGAAAGCCACAGATgc                                                           >  1:811658/1‑139 (MQ=255)
                                                                               gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                                                               <  2:567729/117‑1 (MQ=255)
                                                                               gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                                                               >  1:567729/1‑117 (MQ=255)
                                                                               gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                                                         >  1:726102/1‑139 (MQ=255)
                                                                               gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCAAAAGCCACAGATgcgc                                                         <  2:804385/139‑1 (MQ=255)
                                                                                    ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag                                                    <  1:1178537/139‑1 (MQ=255)
                                                                                      aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg                                                  <  2:4573/139‑1 (MQ=255)
                                                                                         gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt                                               <  2:610778/139‑1 (MQ=255)
                                                                                             cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGcc                                                                                                                       <  2:555445/63‑1 (MQ=255)
                                                                                             cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGcc                                                                                                                       >  1:555445/1‑63 (MQ=255)
                                                                                                  ggCCGATGGTTGGTTTCTAAGGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATGGTTTATCCATTTGCGCCCGGGTTTCAGCGGGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGa                                      <  2:764124/139‑1 (MQ=255)
                                                                                                  ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         >  2:1106468/1‑41 (MQ=255)
                                                                                                  ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         <  1:1106468/41‑1 (MQ=255)
                                                                                                    ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACg                                    >  1:566219/1‑139 (MQ=255)
                                                                                                    ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTCTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACg                                    >  2:883432/1‑139 (MQ=255)
                                                                                                     cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGa                                   <  2:1018528/139‑1 (MQ=255)
                                                                                                       aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAc                                                                                                                    <  1:1213700/56‑1 (MQ=255)
                                                                                                       aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAAc                                                                                                                    >  2:1213700/1‑56 (MQ=255)
                                                                                                       aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATa                                 >  2:1182210/1‑139 (MQ=255)
                                                                                                         ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAAt                               <  2:1178197/139‑1 (MQ=255)
                                                                                                              aTTTCTAACGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAaca                          <  2:230635/139‑1 (MQ=255)
                                                                                                                     aCGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAgcgc                                                                                                                        <  2:424358/38‑1 (MQ=255)
                                                                                                                     aCGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAgcgc                                                                                                                        >  1:424358/1‑38 (MQ=255)
                                                                                                                     aCGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAgc                                                                                                                          <  2:180287/36‑1 (MQ=255)
                                                                                                                     aCGTCGGTGAATTTTCTCGTTTCCCCCCCCGCCAgc                                                                                                                          >  1:180287/1‑36 (MQ=255)
                                                                                                                         cGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGcc                                                                                                                       <  2:138372/35‑1 (MQ=37)
                                                                                                                         cGGTGAATTTTCTCGTTTCCCCCCCCGCCAGCGcc                                                                                                                       >  1:138372/1‑35 (MQ=37)
                                                                                                                          ggTGAATTTTCTCGTTTCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                                                                <  2:1239245/58‑1 (MQ=255)
                                                                                                                          ggTGAATTTTCTCGTTTCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                                                                >  1:1239245/1‑58 (MQ=255)
                                                                                                                          ggTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg                                          >  2:1034348/1‑111 (MQ=255)
                                                                                                                          ggTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg                                          <  1:1034348/111‑1 (MQ=255)
                                                                                                                          ggTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACg              <  2:726102/139‑1 (MQ=255)
                                                                                                                          ggTGAATTTTCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACg              <  2:224143/139‑1 (MQ=255)
                                                                                                                                 tttCTCGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTa       <  1:1182210/139‑1 (MQ=255)
                                                                                                                                   tctcGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACg              >  2:1161921/1‑130 (MQ=255)
                                                                                                                                   tctcGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACg              <  1:1161921/130‑1 (MQ=255)
                                                                                                                                     tcGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc                                         <  1:267696/101‑1 (MQ=255)
                                                                                                                                     tcGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc                                         >  2:267696/1‑101 (MQ=255)
                                                                                                                                      cGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATca  <  1:720404/139‑1 (MQ=255)
                                                                                                                                      cGTTTCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATca  <  1:883432/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                      
GTGAGAGCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATCA  >  NC_000913/3379121‑3379392

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: