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)
               aGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAgcg                                                                                                                         >  2:375999/1‑139 (MQ=255)
                gTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAgcgc                                                                                                                        <  1:1261738/139‑1 (MQ=255)
                    tGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAgcgc                                                                                                                        <  2:389357/135‑1 (MQ=255)
                    tGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAgcgc                                                                                                                        >  1:389357/1‑135 (MQ=255)
                    tGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAAc                                                                                                                    >  1:224143/1‑139 (MQ=255)
                    tGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAAc                                                                                                                    >  1:2782/1‑139 (MQ=255)
                     gCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAAcc                                                                                                                   <  2:1190991/139‑1 (MQ=255)
                      cTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCa                                                                                                                  <  1:375999/139‑1 (MQ=255)
                        gACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAGCGCCAACCATAg                                                                                                               >  1:1134648/1‑139 (MQ=255)
                           tgagggCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCCACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAAcc                                                                                                                   >  2:749926/4‑133 (MQ=255)
                           tGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAAcc                                                                                                                   <  1:749926/133‑1 (MQ=255)
                           tGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGtt                                                                                                             >  2:835578/1‑139 (MQ=255)
                           tGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGtt                                                                                                             >  1:538072/1‑139 (MQ=255)
                                   tCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCAtt                                                                                                     >  2:565587/1‑139 (MQ=255)
                                    cGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCGGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                                                                   <  2:1134648/139‑1 (MQ=255)
                                     gACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                                                                   >  2:720404/1‑139 (MQ=255)
                                              cGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTGTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCGCGGGt                                                                                          >  1:804385/1‑139 (MQ=255)
                                                tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                                                        <  1:835578/139‑1 (MQ=255)
                                                tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                                                        <  2:829415/139‑1 (MQ=255)
                                                  gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         <  1:508098/89‑1 (MQ=255)
                                                  gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         >  2:508098/1‑89 (MQ=255)
                                                        gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                                                                >  1:4573/1‑139 (MQ=255)
                                                        gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                                                                >  1:610778/1‑139 (MQ=255)
                                                          gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTTATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         <  1:621295/81‑1 (MQ=255)
                                                          gaCAGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         >  2:621295/1‑81 (MQ=255)
                                                               ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                                                                  >  1:604948/1‑114 (MQ=255)
                                                               ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                                                                  <  2:604948/114‑1 (MQ=255)
                                                                cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCACCGTGAAGTGGCGa                                                                        <  1:565587/139‑1 (MQ=255)
                                                                           tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                                                             <  2:538072/139‑1 (MQ=255)
                                                                             tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTAATCGTGAAGTGGCGAAAGCCACAGATgc                                                           >  1:811658/1‑139 (MQ=255)
                                                                               gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                                                               <  2:567729/117‑1 (MQ=255)
                                                                               gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                                                               >  1:567729/1‑117 (MQ=255)
                                                                               gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                                                         >  1:726102/1‑139 (MQ=255)
                                                                               gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCAAAAGCCACAGATgcgc                                                         <  2:804385/139‑1 (MQ=255)
                                                                                    ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag                                                    <  1:1178537/139‑1 (MQ=255)
                                                                                      aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg                                                  <  2:4573/139‑1 (MQ=255)
                                                                                         gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt                                               <  2:610778/139‑1 (MQ=255)
                                                                                             cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGcc                                                                                                                       <  2:555445/63‑1 (MQ=255)
                                                                                             cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGcc                                                                                                                       >  1:555445/1‑63 (MQ=255)
                                                                                                  ggCCGATGGTTGGTTTCTAAGGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATGGTTTATCCATTTGCGCCCGGGTTTCAGCGGGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGa                                      <  2:764124/139‑1 (MQ=255)
                                                                                                  ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         >  2:1106468/1‑41 (MQ=255)
                                                                                                  ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                                                         <  1:1106468/41‑1 (MQ=255)
                                                                                                    ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACg                                    >  1:566219/1‑139 (MQ=255)
                                                                                                    ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTCTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACg                                    >  2:883432/1‑139 (MQ=255)
                                                                                                     cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGa                                   <  2:1018528/139‑1 (MQ=255)
                                                                                                       aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAAc                                                                                                                    <  1:1213700/56‑1 (MQ=255)
                                                                                                       aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAAc                                                                                                                    >  2:1213700/1‑56 (MQ=255)
                                                                                                       aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATa                                 >  2:1182210/1‑139 (MQ=255)
                                                                                                         ggTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAAt                               <  2:1178197/139‑1 (MQ=255)
                                                                                                              aTTTCTAACGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAaca                          <  2:230635/139‑1 (MQ=255)
                                                                                                                     aCGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAgcgc                                                                                                                        <  2:424358/38‑1 (MQ=255)
                                                                                                                     aCGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAgcgc                                                                                                                        >  1:424358/1‑38 (MQ=255)
                                                                                                                     aCGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAgc                                                                                                                          <  2:180287/36‑1 (MQ=255)
                                                                                                                     aCGTCGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAgc                                                                                                                          >  1:180287/1‑36 (MQ=255)
                                                                                                                         cGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGcc                                                                                                                       <  2:138372/35‑1 (MQ=37)
                                                                                                                         cGGTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGcc                                                                                                                       >  1:138372/1‑35 (MQ=37)
                                                                                                                          ggTGAATTTTCTCGTTTCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                                                                <  2:1239245/58‑1 (MQ=255)
                                                                                                                          ggTGAATTTTCTCGTTTCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                                                                >  1:1239245/1‑58 (MQ=255)
                                                                                                                          ggTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg                                          >  2:1034348/1‑111 (MQ=255)
                                                                                                                          ggTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg                                          <  1:1034348/111‑1 (MQ=255)
                                                                                                                          ggTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACg              <  2:726102/139‑1 (MQ=255)
                                                                                                                          ggTGAATTTTCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACg              <  2:224143/139‑1 (MQ=255)
                                                                                                                                 tttCTCGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTa       <  1:1182210/139‑1 (MQ=255)
                                                                                                                                   tctcGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACg              >  2:1161921/1‑130 (MQ=255)
                                                                                                                                   tctcGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACg              <  1:1161921/130‑1 (MQ=255)
                                                                                                                                     tcGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc                                         <  1:267696/101‑1 (MQ=255)
                                                                                                                                     tcGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc                                         >  2:267696/1‑101 (MQ=255)
                                                                                                                                      cGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATca  <  1:720404/139‑1 (MQ=255)
                                                                                                                                      cGTTTC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCATTACATca  <  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: ‑