Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,619,269	A→G	V26V (GTA→GTG)	rhsB →	Rhs protein with DUF4329 family putative toxin domain, putative neighboring cell growth inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,619,269	0	A	G	100.0%	76.7 / NA	24	V26V (GTA→GTG)	rhsB	Rhs protein with DUF4329 family putative toxin domain, putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/10);  total (14/10)

CTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCC  >  NC_000913/3619137‑3619326                                                                                                                                     |                                                          cttcGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAATGATTTTACGATGAGCGGCAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGTAtt                                                     <  1:783146/139‑1 (MQ=255)       tttttCGTATGCAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGcc                                               <  2:791681/139‑1 (MQ=255)               aTGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCgg                                       <  2:390073/139‑1 (MQ=255)                       aCTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTg                                                    <  2:549309/118‑1 (MQ=255)                       aCTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTg                                                    >  1:549309/1‑118 (MQ=255)                            cATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCg                          >  2:642358/1‑139 (MQ=255)                            cATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCg                          >  1:821353/1‑139 (MQ=255)                              ttAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGgt                        >  2:280163/1‑139 (MQ=255)                               tAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGgtg                       >  2:1157705/1‑139 (MQ=255)                                   ataataGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGgtg                       >  1:839303/1‑135 (MQ=255)                                   ataataGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGgtg                       <  2:839303/135‑1 (MQ=255)                                       taGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGcc                                               <  1:737337/107‑1 (MQ=255)                                       taGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCa                                           <  1:470398/111‑1 (MQ=255)                                       taGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCa                                           >  2:470398/1‑111 (MQ=255)                                       taGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcgg               >  2:660071/1‑139 (MQ=255)                                       taGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGAGcc                                               >  2:737337/1‑107 (MQ=17)                                            aaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGcc                   >  1:301352/1‑130 (MQ=255)                                            aaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGcc                   <  2:301352/130‑1 (MQ=255)                                             aGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGt         <  2:821353/139‑1 (MQ=255)                                                 ttttACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACg     >  2:290743/1‑139 (MQ=25)                                                 ttttACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACg     >  2:1204649/1‑139 (MQ=255)                                                    tACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcggc              <  2:648624/127‑1 (MQ=16)                                                    tACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcggc              >  1:648624/1‑127 (MQ=16)                                                    tACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTcc  >  1:1189913/1‑139 (MQ=16)                                                                                                                                     |                                                          CTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCC  >  NC_000913/3619137‑3619326

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