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

CCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_vch_tpiA/4604606‑4604891                                                                                                                                             |                                                                                                                                                  ccgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAt                                                                                                                                           >  2:474083/1‑149 (MQ=255)   gacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGATTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAttt                                                                                                                                         >  2:444506/1‑149 (MQ=255)          tccgatctGGCACATTAACCAAAGCCGTACGTTAAGCCAGACGGAGAGAGAACTATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTATTTTTAAAATCAAAACGATCCGTGAACATAAATTTAATGTGa                                                                                                                                  <  1:577261/142‑1 (MQ=255)           aaTCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAAGGTAATATCAAAACGATCCGTGAACATAAATTTAAtgtggt                                                                                                                                 <  2:59322/149‑3 (MQ=255)                       aTTAACCAAAGCCGGACGTTCAGCCCGACGGAGAGAGAAATAAAACGCATTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTggg                                                                                                                     >  1:284009/1‑149 (MQ=255)                       aTTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATg                                                                                                                                >  1:317248/1‑138 (MQ=255)                       aTTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATg                                                                                                                                <  2:317248/138‑1 (MQ=255)                              aaaGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGt                                                                                                              >  1:238771/1‑149 (MQ=255)                                    ggACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGAc                                                                                                        >  1:408850/1‑149 (MQ=255)                                         ttAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAgg                                                                                                   >  2:436887/1‑149 (MQ=255)                                         ttAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAgg                                                                                                   >  2:374770/1‑149 (MQ=255)                                          tAAGCCCGACGGAGAGAGAGAGATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGCTCCGTGAACTGAAAAGTAAGGTTATGTCAAAAATGGGGGTCAGTTCTGACGCAggg                                                                                                  >  2:111352/1‑149 (MQ=255)                                               ccGACGGAGAGAGAAATATTGCGATTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                             <  1:444506/149‑1 (MQ=255)                                                  aCGGAGAGAGAAATATTGCGCTTTAGGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                          <  2:553026/149‑1 (MQ=255)                                                           gaAAGAGTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGa                                                                                 <  2:238771/149‑1 (MQ=255)                                                               tatTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGTAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGACGCAGACAtt                                                                             >  2:340001/1‑149 (MQ=255)                                                               tatTGCGCTGTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATAtt                                                                             <  2:157278/149‑1 (MQ=255)                                                                    cgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGtt                                                                        <  1:436887/149‑1 (MQ=255)                                                                     gcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTta                                                                       >  1:504658/1‑149 (MQ=255)                                                                                       aCTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATAAAAACGATAGGTGAACGTAAATTTTATGAGGTGTCGAAACGGCGCGTCGGTTATGATGTAGGGTGGTTGGACGCCGATGGTCTGCTGTTGTGTATGCGTCATTGCGaa                                                     >  2:19608/1‑149 (MQ=255)                                                                                           aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                                 >  1:461671/1‑149 (MQ=255)                                                                                               gctcaAATCGCCTAGCAGCTTAAACGCAATATCAAAACGATCCTTGAACATAAACTTAATGTGATGCCGAAACTGGGTGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGACACAAGTaa                                             <  2:12282/146‑1 (MQ=255)                                                                                                 acaAATTGCCGAGCAGCTTAAACGCCATATCAAAACGATCCGAAAACAAAAATATAATGTGATGTAGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGaaagg                                           >  1:153549/1‑147 (MQ=255)                                                                                                 acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAAGGTTGTGGGATGCCGTCGATATTATGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                           >  1:353427/1‑149 (MQ=255)                                                                                                   aaaTTGCCGAGCAGCTTAACCGCAATATCAAACCGATCCGTGAACATAAATATAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCGGCAGCACAAAAGTGATTTGTGCGGCATTGCGAAACAAGTAAtgtg                                         >  1:443558/1‑149 (MQ=255)                                                                                                   aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgtg                                         >  1:286980/1‑149 (MQ=255)                                                                                                          cGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTTATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATc                                  <  2:461671/149‑1 (MQ=255)                                                                                                                  ttAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACCAGGGGTTAGTTCTGACGCAGGGTTGTCGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAACCAAGTAATGTGTTGCATCCCTATTaa                          >  2:497228/1‑149 (MQ=255)                                                                                                                       cGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAAATGGGCGTCAGTTCTGACGCAAGGTTCTTGGAGGCCGGAGATATGCTGTTATGTATTCGGCATTGCGAAACAGATAATGTGTTGCATCCCTATTAATCCGc                     >  2:17828/1‑149 (MQ=255)                                                                                                                        gCAAGATCAAAACCATCAGTGAAAATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGTGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGTGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCa                    <  2:245766/149‑1 (MQ=255)                                                                                                                         cAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCat                   <  1:374770/149‑1 (MQ=255)                                                                                                                                 caaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatgt                                                                    >  1:107768/2‑92 (MQ=255)                                                                                                                                 aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatgt                                                                    <  2:107768/92‑1 (MQ=255)                                                                                                                                    cGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGtt        <  1:384265/149‑1 (MQ=255)                                                                                                                                      aTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCg                                                               >  1:389642/1‑92 (MQ=255)                                                                                                                                      aTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCg                                                               <  2:389642/92‑1 (MQ=255)                                                                                                                                      aTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGCTGCCGGGTTTa      >  2:130581/1‑149 (MQ=255)                                                                                                                                          gTGAACATAAATTTAATGTGATGTCGAAACTGGGCATCAGTACTAACGCAGGGATGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTc  >  1:438338/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                  CCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_vch_tpiA/4604606‑4604891

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