Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	349,673	G→C	intergenic (+101/‑339)	prpB → / → prpC	2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	349,673	0	G	C	100.0%	45.3 / NA	24	intergenic (+101/‑339)	prpB/prpC	2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (11/13);  total (11/13)

ACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTT  >  NC_000913_3_pae_tpiA/349526‑349722                                                                                                                                                    |                                                  aCGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCAt                                                 <  3:93161/150‑1 (MQ=25)  cGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCAta                                                <  1:964067/150‑2 (MQ=21)  cGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCAta                                                <  1:85798/150‑2 (MQ=21)    aaGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGc                                              >  1:816141/1‑150 (MQ=21)       agaAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                           <  1:696794/150‑1 (MQ=21)           gCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                                       <  1:768389/150‑1 (MQ=18)               gACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCCCTGAACTGTAGGTCGGATAAGGCGCCTAGCGTCGCATCCGa                                   >  3:344828/1‑150 (MQ=18)                    ccTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                              >  1:97103/1‑150 (MQ=17)                     cTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCg                             >  1:108405/1‑150 (MQ=17)                              cGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                    <  1:556095/150‑1 (MQ=17)                                   ccAGGTGAAATAAAAAACGCCCGTTGATTGTATGCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgcc               <  2:149192/150‑1 (MQ=14)                                     aGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGACCTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc             <  1:562879/150‑1 (MQ=14)                                     aGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc             >  2:771189/1‑150 (MQ=14)                                       gTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg           >  1:787505/1‑150 (MQ=14)                                       gTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg           >  2:558791/1‑150 (MQ=14)                                        tGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa          >  1:668024/1‑150 (MQ=14)                                           aaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACgc       <  2:787505/150‑2 (MQ=14)                                           aaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACgc       <  2:97103/150‑2 (MQ=11)                                           aaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACgc       <  1:558791/150‑2 (MQ=11)                                             tAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa     <  1:614889/150‑1 (MQ=11)                                             tAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa     >  2:764285/1‑150 (MQ=11)                                                aaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt  >  3:126788/1‑150 (MQ=11)                                                  aaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACgc       <  3:149995/143‑2 (MQ=11)                                                  aaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACgc       >  4:149995/1‑142 (MQ=11)                                                                                                                                                    |                                                  ACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTT  >  NC_000913_3_pae_tpiA/349526‑349722

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