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

CTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_pae_tpiA/4604496‑4604777                                                                                                                                         |                                                                                                                                                  cTGAACAACCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTaa                                                                                                                                       <  1:95812/149‑1 (MQ=255)  taaaCAATCAGTGGTACATTAACCAAAGCCGGACGTGAAGCCCGACGGAGAGAGAAATACTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAt                                                                                                                                      <  1:16765/147‑1 (MQ=255)  tGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATGGCTGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAAATTAAATTTaag                                                                                                                                      >  2:281370/1‑148 (MQ=255)     aCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtgtg                                                                                                                                   >  1:332113/1‑149 (MQ=255)         tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGAAGGAGAGAGAAATATTGAGCTTTATGTCACGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATACGTGAACATAAATCTAATGTGATGt                                                                                                                               >  2:297647/1‑149 (MQ=255)         tCAGTGCTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATCTTGCGCTTTATGTCGCGTGGCTTCTCCATGACACAAATTGCCGAGCAACTTAAACGCAGTATCAAAACGATCCGTGAACATAAATTTAACGTGATGt                                                                                                                               >  2:297552/1‑149 (MQ=255)               gTACATTAACCAAAGCCGGACGTTAACCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACCCAAATTGCCGAGCAGCTTAAACGCAATCTCAAAAGACTCCGTGAAACTAAATTTATTATGATGTTGAAAc                                                                                                                         >  1:2432/1‑149 (MQ=255)                 aCATTAACCAAAGGCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTg                                                                                                                       <  1:266747/149‑1 (MQ=255)                        ccAAAGCCGGAAGTTAAGCCCGACGGAGAGAGAAATAGTGGGGTTTATGTCGCGTGGTGACTCAATGACACAAAGTGCCGGGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCa                                                                                                                <  1:146029/149‑1 (MQ=255)                                               cGGAGAGAGAAATATTGCGCTTTATGTCGGGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCCAACAAGGGTGTCAGTTCTGACGCAGGGTTGTTGGAg                                                                                         <  1:492361/149‑1 (MQ=255)                                                       gaAATATTGCGCTTTATGTCGCGTGGCTACACAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGa                                                                                 <  1:327651/149‑1 (MQ=255)                                                         aaTATTTCGCTTTATGTCGCGTGGTTATTCAATGAAAAAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGata                                                                               <  1:251832/149‑1 (MQ=255)                                                         aaTATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGata                                                                               <  1:211554/149‑1 (MQ=255)                                                               gcgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATTCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGt                                                                         >  2:314473/1‑149 (MQ=255)                                                                   tttATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatg                                                                     >  2:262369/1‑149 (MQ=255)                                                                    ttATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatgt                                                                    <  1:276573/149‑1 (MQ=255)                                                                             gTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCATGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                                                                                             >  2:169997/1‑99 (MQ=255)                                                                             gTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCATGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                                                                                             <  1:169997/99‑1 (MQ=255)                                                                                    cTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATTCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaaa                                                    <  1:314473/149‑1 (MQ=255)                                                                                     tCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGGCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCAATGCGAAAc                                                   >  1:440876/1‑149 (MQ=255)                                                                                       aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                                 >  1:177954/1‑149 (MQ=255)                                                                                            cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAt                                            >  2:194517/1‑149 (MQ=255)                                                                                              caAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgt                                          <  1:388478/149‑1 (MQ=255)                                                                                                             cTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTa                           >  2:441936/1‑149 (MQ=255)                                                                                                             cTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTa                           >  2:69053/1‑149 (MQ=255)                                                                                                                    gCAATATAAAAACGATCCGTGAACAGAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGCCGCAGGGTTGCTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCa                    <  2:2432/149‑1 (MQ=255)                                                                                                                     aaaTATCAAAACGTTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCat                   <  2:217946/148‑1 (MQ=255)                                                                                                                       atatCAAAACGACCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGTatga                 <  1:455691/149‑1 (MQ=255)                                                                                                                        tatCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCa                                    <  2:32267/129‑1 (MQ=255)                                                                                                                        tatCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCa                                    >  1:32267/1‑129 (MQ=255)                                                                                                                          tCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGc              >  2:88140/1‑149 (MQ=255)                                                                                                                                 gATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGttt       <  1:332095/149‑1 (MQ=255)                                                                                                                                      gTGAACATAAATTTAATGTAATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATACGCGGTTTACTTc  >  2:333832/1‑149 (MQ=255)                                                                                                                                         |                                                                                                                                                  CTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_pae_tpiA/4604496‑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: ‑