Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,153,441	G→A	intergenic (‑343/+255)	ppc ← / ← argE	phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,153,441	0	G	A	100.0%	78.9 / NA	23	intergenic (‑343/+255)	ppc/argE	phosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/14);  total (9/14)

AAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCC  >  NC_000913/4153207‑4153685                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     aaaaGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCAt                                                                                                                                                                                                                                       <  1:273848/250‑1 (MQ=255)           aCGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCa                                                                                                                                                                                                                            <  1:96948/251‑1 (MQ=255)                    gctgcGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTgtcggtc                                                                                                                                                                                                                   <  1:109715/251‑1 (MQ=255)                          aaaTCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATa                                                                                                                                                                                                              >  1:184092/1‑250 (MQ=255)                                       aaaCGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTAt                                                                                                                                                                                                <  1:305148/251‑1 (MQ=255)                                               aaataGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTGGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCg                                                                                                                                                                                        <  2:123126/251‑1 (MQ=255)                                                aataGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCg                                                                                                                                                                                        >  1:34127/1‑250 (MQ=255)                                                   aGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATTTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGGCGGTCGGATAAGATGCGCAAGTATCGCATCCGACat                                                                                                                                                                                    <  2:351485/251‑1 (MQ=255)                                                                      aaaTTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGttt                                                                                                                                                                 <  1:142942/251‑1 (MQ=255)                                                                                                    tATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGAGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGAGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc                                                                                                                                   >  2:348414/1‑251 (MQ=255)                                                                                                            tGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGc                                                                                                                                                      <  1:186840/224‑1 (MQ=255)                                                                                                            tGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGc                                                                                                                                                      >  2:186840/1‑224 (MQ=255)                                                                                                                  aaaaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACg                                                                                                                     >  2:82136/1‑251 (MQ=255)                                                                                                                              aGTAACACAAAAGGTGTGGGTCGGATAAGCTGCGCAAGTATCGCAACCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGt                                                                                                         <  2:78246/251‑1 (MQ=255)                                                                                                                                tAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGttt                                                                                                       >  1:351222/1‑251 (MQ=255)                                                                                                                                               aGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGGTGGGCATCGTTTGAGTCACTGTAGGt                                                                                         <  1:111454/250‑1 (MQ=255)                                                                                                                                                                  aGTAGCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAgcat                                                                     <  2:285632/251‑1 (MQ=255)                                                                                                                                                                                aatgattGTGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGAGTGATCCTACCTACAGGTGGGCATTGTTTGAGTCACTGTCGGTCGGATAAGAGGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGTATAAGATGCGCAAGCATCGCATCCGACatta                                                       <  2:306846/247‑1 (MQ=18)                                                                                                                                                                                   tattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTTGGTATCGTTAGAGTCACTGTAGTTCCGATAAGATGCGCTAGCATCGCGTCCggacatattg                                                    >  2:126954/1‑242 (MQ=17)                                                                                                                                                                                               ggAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGt                                        >  1:240385/1‑251 (MQ=255)                                                                                                                                                                                                               gccggtggcggcgCGATCGCCTTATCCGACCTACAGTTGGGAATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCAGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCgg                        <  2:249174/244‑1 (MQ=255)                                                                                                                                                                                                                    aTGCGGCGCGAGCGCCTTATCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATAGTGTGAGTCACTGTAGGTCGGATAAGATGCGCACGCATCGCATCCGACATTATTGCGGCCCTGGAGTTTGGCAACAGCGCCGGATgcg                   >  2:74064/1‑251 (MQ=255)                                                                                                                                                                                                                                      aTCCGACCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATcc  <  1:6463/250‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     AAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCC  >  NC_000913/4153207‑4153685

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