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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,764	1	.	C	96.7%	99.4 / ‑3.1	30	coding (234/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (15/14);  total (16/14)
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

CCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCG‑CCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCT  >  NC_000913/3379665‑3379856                                                                                                     |                                                                                              ccGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCg                                                                                                <  2:1134400/100‑1 (MQ=255)    cagcaTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCG‑ccc                                                                                            <  2:1029624/100‑1 (MQ=255)       caTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGTCAGGTTTTCCCGCGTCCTACACCG‑ccc‑cca                                                                                        <  2:146649/100‑3 (MQ=255)       caTAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCG‑ccc‑cca                                                                                        >  2:1007396/1‑98 (MQ=255)         tAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCG‑ccc‑ccaca                                                                                      >  2:955093/1‑96 (MQ=255)          aaaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCTGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCG‑ccc‑ccacat                                                                                     >  2:1122484/1‑95 (MQ=255)           aaaCGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACata                                                                                     >  1:234292/1‑100 (MQ=255)              cGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACatatat                                                                                  >  1:1184124/1‑100 (MQ=255)                 taaaaGTGCAGGCACTGTTTCCGCTCTCCCGACAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAg                                                                               >  1:878218/3‑100 (MQ=255)                 tGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAg                                                                               <  2:921035/100‑1 (MQ=255)                      gTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACg                                                                          >  2:750051/1‑100 (MQ=255)                               cTGTTTCCGCTCTTCCGGCAGGCTTTGATTGAAAAGGTCCATCAGCCAGGTTTTCCTGTGTTCTACATTT‑CCCTCACATATATAAGCCACGCACTGGc                                                                   >  2:868840/1‑98 (MQ=255)                                tggttCCGCACTCCCGGCAGGCCTTGTTAGAAAAGCGCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCGCTGGCGTa                                                                <  2:936567/97‑1 (MQ=255)                                 gTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTAt                                                               <  1:325988/100‑1 (MQ=255)                                 gTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTAt                                                               <  2:1060725/100‑1 (MQ=255)                                    tCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCt                                                            >  2:870944/1‑100 (MQ=255)                                      cGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCttt                                                          <  2:307645/100‑1 (MQ=255)                                       gCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTtg                                                         <  2:667986/100‑1 (MQ=255)                                          ctcCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTtgtg                                                       <  1:402765/99‑1 (MQ=255)                                             ccGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCtt                                                   >  1:538671/1‑100 (MQ=255)                                             ccGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCtt                                                   >  1:538249/1‑100 (MQ=255)                                               ggCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTcg                                                 >  1:525433/1‑100 (MQ=255)                                                   ggCTTTGATAGAAAAGGTCCATCAGCCAGGTTATCCCGCCTCATACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGtt                                             <  2:914996/100‑1 (MQ=255)                                                    gCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGttt                                            >  2:279459/1‑100 (MQ=255)                                                        tGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTAccc                                        >  2:442113/1‑100 (MQ=255)                                                                     ccATCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGa                           >  1:554833/1‑100 (MQ=255)                                                                       aTCAGCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGAc                          >  2:423278/1‑99 (MQ=255)                                                                           gCCAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCcgcg                     >  1:1058749/1‑100 (MQ=255)                                                                            ccAGGTTTTCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCcgcn                     <  2:480538/99‑2 (MQ=255)                                                                                  tttCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAg              <  1:870950/100‑1 (MQ=255)                                                                                    tCCCGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGcc            <  2:328857/100‑1 (MQ=255)                                                                                     cccGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGccc           >  1:156462/1‑100 (MQ=255)                                                                                     cccGCGTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGccc           <  1:868846/100‑1 (MQ=255)                                                                                         cgTCCTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTc       <  1:955101/100‑1 (MQ=255)                                                                                            ccTACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTc    <  2:435325/100‑1 (MQ=255)                                                                                              tACACCGCCCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCt  >  2:76438/1‑100 (MQ=255)                                                                                                     |                                                                                              CCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAGGTTTTCCCGCGTCCTACACCG‑CCC‑CACATATATAAGCCACGCACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCT  >  NC_000913/3379665‑3379856

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: ‑