Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,379,594	G→A	19.0%	A135V (GCC→GTC)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,594	0	G	A	19.0%	45.7 / 5.6	21	A135V (GCC→GTC)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand):  ref base G (9/8);  new base A (2/2);  total (11/10)
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.78e-01

CCGGCGGAATATTTGACGTCGCTACCAGGGTAATACCGCGAGCGAACAGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAG  >  NC_000913/3379463‑3379743                                                                                                                                    |                                                                                                                                                      ccGGCGGAATATTTGACGGCGCTACCAGGGTAATCCCGCGAGCGAACTGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCaa                                                                                                                                     <  1:245049/150‑1 (MQ=255)        aaTATTTGACGTCGCTACCAGGGTAATACCGCGAGCGAACAGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTc                                                                                                                              <  1:62089/150‑1 (MQ=255)            tttGACGTCGCTACCAGGGTAATACCGCGAGCGAACAGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGACTTTAAAGCGATCGGCAATAATTTCCAGc                                                                                                                          >  1:51466/1‑150 (MQ=255)            tttGACGTCGCTACCAGGGTAATACCGCGAGCGAACAGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGACTTTAAAGCGATCGGCAATAATTTCCAGc                                                                                                                          >  1:51465/1‑150 (MQ=255)                 cGTCGCTACCAGGGTAATACCGCGAGCGAACAGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATc                                                                                                                     <  1:81629/150‑1 (MQ=255)                                    ccGCGAGCGAACAGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCa                                                                                                  <  2:93963/150‑1 (MQ=255)                                      gcgagcgaACAGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGt                                                                                                <  2:43742/150‑1 (MQ=255)                                          gcgaACAGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGc                                                                                            >  1:73766/1‑150 (MQ=255)                                            gaACAGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGctc                                                                                          >  2:177048/1‑150 (MQ=255)                                                aGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCa                                                                                      >  1:61090/1‑150 (MQ=255)                                                  ggCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATg                                                                                    >  1:202766/1‑150 (MQ=255)                                                        catcaGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTACGGCAGTTAGCTCTTCATGCACCCg                                                                              >  1:21775/1‑150 (MQ=255)                                                        catcaGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCg                                                                              >  1:39917/1‑150 (MQ=255)                                                                 gCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCGCAGCATaaa                                                                     >  2:11387/1‑150 (MQ=255)                                                                                 tCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGACTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAgg                                                     <  2:51465/150‑1 (MQ=255)                                                                                 tCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGACTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAgg                                                     <  2:51466/150‑1 (MQ=255)                                                                                            aGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCg                                          <  2:61090/150‑1 (MQ=255)                                                                                             gAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGc                                         <  1:89904/150‑1 (MQ=255)                                                                                                  aaaaaaTTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTcc                                    >  1:11611/1‑150 (MQ=255)                                                                                                      aaTACGACAAAACAGAGCACGACAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCAGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACCGTGAAAGTGCAGGCGCTGTTTCCCCACTcccccc                                <  2:202766/150‑4 (MQ=255)                                                                                                                           tCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTcc           >  1:206999/1‑150 (MQ=255)                                                                                                                                    gCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCTCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAg  <  1:177048/150‑1 (MQ=255)                                                                                                                                    |                                                                                                                                                      CCGGCGGAATATTTGACGTCGCTACCAGGGTAATACCGCGAGCGAACAGGGCTTTCATCAGACCGCCAAGTAGCATGGCATCGGTAATATCAGAAACAAAAAATTCGTCAAAACAGAGCACGTCAGTTTCGGCTTTAAAGCGATCGGCAATAATTTCCAGCGGATCGGTCTGCCCCTGTAAGGCAGTTAGCTCTTCATGCACCCGCAGCATAAAACGGTGAAAGTGCAGGCGCTGTTTCCGCTCTCCCGGCAGGCTTTGATAGAAAAGGTCCATCAGCCAG  >  NC_000913/3379463‑3379743

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑