Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_vch_tpiA	4,604,746	C→A	A186E (GCA→GAA)	bglJ →	bgl operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_vch_tpiA	4,604,746	0	C	A	100.0%	64.8 / NA	23	A186E (GCA→GAA)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (11/12);  total (11/12)

GACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTT  >  NC_000913_3_vch_tpiA/4604608‑4604890                                                                                                                                            |                                                                                                                                                 gacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGAACATaa                                                                                                                                            <  1:307190/145‑1 (MQ=255) gacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGAACATaa                                                                                                                                            >  2:307190/1‑145 (MQ=255)   cTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAACCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAACTCGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGAACATAAATTTaa                                                                                                                                      >  1:349497/1‑149 (MQ=255)      aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGAACAt                                                                                                                                              <  1:38472/138‑1 (MQ=255)      aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGAACAt                                                                                                                                              >  2:38472/1‑138 (MQ=255)      aaCAATCAGTGGTACACTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGACCTTTATGTCGCATGGCTACTCAATTACACAAAGTGCCGAGCAGCTTTAACGCAA‑TATCAAAACGATCGGTGAACATAAATTTAAttt                                                                                                                                   <  2:173568/149‑3 (MQ=255)           tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGAACATAAATTTAATGTGATGt                                                                                                                              <  2:438259/149‑1 (MQ=255)               tGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCACGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGaa                                                                                                                          >  1:325375/1‑149 (MQ=255)                     aTTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTCAACACAACCATC‑CAATGATCTGGGAACATCATTTTAATGTTATGTCGAAACTggg                                                                                                                    >  2:194947/1‑149 (MQ=255)                                      gTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGTTTAAACGCAA‑TATCAAAACGATCAGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAg                                                                                                   <  1:14343/149‑1 (MQ=255)                                             ccGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAAAGCAA‑TATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                            >  1:271987/1‑149 (MQ=255)                                                                 gcgaTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGTTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGt                                                                        <  1:533479/149‑1 (MQ=255)                                                                    cTTTATGTCGTGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat                                                                     <  2:268774/149‑1 (MQ=255)                                                                       tATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGCCCGCAGATATTCTGTtatgta                                                                  >  2:86183/1‑149 (MQ=255)                                                                                         aaTGACACAAATTGCCGAGCAGATTATACGCAG‑CACCAACACCGTCCGTGAACATAAATTTAATGTGATGTCGCAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGACACAGATATTCTGTTATGTGTGCAGCATGGCGAGACaa                                                >  2:5214/1‑149 (MQ=255)                                                                                           tGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAg                                               <  2:272334/148‑1 (MQ=255)                                                                                           tGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAg                                               >  1:272334/1‑148 (MQ=255)                                                                                             tcacaAATTGCGGAGAGGCTTAAACGCAA‑TATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTTTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTaa                                            <  2:305371/148‑1 (MQ=255)                                                                                                caatgtGCCGAGCAGCTTAAATGTAA‑TATCAAATCGTTTCGTGAACATGAATTTAATGTGATGTCGATACTGGGCGTCAGTTCTGACGTAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgt                                         <  1:244319/144‑1 (MQ=255)                                                                                                     tGCCGGGCAGCTTAAAGGCAA‑TAGCAAAACGATCCGTGAACATAAATTTAAGGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGc                                    <  2:60749/149‑1 (MQ=255)                                                                                                           gcagcTTAAACGCAA‑TATCAAAACGATCCGTGAACATAAATTTAATGTGATGTAGAAAATGGACGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGGAGATATTCAGTTATGTATGCGGCATTACGAAACAAGTAATGTGTTGCATACCt                              >  1:41046/1‑149 (MQ=255)                                                                                                                 tAAACGCAA‑TATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAAt                        >  2:420440/1‑149 (MQ=255)                                                                                                                                        cGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACtt  <  1:420440/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                 GACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAA‑TATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTT  >  NC_000913_3_vch_tpiA/4604608‑4604890

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