Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,379,266 (C)7→9 coding (732/1128 nt) zapE ← cell division factor ZapE with ATPase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,379,2591.C89.7% 85.1 / 4.0 30coding (739/1128 nt)zapEcell division factor ZapE with ATPase activity
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (11/15);  total (13/17)
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
*NC_0009133,379,2592.C90.0% 90.8 / 4.0 30coding (739/1128 nt)zapEcell division factor ZapE with ATPase activity
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (12/15);  total (13/17)
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

GTGAGAGCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCA  >  NC_000913/3379121‑3379384
                                                                                                                                             ||                                                                                                                                
gTGAGAGCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c                                                                                                                                <  2:315961/140‑1 (MQ=255)
    gagCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccc                                                                                                                            <  2:714744/140‑1 (MQ=255)
       cgcAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccc‑‑ccc                                                                                                                       >  1:682936/1‑139 (MQ=255)
       cgcAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccc‑‑ccc                                                                                                                       >  1:241697/1‑139 (MQ=255)
         cAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccc‑‑ccccg                                                                                                                     <  1:10128/140‑4 (MQ=255)
          aaTATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccc‑‑ccccgc                                                                                                                    >  1:31395/1‑136 (MQ=255)
          aaTATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccc‑‑ccccgc                                                                                                                    >  1:490469/1‑136 (MQ=255)
                gTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAgcgc                                                                                                                <  1:596813/140‑1 (MQ=255)
                     gCTGACTGCGGGCGTCGACCCCCAGCGTAGTAAAAGAGACCGCCAGCGTTTGGTTCTCGACGCCCTTTTTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAAcc                                                                                                           <  2:196294/140‑1 (MQ=255)
                      cTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGATCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCa                                                                                                          >  1:687937/1‑140 (MQ=255)
                         aCTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccc‑‑ccccgc                                                                                                                    <  1:226137/125‑5 (MQ=255)
                         aCTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccc‑‑ccccgc                                                                                                                    <  1:74224/125‑5 (MQ=255)
                         aCTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccc‑‑ccccgc                                                                                                                    >  2:226137/1‑121 (MQ=255)
                         aCTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccc‑‑ccccgc                                                                                                                    >  2:74224/1‑121 (MQ=255)
                          cTGCGGGCGTCGCCGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGt                                                                                                      <  2:712478/140‑1 (MQ=255)
                          cTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGt                                                                                                      <  2:235295/140‑1 (MQ=255)
                                   tCGACCCACGGTGTAGTAAAAGAGACCGCCGGCGTCTGGTTCTCGACGCCAATTGTTGCTAATGGCCGAGGGTTGATTTCTAACGCCGGTGATTTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCAtt                                                                                             <  2:31395/140‑1 (MQ=255)
                                       cGCACAGCGTAGTAAAAAAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TTCCCCCC‑‑CCCCCCGCGCCCACCATAGTTTATCCCTTTgcg                                                                                         >  2:205876/1‑140 (MQ=255)
                                                           aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                                                     >  2:425341/1‑140 (MQ=255)
                                                           aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CCCCCCACCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                                                     >  1:474429/1‑140 (MQ=255)
                                                             cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                                                 <  2:382058/126‑1 (MQ=255)
                                                             cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                                                 >  1:382058/1‑126 (MQ=255)
                                                             cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                                                   >  2:622554/1‑140 (MQ=255)
                                                             cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                                                   >  2:131949/1‑140 (MQ=255)
                                                                 ggCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                                                               <  1:556202/139‑1 (MQ=255)
                                                                     tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                                                           >  2:155686/1‑140 (MQ=255)
                                                                            ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CCCCCCACCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                                                    <  2:474429/140‑1 (MQ=255)
                                                                                aCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                                                <  2:687937/140‑1 (MQ=255)
                                                                                    ctaTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag                                            <  2:312178/138‑1 (MQ=255)
                                                                                    ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag                                            >  1:429279/1‑140 (MQ=255)
                                                                                     cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt                                           <  1:205876/140‑1 (MQ=255)
                                                                                       ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCTGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga                                         >  2:10777/1‑140 (MQ=255)
                                                                                                aTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc                                 >  1:571376/1‑140 (MQ=255)
                                                                                                         ggTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAAt                       >  1:132249/1‑140 (MQ=255)
                                                                                                                ttCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACAcc                <  1:622554/140‑1 (MQ=255)
                                                                                                                     aCGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGt           >  2:151101/1‑140 (MQ=255)
                                                                                                                      cGTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTc          <  1:425341/140‑1 (MQ=255)
                                                                                                                       gTCGGTGAATTTTCTCGTTCCC‑‑‑CCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTcc         <  2:682936/140‑1 (MQ=255)
                                                                                                                       gTCGGTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTcc         <  2:241697/140‑1 (MQ=255)
                                                                                                                          ggTGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACg      <  2:132249/140‑1 (MQ=255)
                                                                                                                            tGAATTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGtt    <  1:131949/140‑1 (MQ=255)
                                                                                                                              aaTTTTCTCGTT‑‑TCCCCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCa  <  1:151101/140‑1 (MQ=255)
                                                                                                                                             ||                                                                                                                                
GTGAGAGCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CCCCC‑‑CCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTTCA  >  NC_000913/3379121‑3379384

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