Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,296,381	+GC	intergenic (+587/+55)	gltP → / ← yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,380	1	.	C	100.0%	59.3 / NA	18	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (11/7);  total (11/7)
*	NC_000913	4,296,380	2	.	G	100.0%	58.9 / NA	18	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (11/7);  total (11/7)

CTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTA‑‑CCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCAGCGCC‑‑TTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATT‑CCCG‑CCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGTT  >  NC_000913/4296310‑4296666                                                                            ||                                                                                                                                                                                                                                                                                                                cTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGt                                                                                                                                                                   <  1:894854/202‑1 (MQ=255)  ttATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCGTTTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTCTTTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCAGCGCC‑‑TTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTaa                                                                                                                >  1:1194447/1‑252 (MQ=255)  ttATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCAGCGCCT‑TTTTGCTTGTTCTTCTCGATAAAAcc                                                                                                                               >  1:1004064/1‑236 (MQ=255)  ttATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCAGCGCCT‑TTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGt                                                                                                                    >  1:304455/1‑247 (MQ=255)  ttATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCAGCGCC‑‑TTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATT‑CCCG‑CCCAGTACTCGGAATAACCGGTTCGGGAAAtt                                                                >  1:894856/1‑296 (MQ=255)      cAGGCCTACGCCAGACAGCGCTAATAGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCGTTTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑‑TTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TACCAGCTCAGGTTCAACCAGTg                                                                                                                                                                  >  1:894858/1‑200 (MQ=18)              cGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCGTTTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTT‑A‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCCTTTCCACTA‑CAGCTCAGGTTCAACCAg                                                                                                                                                                    >  1:1035980/1‑190 (MQ=18)                        cgc‑AATAGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑ttttcaa                                                                                                                                                                                                                             >  1:304457/1‑124 (MQ=14)                        cgc‑AATAGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCAGCGCC‑‑TTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATT‑CCCG‑CCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAgg                                       >  1:304456/1‑299 (MQ=255)                               aGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCCTTATCTTAACCG‑TTGGTTAA‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTTAACCAGTGCAGCGCC‑‑TTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATT‑CCCG‑CCCAGTACTCGGAATAACCGGTTCGGGaa                                                                   <  1:1007827/266‑1 (MQ=255)                                gCCTGATTTAGCGTGATTTTGTTAGG‑TCGGATAAGGCGTTTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGCTTTGCCGCGTCTTTATCAGGCC‑TATCTTAACCGTTTGGTTAATTTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCAGCGCCCTTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTaa                                                                                                                >  1:1191064/1‑229 (MQ=255)                                  tgattttAGCGTGATTTTG‑TAGG‑TCGGATAAGGCGTTTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTT‑A‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCACTA‑CAGCTCAGGTTCAACCAGTGCAGCGCC‑‑TTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcatc                                                                                                         >  1:894861/5‑227 (MQ=37)                                      tttAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑‑TGGTTAA‑TTTTTC‑AAACTCTTCACACCCCGGTATCAAACCCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCAGCGCC‑‑TTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCCGCCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGcccccacgccg  <  1:894859/326‑9 (MQ=38)                                       ttAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑TTTTTCAAAACTCTTCACACCCCGGTATCAAACCCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCAGCGCC‑‑TTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCCGCCCCAGTACTCGGAATAACCGGTTCGGGaaaattg                                                              <  1:894860/267‑5 (MQ=37)                                        tAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCa                                                                                                                                                                >  1:304459/1‑167 (MQ=255)                                                ttttG‑TAGGTTCGGATAAGCTG‑TTGACGCCGCTATCCGACATCAATGACCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑‑TTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTc                                                                                                                                                                                <  1:304458/145‑1 (MQ=21)                                                      tAGG‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTA‑‑‑TTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTc                                                                                                                                                                                            <  1:994975/126‑1 (MQ=18)                                                         g‑TCGGATAAGGCG‑TTTACGCCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCAGCGCC‑‑TTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGtt                                                                                                                   <  1:894862/194‑1 (MQ=255)                                                                            ||                                                                                                                                                                                                                                                                                                                CTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTG‑TAGG‑TCGGATAAGGCG‑TTTA‑‑CCGC‑ATCCGACATCAATG‑CCTGATGCGACGC‑TTGCCGCGTC‑TTATCAGGCC‑TATCTTAACCG‑TTGGTTAA‑TTTTTC‑AAACTCTTCACA‑CCCGGTATCAAA‑CCC‑TTCCA‑TA‑CAGCTCAGGTTCAACCAGTGCAGCGCC‑‑TTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATT‑CCCG‑CCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGTT  >  NC_000913/4296310‑4296666

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