Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,707,947	C→A	intergenic (‑242/+669)	dppA ← / ← proK	dipeptide ABC transporter periplasmic binding protein/tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,707,947	0	C	A	100.0%	46.3 / NA	14	intergenic (‑242/+669)	dppA/proK	dipeptide ABC transporter periplasmic binding protein/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (8/6);  total (8/6)

CGCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCA  >  NC_000913/3707701‑3708151                                                                                                                                                                                                                                                       |                                                                                                                                                                                                             cGCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGCTGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACAGGGATTTACACCAACGGTGAGAATCcacaaaca                                                                                                                                                                                                           >  1:275359/1‑250 (MQ=255)             cTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGCGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGAGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTaaaaaa                                                                                                                                                                                               >  1:115506/1‑250 (MQ=255)                    gctgtttgaagttgttTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCaa                                                                                                                                                                                        >  2:216667/11‑250 (MQ=255)                         ttgaagttgttTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGc                                                                                                                                                                                   <  1:216667/245‑1 (MQ=255)                                      aaCCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCaga                                                                                                                                                                      <  1:8180/250‑1 (MQ=255)                                                 gTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCt                                                                                                                                                          >  1:243836/1‑251 (MQ=255)                                                                                tttgCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTaa                                                                                                                            >  1:118289/1‑250 (MQ=255)                                                                                tttgCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTaa                                                                                                                            <  2:118289/250‑1 (MQ=255)                                                                                                      aTGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGAtt                                                                                                      >  1:153714/1‑250 (MQ=255)                                                                                                             ggggATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACCCCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGTCCAGATCAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTa                                                                                              <  1:266518/251‑1 (MQ=255)                                                                                                                               aaCATCAGTGTAACAATGCGGATCCTCCATAAAATGCCCCTAGTGACCTACAATCTGTCAAAAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGACTCAACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGAAACATTTGGTGATTAAGTTTTATGCATTTAATGaaaaaa                                                                             >  2:246624/1‑250 (MQ=255)                                                                                                                                                                     ccTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCAt                                      <  1:169640/251‑1 (MQ=255)                                                                                                                                                                                                  gtgAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGcc         >  2:185018/1‑251 (MQ=255)                                                                                                                                                                                                   tgAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGcc         >  2:110419/1‑250 (MQ=255)                                                                                                                                                                                                         cGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAATAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATGAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGgcca  <  1:266813/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                             CGCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCA  >  NC_000913/3707701‑3708151

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