Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,769,447	G→A	E35K (GAA→AAA)	yibV →	PF15596 family protein YibV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,769,447	0	G	A	100.0%	29.3 / NA	13	E35K (GAA→AAA)	yibV	PF15596 family protein YibV
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/5);  total (8/5)

GACAATTTCAGAGAGAGTTAGATTTGGCAGACAAAGCAATATTATGGGCACTAATATCAGCCAGTACTAAAGAGGGACGAAAAGCGTGTTCACTAAGCTATTTTGCCTGTAAGGCCGCAGAGGCTGAACTCGGACTGGCATATATGGCTGCGAATGATAACAAGGAATTTCTCACATCCTTATCTAATATAATGAGGTACAAAATAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAAC  >  NC_000913/3769322‑3769640                                                                                                                              |                                                                                                                                                                                                  gACAATTTCAGAGAGAGTTAGATTTGGCAGACACAGCAATATTATGGGCAATAATATCCGCCAGTACTCAAGAGGGACGACAAGCGTGTTCACTAAGCTATTTTGCCTGTCAGGCTGCAGAGGCTAAACTCGGACTGGCATATATGACTGCGAATGATAACAAGGAATTTCTCACATCCTTATCTAATATAATGAGGTACAAAATAg                                                                                                                  >  1:317761/1‑207 (MQ=255)                          ggCAGACAAAGCAATATTATGGGCACTAATATCAGCCAGTACTAAAGAGGGACGAAAAGCGTGTTCACTAAGCTATTTTGCCTGTAAGGCCGCAGAGGCTAAACTCg                                                                                                                                                                                             <  1:176412/107‑1 (MQ=255)                          ggCAGACAAAGCAATATTATGGGCACAAATATCAGCCAGTACTAAAGAGGGACGAAAAGCGTGTTCACTAAGCTATTTTGCCTGTAAGGCCGCAGAGGCTAAACTCGGACTGGCATATATGGCTGCGAATGATAACAAGGAATTTCTCACATCCTTATCTAATATAATGAGGTACAAAATAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTg                                              <  1:249884/250‑1 (MQ=255)                              gACAAAGCAATATTATGGGCACTAATATCAGCCAGTACTAAAGAGGGACGAAAAGCGTGTTCACTAAGCTATTTTGCCTGTAAGGCCGCAGCGGCTAAAc                                                                                                                                                                                                >  1:77548/1‑100 (MQ=255)                                cAAAGCAATATTATGGGCACTAATATCAGCCAGTACTAAAGAGGGACGAAAAGCGTGTTCACTAAGCTATTTTGCCTGTAAGGCCGCAGAGGCTAAACTCGGACTGGCATATATGGCTGCGAATGATAACAAGGAATTTCTCACATCCTTATCTAATATAATGAGGTACAAAATAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATc                                       >  1:328105/1‑251 (MQ=255)                                               ngCACTAATATCAGCCAGTACTACAGAGGGACGAAAAGCGTGTTCACTAAGCTATTTTGCCTGTAAGGCCGCAGAGGCTAAACTCGGACTGGCATATATGGCTGCGAATGATAACAAGGAATTTCTCACATCCTTATCTAATATAATGAGGTACAAAATAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCat                                                   >  1:115466/2‑224 (MQ=255)                                                   ntAATATCAGCCAGTACTAAAGAGGGACGAAAAGCGTGTTCACTAAGCTCTTTTGCCTGTACGGCCGCAGCGGCTAAACTCGGACTGGCATATATGGCTGCGAATGATAACACGGAATTTCTCACATCCTTATCTACTATAATGAGGTACAAAATAGATGCTGAAATTTCCGAATCCTATACCTGCTATTTATTACGTAAGGGAAAGGTCATAAGACCATATCTGAAAAATCTAAATCTACTcc                           >  1:266016/2‑244 (MQ=255)                                                     aaTATCAGCCAGTACTAAAGAGGGACGAAAAGCGTGTTCACTAAGCTATTTTGCCTGTAAGGCCGCAGAGGCTAAACGCGGACTGGCATATATGGCTGCGAATGATAACAAGGAATTGCTCACATCCTTATCTAATATAATGAGGTACAAAAGAGATGCTGGACTTTCCGAATCCGATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGc                    <  1:294931/249‑1 (MQ=255)                                                         tCAGCCAGTACTAAAGAGGGACGAAAAGCGTGTTCACTAAGCTATTTTGCCTGTAAGGCCGCAGAGGCTAAACTCGGACTGGCATATATGGCTGTGCATGATAACAAGGAATTTCTCACATCCTTATCTAATATAATGAGGTACAAAACAGATGCTGGa                                                                                                          >  1:452558/1‑159 (MQ=255)                                                                gTACTAAAGAGGGACGAAAAGCGTGTTCACTAAGCTATTTTGCCTGTAAGGCCGCAGAGGCTAAACTCGGACTGGCATATATGGCTGCGAATGATAACAAGGAATTTCTCACATCCTTATCTAATATAATGAGGTACAAAATAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACACCAACTTGCCGCTGATTGCAtt       >  1:204913/1‑251 (MQ=255)                                                                            gACGAAAAGCGTGTTCACTAAGCTATTTTGCCTGTAAGGCCGCAGAGGCTAAACTCGGACTGGCATATATGGCTGCGAATGATaa                                                                                                                                                                 >  1:365484/1‑85 (MQ=255)                                                                                    ncNTGNNCACTNAGCTATTTTGCCTGTAAGGCCGCANAGGCTAAACTCGGACTGGCATATATTGCTGCNANTGATAACAAGGAATTTCTCACATCCTTATCTAATATAATGAGGTACAANATAGATGCTGNACTNNCCGANTNCNATACATGCTATTTATTAAGTAAGGGAAAGNTCATAAGACCATATCTGAAAAATNNAANNNCANNCCAAc                        <  2:328105/213‑1 (MQ=255)                                                                                                                  gCCGCAGAGGCTAAACTCTGACGTGCATATATGTCTGCGAATGATTACAAGGAATTTCTCTCATCCGTATCTAATATAATGTAGTATAAAATAGATGCTGAACTTGTCTAATCCCATACATGCTATGTATTAAGTAATTGAAAGAGCATATGACCATATCTGAAAAATTCAAATCCACTCCAACTTGCCGCTGATTGCATTGAAAc  <  1:301521/206‑1 (MQ=255)                                                                                                                              |                                                                                                                                                                                                  GACAATTTCAGAGAGAGTTAGATTTGGCAGACAAAGCAATATTATGGGCACTAATATCAGCCAGTACTAAAGAGGGACGAAAAGCGTGTTCACTAAGCTATTTTGCCTGTAAGGCCGCAGAGGCTGAACTCGGACTGGCATATATGGCTGCGAATGATAACAAGGAATTTCTCACATCCTTATCTAATATAATGAGGTACAAAATAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAAC  >  NC_000913/3769322‑3769640

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