Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,379,266:1 (C)7→8 100% coding (732/1128 nt) zapE ← divisome ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,379,2591.C92.8% 322.6 / 19.5 97coding (739/1128 nt)zapEdivisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (6/1);  new base C (44/46);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353
                                                                                                      |                                                                                                    
cGCACAGCGTAGTCAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGtt‑‑t                                                                                                       >  2:1614570/1‑100 (MQ=255)
cGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGtt‑‑t                                                                                                       <  1:2423366/100‑1 (MQ=255)
 gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c                                                                                                   <  1:37225/100‑1 (MQ=255)
  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑cc                                                                                                 <  2:1453652/100‑1 (MQ=255)
   acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccc                                                                                                >  1:205511/1‑100 (MQ=255)
    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑cccc                                                                                              >  2:317233/1‑100 (MQ=255)
    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑cccc                                                                                              >  1:1932518/1‑100 (MQ=255)
     aGCGTAGTAAAAGAGACCGCCAGCTTCTTGTTCTCGACTTCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccc                                                                                            <  1:1517765/100‑1 (MQ=255)
      gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑cccccc                                                                                           >  1:2334714/1‑100 (MQ=255)
           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccccccgcc                                                                                     <  1:675782/100‑5 (MQ=255)
           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccccccgcc                                                                                     <  1:1688696/100‑5 (MQ=255)
            taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CCCCACCCGCCa                                                                                     >  1:890282/1‑100 (MQ=255)
               tccgatctgCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAgcg                                                                                  <  2:1814796/92‑1 (MQ=255)
               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAgcg                                                                                  <  1:1128769/100‑1 (MQ=255)
               aaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAgcg                                                                                  <  2:729308/100‑1 (MQ=255)
                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAgcgc                                                                                 <  1:348406/100‑1 (MQ=255)
                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGcc                                                                                >  2:665319/1‑100 (MQ=255)
                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGcc                                                                                >  2:657109/1‑100 (MQ=255)
                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCCGCGCCa                                                                               >  2:2562652/1‑99 (MQ=255)
                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCaa                                                                              >  2:1886201/1‑100 (MQ=255)
                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAAcc                                                                            <  1:1306962/100‑1 (MQ=255)
                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCa                                                                           <  1:1358294/100‑1 (MQ=255)
                      cGCCAGCGTCTGGTTCTCGACCCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGGGCCAACCa                                                                           >  2:359030/1‑100 (MQ=255)
                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAAcc                                                                            >  2:2410073/1‑99 (MQ=255)
                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAgcgagatcgga                                                                          >  1:1814805/1‑92 (MQ=255)
                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCAt                                                                          >  2:52046/1‑100 (MQ=255)
                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CCCCACCGCGCCAGCGCCAACCa                                                                           >  2:415582/1‑100 (MQ=255)
                         taGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAg                                                                        <  1:1898022/99‑1 (MQ=255)
                          aGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAga                                                                       >  2:1898013/1‑99 (MQ=255)
                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGtt                                                                      <  2:888949/100‑1 (MQ=255)
                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGttt                                                                     >  2:2304867/1‑100 (MQ=255)
                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGttt                                                                     >  2:2485562/1‑100 (MQ=255)
                             gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTa                                                                    <  1:2393495/100‑1 (MQ=255)
                             gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CCCACCCCGCCAGCGCCAACCATAGTTTa                                                                    >  1:875614/1‑100 (MQ=255)
                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTAt                                                                   <  2:293471/100‑1 (MQ=255)
                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATCTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTa                                                                    >  2:202505/1‑99 (MQ=255)
                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATcc                                                                 >  1:2592057/1‑100 (MQ=255)
                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGCTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATcc                                                                 <  2:1214669/100‑1 (MQ=255)
                                   ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCAtt                                                              >  1:918744/1‑100 (MQ=255)
                                   ttCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCAtt                                                              <  1:2479849/100‑1 (MQ=255)
                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  2:454509/100‑1 (MQ=255)
                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                            <  1:1942185/100‑1 (MQ=255)
                                       cGACGCCGATTGTTGCTAATGGCCGATGGTTGATTTCTCACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CGCTCCACGCCAGCGCGAACCATATTTTATCCATTTg                                                            >  2:2519812/1‑98 (MQ=255)
                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  1:1208276/1‑100 (MQ=255)
                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  1:255754/1‑100 (MQ=255)
                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACAATAGTTTATCCATTTgcg                                                          >  2:1726797/1‑100 (MQ=255)
                                          cGCCCATTGTTGCTAATGGCCGATGGTTCAGTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACACTAGCTCATCCCTTTGcccc                                                        >  2:353333/1‑99 (MQ=255)
                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      >  1:1535622/1‑100 (MQ=255)
                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     <  1:1629583/100‑1 (MQ=255)
                                              cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                                   <  1:292564/100‑1 (MQ=255)
                                                    tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  1:176417/1‑100 (MQ=255)
                                                      cTAATTGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCgg                                           <  2:1161199/100‑2 (MQ=255)
                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           <  1:1846653/100‑1 (MQ=255)
                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       <  2:1882986/100‑1 (MQ=255)
                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       >  2:1734504/1‑100 (MQ=255)
                                                           ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAgn                                      <  1:283876/100‑2 (MQ=255)
                                                            gCCAATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     <  2:593934/100‑1 (MQ=255)
                                                             ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    >  2:2556987/1‑100 (MQ=255)
                                                             ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    >  1:2428043/1‑100 (MQ=255)
                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  1:1060578/100‑1 (MQ=255)
                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCGTAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:733383/100‑1 (MQ=255)
                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:1208576/100‑1 (MQ=255)
                                                               gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:1719179/100‑1 (MQ=255)
                                                                  ggTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                               >  1:2558497/1‑100 (MQ=255)
                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             >  2:822838/1‑100 (MQ=255)
                                                                     tGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            <  2:31790/100‑1 (MQ=255)
                                                                     tGATTTCTAACGTCGGGGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCAGGGTTTCATCGTGAAGTGGCGAAAGcc                            <  1:1546615/100‑1 (MQ=255)
                                                                      gATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  2:464994/1‑100 (MQ=255)
                                                                      gATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CCCTCCATTCCAGCGCTAAACATAGTTTATTCATTTGCGCCCAGGTTTCATCGTGAAGTGGCGAAAGCca                           >  2:608983/1‑100 (MQ=255)
                                                                        tttCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          <  1:2410086/100‑1 (MQ=255)
                                                                             aaCGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    >  2:2327892/1‑100 (MQ=255)
                                                                              aCGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     >  2:488444/1‑100 (MQ=255)
                                                                               cGTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                  >  1:2294644/1‑100 (MQ=255)
                                                                                gTCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 >  2:23077/1‑100 (MQ=255)
                                                                                 tCGGTGAATTTTCTCGTTCCC‑‑‑CCCCCCCGCCCGCCCCCACCATTGTTTTTCCCTTTTCGCCCCGGCCTCCTCCCTCAGTGGCGaaaac                             >  2:1057847/1‑86 (MQ=255)
                                                                                 tCGGTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCt                <  2:1493587/100‑1 (MQ=255)
                                                                                   ggTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga              <  1:23077/100‑1 (MQ=255)
                                                                                   ggTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               >  2:53059/1‑99 (MQ=255)
                                                                                   ggTGAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCATGAAGTGGCGAAAGCCACAGATGCGCCtga              <  1:2588276/100‑1 (MQ=255)
                                                                                     tgtctttcCTCGTT‑‑CC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            <  1:125682/96‑1 (MQ=255)
                                                                                      gAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  1:188666/100‑1 (MQ=255)
                                                                                      gAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  1:88321/100‑1 (MQ=255)
                                                                                      gAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  2:2002570/100‑1 (MQ=255)
                                                                                      gAATTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  1:692866/100‑1 (MQ=255)
                                                                                       aaTTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga          >  2:592321/1‑100 (MQ=255)
                                                                                        aTTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:2437351/100‑1 (MQ=255)
                                                                                        aTTTTCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  2:820100/1‑100 (MQ=255)
                                                                                         ttttCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:1535613/100‑1 (MQ=255)
                                                                                         ttttCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGCTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  2:2115787/1‑99 (MQ=255)
                                                                                          tttCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  1:982109/100‑1 (MQ=255)
                                                                                          tttCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATACATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       >  1:247866/1‑100 (MQ=255)
                                                                                           ttCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  1:797081/1‑100 (MQ=255)
                                                                                           ttCTCGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  2:963835/1‑100 (MQ=255)
                                                                                            tctcGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  1:2415128/100‑1 (MQ=255)
                                                                                             ctcGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    >  1:293723/1‑100 (MQ=255)
                                                                                             ctcGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    >  2:1871254/1‑100 (MQ=255)
                                                                                             ctcGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    >  1:961753/1‑100 (MQ=255)
                                                                                              tcGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATACTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  1:908441/100‑1 (MQ=255)
                                                                                               cGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  1:384012/1‑100 (MQ=255)
                                                                                               cGTT‑‑TC‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATGTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  2:535449/100‑1 (MQ=255)
                                                                                                      |                                                                                                    
CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑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: