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

GACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTT  >  NC_000913_3_pae_tpiA/4604490‑4604771                                                                                                                                               |                                                                                                                                            gaccgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATaa                                                                                                                                       >  1:183210/1‑149 (MQ=255)    cgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAGTATCAAAACGATCCGTGAACATaaagt                                                                                                                                    >  1:143326/1‑147 (MQ=255)     gacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAttt                                                                                                                                   >  1:132356/1‑149 (MQ=255)     gacTGAACAATCAGTGGTACATTAACCAAACCCGGACGTTAAGACCGACGGACAGACAAATATTGCGCTTTATGTCGCGTGGCTCCTCAATGACACAAATTCCCCAGCAGCTTAAGCAGAATATCCAACCGCTCAGAGAACATAAAttt                                                                                                                                   >  2:464522/1‑149 (MQ=255)        tGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCCAAACGATCCGTGAACATAAATTTAAt                                                                                                                                <  2:136176/149‑1 (MQ=255)         gAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtg                                                                                                                               >  1:399741/1‑149 (MQ=255)           aCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtgtg                                                                                                                             <  1:353308/149‑1 (MQ=255)            cAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATGCGTGAACATAAATTTAAtatta                                                                                                                            >  2:74159/1‑145 (MQ=255)             aaTCAGTGGTACATTAACCAAAGCCGGACGTTAAGGTAGAGGGAGAGAGAAGTGTTGCGCTTTATGTCGCGTGGCTACTCAATGACAGAAATTGCTGCGCAGCTTAAACGCAATATCGACACGATCAGTGAAAATAAATTTAATGTGAt                                                                                                                           >  2:96798/1‑149 (MQ=255)                        cATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTgg                                                                                                                <  1:54593/149‑1 (MQ=255)                         aTTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTggg                                                                                                               >  2:518025/1‑149 (MQ=255)                          ttAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                              <  2:151011/149‑1 (MQ=255)                             aCCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTgcgctgc                                                                                                           >  2:228921/1‑146 (MQ=255)                                aaaGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGt                                                                                                        >  1:46348/1‑149 (MQ=255)                                 aaGCCGGACTTTAATCCTGACGGAGAGAAAAATTTTGTGCTTAATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCATCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                                                                                       <  2:246998/149‑1 (MQ=255)                                   gCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATTAAGTTAATGTGATGTCGAAACTGGGCGTCAGTTCt                                                                                                     <  2:122160/149‑1 (MQ=255)                                     cGGACGTTAATCCCGACGGAGAGAGAAATACTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                                                                                   <  1:52686/149‑1 (MQ=255)                                     cGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                                                                                   >  1:207521/1‑149 (MQ=255)                                          gTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAg                                                                                              <  1:228921/149‑1 (MQ=255)                                                            agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCGGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAg                                                                            >  1:453609/1‑149 (MQ=255)                                                                 tatTGCGCTTTATGTCGCGAGGCTAAGCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATAtt                                                                       <  1:529536/149‑1 (MQ=255)                                                                     gcgcTTTATGGCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAGACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGt                                                                   <  2:516591/149‑1 (MQ=255)                                                                                    gtgCTACTCAATGACACAAATTGCCTAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCAtt                                                    <  2:77898/147‑1 (MQ=255)                                                                                         aCTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaa                                               <  2:164294/149‑1 (MQ=255)                                                                                             aaTGACACAAATTGCCGAGCAGCTTAAACGAAATATCAAAAGGATCCGTGAACAGTCAATCAAGGAGATGTCGAAACTGGGAGTCAGTTCTGACGCAGGGGTGTTGGACGCCGCAGCTATTCTGTTATGTATGCGGCATTGCGCAACac                                           >  2:492006/1‑148 (MQ=255)                                                                                                 acacAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTaa                                       <  2:453609/149‑1 (MQ=255)                                                                                                 acaaAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTaa                                       <  2:160079/149‑1 (MQ=255)                                                                                                  cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatg                                                               <  1:243237/124‑1 (MQ=255)                                                                                                  cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatg                                                               >  2:243237/1‑124 (MQ=255)                                                                                                  cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGGGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGTGAAACAAGTAAt                                      <  1:96798/149‑1 (MQ=255)                                                                                                    caAATTGCCGAGCAGCTTAAACCCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCAGACGCAGGGTTGTTGGAGGCCGCAGATATGCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgt                                    <  1:135138/149‑1 (MQ=255)                                                                                                            cGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATc                            >  1:359796/1‑149 (MQ=255)                                                                                                               ggagtttAAAAGCAATATCAAAACGTTCCGTGAACATAAATTTCATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCt                         <  2:417252/144‑1 (MQ=255)                                                                                                                     tAAACGCAATATCAAAAAGATCCGGGAACCTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAAt                   <  2:185527/149‑1 (MQ=255)                                                                                                                       aaCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATcc                           <  1:384940/139‑1 (MQ=255)                                                                                                                       aaCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATcc                           >  2:384940/1‑139 (MQ=255)                                                                                                                          gCATTATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAATCAAGTAATGTGTTGCATCACTATTAATCCGCa              >  1:114496/1‑149 (MQ=255)                                                                                                                          gCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaa                                               >  2:143750/1‑116 (MQ=255)                                                                                                                          gCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaa                                               <  1:143750/116‑1 (MQ=255)                                                                                                                                tCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGc        >  1:568010/1‑149 (MQ=255)                                                                                                                                   aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGGTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAACTAATGTGTTGCATCCCTATTAATCCGCATGAGTCCgg     >  1:459331/1‑149 (MQ=255)                                                                                                                                    aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCggg    >  1:377622/1‑149 (MQ=255)                                                                                                                                      cGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGtt  <  2:399741/149‑1 (MQ=255)                                                                                                                                      cGATCCGCGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACTCAGGGTTGTTGGAGGCCTCAGATAATCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGtt  <  1:138077/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                            GACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTT  >  NC_000913_3_pae_tpiA/4604490‑4604771

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