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%	123.1 / NA	38	A186E (GCA→GAA)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/22);  total (16/22)

CTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTAC  >  NC_000913_3_pae_tpiA/4604496‑4604774                                                                                                                                         |                                                                                                                                               cTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTaa                                                                                                                                    <  2:361623/149‑1 (MQ=255)  tGAACAATCAGTGGTACAGTAACCAAAGCCGGACGTTAATCCCGACGCAGAGAGAAATATTGCGCTTTATGTCGCGTGGATACTCAATGACACAAATTGCCGAGCGGCTTAAACGCAATATCATAACGATCCGTGAACATAAATTTAAt                                                                                                                                   <  1:519537/149‑1 (MQ=255)         tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGt                                                                                                                            >  2:455237/1‑149 (MQ=255)           aGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCg                                                                                                                          >  1:352393/1‑149 (MQ=255)           aGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAAAATAAATTTAATGTAATGTCg                                                                                                                          >  1:288594/1‑149 (MQ=255)            gTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGa                                                                                                                         >  1:197443/1‑149 (MQ=255)             tGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTCAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGaa                                                                                                                        >  1:98890/1‑149 (MQ=255)                    ttAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                                 <  1:150138/149‑1 (MQ=255)                        ccAAAGCCGGAGGTTAAGCCTGCTGTGGAGAGAATTATTGGGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATTAAGTTAATGTGATGTCGGAACTGGGCGTCa                                                                                                             <  2:315112/149‑1 (MQ=255)                                      tAAGCCCGACGAAGAGATATAGATTGTGCTTTATGGCGCGTTTTTACTCAATGACACAAATTGCCTAGCAGCTTAAACGAAATAGCAAATCGAGCCTTGAACATAAATTTATTGGGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                                                                                               <  2:328450/149‑1 (MQ=255)                                         gCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttg                                                                                            <  1:564149/149‑1 (MQ=255)                                             gACGGATGGAGAAATATTGCGCTTTATGTCGCTTGGCTACTCCATGACACAAATTGGCGAGCAGCTTAAACGCAATATCAAAACGATCCGGGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTgg                                                                                        <  2:237758/149‑1 (MQ=255)                                             gACGGAGAGAGAAATATTTCGCTTTACGGCGCGTGGCTACTCAATGACACAAAGTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTgg                                                                                        <  1:643283/149‑1 (MQ=255)                                             gACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTgg                                                                                        >  2:352161/1‑149 (MQ=255)                                                ggAGAGAGAAATATTGCGCTTTATGGCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                                                                                     <  2:173532/149‑1 (MQ=255)                                                      agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGCTCTGACGCAGGGGTGTTGGAGGCCGCAg                                                                               >  1:261960/1‑149 (MQ=255)                                                        aaaTATTGCGCTTTATGTCGCGTGTCTACGCAAGGACACAAATTGCCGAGCAGCCTCAACGCAACATCAAAACGTTCCGATAGCATAATTTTAATGTGATGTCGAACCTGGTCGTCAGTTCTGACGCAGGGTTGTTGGTGGCCGCAGat                                                                             <  2:248319/149‑1 (MQ=255)                                                              tGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTg                                                                       <  1:341754/149‑1 (MQ=255)                                                                cgttttATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGtt                                                                     <  2:491554/146‑1 (MQ=255)                                                                        gTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGc                                                             <  1:462537/149‑1 (MQ=255)                                                                              tGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCAtt                                                       <  1:216948/149‑1 (MQ=255)                                                                                    cTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaaa                                                 >  2:421794/1‑149 (MQ=255)                                                                                    cTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaaa                                                 <  1:377089/149‑1 (MQ=255)                                                                                       aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAATTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                              <  1:372978/149‑1 (MQ=255)                                                                                             acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                        >  1:489470/1‑149 (MQ=255)                                                                                                aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATgtgt                                     <  1:207893/149‑1 (MQ=255)                                                                                                            gcTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCtttt                         >  2:678110/1‑149 (MQ=255)                                                                                                             cTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTa                        >  1:140755/1‑149 (MQ=255)                                                                                                               tAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAAt                      >  1:54904/1‑149 (MQ=255)                                                                                                                aaaCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATc                     >  1:203090/1‑149 (MQ=255)                                                                                                                       atatCAAAACGATCCGTGAACATAAATTTAATGTGATGCCGAAACTGCGCTTCAGTTCTGACGCAGGGTTGTTGGAGGCCGAAGATATTCTTTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatga              <  2:379834/149‑1 (MQ=255)                                                                                                                         atCAAAACGCTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatgatg            <  1:455319/149‑1 (MQ=255)                                                                                                                            aaaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCCAAACAAGTAATGTGTTGCATCCCTATTAATCCGCCTGATGCCg         >  2:36515/1‑149 (MQ=255)                                                                                                                             tcaCGCTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg        <  2:140755/147‑1 (MQ=255)                                                                                                                             aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGTGGCCGCAGATATTCTGTTATGTATGCGGCCTTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatgcttccag        >  2:262278/1‑142 (MQ=255)                                                                                                                             aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg        >  2:219763/1‑149 (MQ=255)                                                                                                                                  aTCCGGGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGCAATGTGTTGCATCCCTCTTCTGCCGCATGATGCCGGGTTTa   <  2:179215/149‑1 (MQ=255)                                                                                                                                   tcggTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTAc  <  2:197443/146‑1 (MQ=255)                                                                                                                                         |                                                                                                                                               CTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTAC  >  NC_000913_3_pae_tpiA/4604496‑4604774

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: ‑