Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,282,332	(A)5→6	coding (744/1293 nt)	yjcF ←	pentapeptide repeat‑containing protein YjcF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,282,332	1	.	A	100.0%	16.3 / NA	7	N248? (AAT→AAN)	yjcF	pentapeptide repeat‑containing protein YjcF
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (2/5);  total (2/5)

GAGTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTT  >  NC_000913/4282096‑4282520                                                                                                                                                                                                                                              |                                                                                                                                                                                             gagTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCaa                                                                                                                                                                                 <  1:210653/251‑1 (MQ=255)         aCATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAATGTCGGCTCGAAATAAATTAATATATGGCTCTATTACTTTATAAATATTTACCGCAGCCAGTTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTa                                                                                                                                                                         <  1:343621/251‑1 (MQ=255)           atatAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATGCAGGTTTACACGAGACACATCAGCATTTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCCAAATAAATTAATAGATAGCGCTATTAATTTATAAAGATTTGCCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTAtt                                                                                                                                                                       <  2:272827/251‑1 (MQ=255)                                                                                                              gtgaCAGATCAGCATGTTCAAAGTTAATTAAATCAAGGTTAGAGAAGCTAAGGTCGGCTATAAATAAATTAATATATGGCGCTATTACTTTAAAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATGTTTAATAGTGGTTCTATTCATTGTTGCATTGTAGAACTGAACATGCATAGAGGAACATGAATATAGAGTTCCGTCAGTCAGTTTTATCTGAGCGAAAtt                                                                    <  2:150510/249‑1 (MQ=255)                                                                                                                                                 aGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTAttt                                 >  1:303212/1‑251 (MQ=255)                                                                                                                                                               aaGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCg                   >  2:284918/1‑251 (MQ=255)                                                                                                                                                                                aatTAATATATGGCGATATTAATTTATAAAGATTTACCTTAGCAATGTATGCCATATAAAAAATTTGAATAATCAAGGAAGGTATTGTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAAGTTAGAGTTCCGTCAGTCATTTTTATCCGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTAGTTTTTCCATTATGGCGAAATTAAAAGTCggtgt  <  2:176714/251‑5 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                             GAGTCGGTACATATAACCATTTCGAGGAACGTATTTGTTAACCGCGTAAAAAAGAGTTTGCTATCAATTAAGTTTATATTCTGGAGGGTTGCTTTATTCAGGTTGACACGAGACAGATCAGCATGTTCAAAGTTAATTAAATCAAGTTTAGAGAAGCTAAGGTCGGCTCTAAATAAATTAATATATGGCGCTATTACTTTATAAAGATTTACCTCAGCCATGTATGCCATATAAAAA‑TTTGAATAATCAAGGAAGGTATTTTTAATATTGGCTCTATTCATTGTTGCATTGTAGAATTGAACATGCATAGCGGAACATGAATTTAGAGTTCCGTCAGTCATTTTTATCTGAGCGAAATTTGAGTCATCGAGAATACAATTATCAAAAAGTATTTTTTCCATTATGGCGAAATTAAAATTCGTTTT  >  NC_000913/4282096‑4282520

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