Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_pgi	4,153,200	G→A	intergenic (‑102/+496)	ppc ← / ← argE	phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_pgi	4,153,200	0	G	A	100.0%	90.8 / NA	29	intergenic (‑102/+496)	ppc/argE	phosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/13);  total (16/13)

ATTACTACGCAATGCGGAATATTGTTCGTTCATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGT  >  NC_000913_3_pae_pgi/4153066‑4153334                                                                                                                                       |                                                                                                                                       aTTACTACGCAATGCGGAATATTGTTCGTTCATATTACCCCAGACACCCCATCTTTTCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCgg                                                                                                                          <  2:140204/149‑1 (MQ=255)    actacGCAATGCGGAATATTGTTCGTTCATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGa                                                                                                                       <  1:123242/149‑1 (MQ=255)      tacGCAATGCGGAATATTGTTCGTTCATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACg                                                                                                                     <  1:178259/149‑1 (MQ=255)                ggAATATTGTTCGTTCATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAAtgctg                                                                                                           <  1:359524/149‑1 (MQ=255)                       tGTTCGTTCATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATc                                                                                                    <  2:25314/149‑1 (MQ=255)                         ttcgttcATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGAc                                                                                                                      <  1:100086/129‑1 (MQ=255)                         ttcgttcATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGAc                                                                                                                      >  2:100086/1‑129 (MQ=255)                           cgttcATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCtt                                                                                                <  1:105208/149‑1 (MQ=255)                             ttcATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCa                                                                                              >  1:116278/1‑149 (MQ=255)                                                    tCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAAt                                                                       <  1:476657/149‑1 (MQ=255)                                                           gTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTCCAGCAAACGAATAAATAGCAGGAATTTACGTc                                                                >  2:21633/1‑149 (MQ=255)                                                            tttGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGCAATCGCTTCAGCAAACGAATAAATAGCCGGAATTTCCGTCa                                                               >  1:198068/1‑149 (MQ=255)                                                            tttGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGACTTTACGTCa                                                               >  2:7344/1‑149 (MQ=255)                                                            tttGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCa                                                               >  1:52832/1‑149 (MQ=255)                                                              tGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCAtt                                                             <  1:254173/149‑1 (MQ=255)                                                                    cccTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAAtgctgc                                                                                                          >  1:462104/1‑98 (MQ=255)                                                                    cccTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAAtgctgc                                                                                                          <  2:462104/98‑1 (MQ=255)                                                                      cTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTCCGTCATTAAATTAac                                                     >  1:158984/1‑149 (MQ=255)                                                                      cTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCac                                                     >  2:423711/1‑149 (MQ=255)                                                                              ttCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACCCttt                                             >  2:433694/1‑149 (MQ=255)                                                                                        aTTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTTAATAAACGt                                   >  2:168274/1‑149 (MQ=255)                                                                                                          gggCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTaa                 <  1:384134/149‑1 (MQ=255)                                                                                                           ggCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCgg                                                                                                                          >  1:386405/1‑43 (MQ=255)                                                                                                           ggCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCgg                                                                                                                          <  2:386405/43‑1 (MQ=255)                                                                                                           ggCTTTGCTTTTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTaaa                <  2:116278/149‑1 (MQ=255)                                                                                                                tgctttTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGAGGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATc           >  1:253347/1‑149 (MQ=255)                                                                                                                tgctttTCGTCGTCTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATc           >  1:62577/1‑149 (MQ=255)                                                                                                                   ttttCGTCGTCTTTTATAAAACCACGTAAAAGAGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTACATAAGCGTAACTTATGGAAATGTTAAAAAATCGcc        >  1:369939/1‑149 (MQ=255)                                                                                                                         tcgtcTTTTATAAAACCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGt  >  1:27888/1‑149 (MQ=255)                                                                                                                                       |                                                                                                                                       ATTACTACGCAATGCGGAATATTGTTCGTTCATATTACCCCAGACACCCCATCTTATCGTTTGATAGCCCTGTATCCTTCACGTCGCATTGGCGCGAATATGCTCGGGCTTTGCTTTTCGTCGTCTTTTATAAAGCCACGTAAAAGCGGTGACGTCAAATGCTGCGAAATCGCTTCAGCAAACGAATAAATAGCAGGAATTTACGTCATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGT  >  NC_000913_3_pae_pgi/4153066‑4153334

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