Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	4,604,632	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_pae_tpiA	4,604,632	0	C	A	100.0%	85.5 / NA	27	A186E (GCA→GAA)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (14/13);  total (14/13)

GACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_pae_tpiA/4604494‑4604777                                                                                                                                           |                                                                                                                                                  gacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGAGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGGAGCTTAAACGAAATATCAAAACGATCCGTGAACATAAAttt                                                                                                                                         >  1:457265/1‑149 (MQ=255)         aaTCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATCTTGCGCTTTATGTCGCGTAGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCAGTGAACATAAATTTAATGTGAt                                                                                                                                 >  2:426517/1‑149 (MQ=255)         aaTCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGAt                                                                                                                                 >  2:224648/1‑149 (MQ=255)                  tACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACt                                                                                                                        >  2:241433/1‑149 (MQ=255)                   aCATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTg                                                                                                                       <  2:376338/149‑1 (MQ=255)                               gCCGGACGTTAAGACCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCt                                                                                                           <  1:451807/149‑1 (MQ=255)                                      gTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGGCGAGCAGCTTAAACGCAATATCAAAACGATCAGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAg                                                                                                    >  2:588729/1‑149 (MQ=255)                                       ttAAGGCCGACGGAGAGAGCAATCTTGCGCTTTAGGTCGCGTGGCTACTCAGTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAGTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAgg                                                                                                   <  2:386504/149‑1 (MQ=255)                                            cccGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgt                                                                                              >  1:379810/1‑149 (MQ=255)                                              cGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                                                                                            <  1:549057/149‑1 (MQ=255)                                              cGACGGAGAGAGAAATATTGCGCTTCATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                                                                                            <  1:157380/149‑1 (MQ=255)                                                  ggAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCCGCGTTGTTGGAgg                                                                                        >  1:185101/1‑149 (MQ=255)                                                      agagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGc                                                                                    >  1:374806/1‑149 (MQ=255)                                                            atatTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGatat                                                                              <  1:224648/149‑1 (MQ=255)                                                              atTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTc                                                                            <  2:379810/149‑1 (MQ=255)                                                                   gcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGGAATATCAAAACGATCCGTGAACAGAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTta                                                                       <  1:347571/149‑1 (MQ=255)                                                                     tttATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCGGGGTTGTTGGAGGCCGCAGATATTCTGTtatg                                                                     >  1:36038/1‑149 (MQ=255)                                                                      tgaTGTTGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatgt                                                                    <  2:81420/147‑1 (MQ=255)                                                                           tCGCGTGGCTACTCCATGACACAAATTGCCGAGCAGCTTAAACGCAATATAAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTGGGAGGCCGCGGATATTCTGTTATGTATGCg                                                               >  2:229697/1‑149 (MQ=255)                                                                                       tCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGAGATTCTGTTATGTATGCGGCATTGCTAAAc                                                   <  1:588729/149‑1 (MQ=255)                                                                                                 aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCGGCGGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgtg                                         >  1:50059/1‑149 (MQ=255)                                                                                                        cGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATc                                  <  1:68651/149‑1 (MQ=255)                                                                                                        agAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATCTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGATGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATc                                  <  1:364809/148‑1 (MQ=255)                                                                                                                   aaCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATcc                       >  1:376914/1‑149 (MQ=255)                                                                                                                               aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg           >  2:358358/1‑149 (MQ=255)                                                                                                                               aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg           <  1:241433/149‑1 (MQ=255)                                                                                                                                        gTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTc  >  1:85871/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                  GACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_pae_tpiA/4604494‑4604777

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