Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,641,929	T→A	I306F (ATC→TTC)	rodZ ←	cytoskeleton protein RodZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,641,929	0	T	A	95.8%	72.8 / ‑6.4	25	I306F (ATC→TTC)	rodZ	cytoskeleton protein RodZ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (14/9);  minor base C (0/1);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGATCTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCG  >  NC_000913/2641707‑2642153                                                                                                                                                                                                                                |                                                                                                                                                                                                                                 aTGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTc                                                                                                                                                                                                      >  1:20134/1‑251 (MQ=255)    gACTGTACGGCGATGGGAGCTCCACCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTACGCATGAAAAAGCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCGTGATACT‑GGAACAGTACGGCGGCTGGCGCACCAATTTTCAGt                                                                                                                                                                                                   <  2:52249/251‑1 (MQ=255)        gTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGttt                                                                                                                                                                                                 <  2:20134/249‑1 (MQ=255)         tACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAACTACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACACGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGt                                                                                                                                                                                               >  1:205682/1‑250 (MQ=255)          aCGGCGATGGGAGCACCTTCGCCAATTGGCACAGTCCCAACTTAAATACGTTTTGATGTGCTACTGTGAATTGGAGCCTGGTTATGCATGGAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCTGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCGGGATACT‑GGAACTGTACTGCGGCTGGCGCACCTTTTTTCAGTTTGTAc                                                                                                                                                                                             <  2:282439/251‑1 (MQ=255)          aCGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGAGCGACAGGTTTCCCTTGAGACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGATTGTAc                                                                                                                                                                                             <  2:98108/251‑1 (MQ=255)                   ggAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTgg                                                                                                                                                                                    >  1:172290/1‑251 (MQ=255)                        aCCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGt                                                                                                                                                                               <  2:51591/251‑1 (MQ=255)                                                          cGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCATGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCg                                                                                                                                             >  2:183440/1‑251 (MQ=255)                                                                        aCGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGGTACCCGTCTGTCAATGCGCCGGTGATTGTTCGGCATTGAGAGTCAGACGGGCAAACTGGTTAGTTCTGATAAAACGACTGAGATAGACAGGTTTCCCTTGATACT‑GGAACTGAACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTACCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAAtttttt                                                                                                                               >  1:271650/1‑251 (MQ=255)                                                                             gAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCgg                                                                                                                          >  2:248394/1‑251 (MQ=255)                                                                                                      aaaTCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACTACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCTCCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAGCAAGGTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCa                                                                                                 <  1:206942/251‑1 (MQ=255)                                                                                                                 gTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAAGGACAGACGCGCAACCTGGTTAGATCTGATAAAACGACTCAGATCGACAGGATTCCCTTGATAAT‑GGAACTGTACTGCGGCTGGCGTACCCATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTAACGTCTTTACGCGGCATACCGCTAAACAATTTTTTACCGGTAGCATAAGTGACATCCAGCCAGCAATCGGCAGTa                                                                                       >  1:185621/1‑250 (MQ=255)                                                                                                                    aCCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAg                                                                                    <  2:74841/250‑1 (MQ=255)                                                                                                                                       ggTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCgg                                                                >  1:202143/1‑251 (MQ=255)                                                                                                                                                   gagctGAGGGGCAGACGCGCAACCTTGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCGCTGGATACT‑GGACCTGTACTGCGTCTTGCGCACCAATTTTCAGTTTGTACGTTGCTTTGCCGGTTAAGTTCAAATTTCCGTCTGTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCAGCAGTTACCTCCAGCCACCAATCGGCAGTAAATTTCTGCACCAGCGCACGCGGATCAGCCACCgg                                                    <  2:281580/247‑1 (MQ=255)                                                                                                                                                        gAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAACAGACTCAGATCGACAGGCTTAACTTGATAAT‑GGAACTGTACTGAGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCATGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAAAAATTTTTTAACGGTAGCATCAGTGAACTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCgg                                                    >  1:281580/1‑246 (MQ=255)                                                                                                                                                        gAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATAGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTAAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTgg                                               >  1:220501/1‑251 (MQ=255)                                                                                                                                                                    cgccaCCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCCCACCCATTTTCAGTTTGTACGGAGCCTGGCCGGTTAAGTTAAAATTACCGTCTTTACGCTGCATACCGATAAACAATTTTTTACAGGTAGCATCAGTGACCTCCAGCAAGCCATCGGCAGTAAAGTTCATCAACAGCGCATTCGGAACAGCCACCGGCGCGGTCACGCCAGCCt                                   >  2:264062/5‑251 (MQ=255)                                                                                                                                                                            ggTTAGTTCTGATAACACGACTCAGAGCGACAGGTTGCCCGTGATACTTGGA‑CTGTACTTCTGCTGCAGCTCCTATTTTCTGGTTGTACGGGTCCTGGCCGGTTAAGGTCAAATTACCGTTTTGACGCTGCATACCGCTAAACAAATTTTTACCGGTAGCATCAGTGACCGCCAGCCAGCAATCGGCATTAAAGGTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCCGCCTGATCGGtt                           <  2:11484/251‑1 (MQ=255)                                                                                                                                                                            ggTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGtt                           <  2:215750/251‑1 (MQ=255)                                                                                                                                                                                       aTAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGgcgc                >  1:111310/1‑251 (MQ=255)                                                                                                                                                                                       aTAAAACGACTCAGATCGACAGGTTTCCATTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGgcgc                >  2:157069/1‑251 (MQ=255)                                                                                                                                                                                         aaaaCGACTCAGATCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGAGCGGTTGGCAACGGCgcag              <  2:57583/251‑1 (MQ=255)                                                                                                                                                                                                     aTCGACAGGTTTCCCTTGATACT‑GGAACTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCg  >  2:262444/1‑251 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                 ATGGACTGTACGGCGATGGGAGCACCATCGCCAATCGGCACATTCCCAACGTAAATACGTGTTGATTTTCTACGTTGAATTGGAGCCTGGTTATGCATGAAAAATCTCCCGCGTTACCCGTCTGTTACTGCGCCGGTGATTGTTCGGCATTGAGGGTCAGACGCGCAACCTGGTTAGTTCTGATAAAACGACTCAGATCGACAGGTTTCCCTTGATACT‑GGATCTGTACTGCGGCTGGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCG  >  NC_000913/2641707‑2642153

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑