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

CCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCG  >  NC_000913_3_vch_tpiA/4604606‑4604881                                                                                                                                             |                                                                                                                                        ccgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTGAAACGCAATATCAAAACGATCCGTGAACATAAAt                                                                                                                                 >  1:151038/1‑149 (MQ=255)    acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATaa                                                                                                                                   <  1:143896/144‑1 (MQ=255)    acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATaa                                                                                                                                   >  2:143896/1‑144 (MQ=255)      tGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAt                                                                                                                            <  2:267431/149‑1 (MQ=255)       gAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtg                                                                                                                           >  1:321699/1‑149 (MQ=255)       gAACAATCAGTGGTACAGTAACCAAAGCCGGACGTTAAGCCCGACGGATAGAGAAATATTGCGCTTTACGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtg                                                                                                                           <  2:284292/149‑1 (MQ=255)        aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtgt                                                                                                                          >  1:429334/1‑149 (MQ=255)        aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATACGTGAACATAAATTTAAggt                                                                                                                          >  2:162019/1‑149 (MQ=255)          cAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGAGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGa                                                                                                                        >  1:76094/1‑149 (MQ=255)                  ggTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGTTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGaaa                                                                                                                <  1:339853/149‑1 (MQ=255)                       aTTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTggg                                                                                                           >  2:461015/1‑149 (MQ=255)                          aaCCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTg                                                                                                             <  2:237432/144‑1 (MQ=255)                          aaCCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTg                                                                                                             >  1:237432/1‑144 (MQ=255)                            gcAAAGCCGGGCGTGAAGCACGACGGCGAGAGCAATGTTGCGCGTTATGTCGCGTGTCTACTCAATGAACCAAATTGGTGAGCAGCTTAAACGCAATATCAAAACGAGCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCa                                                                                                      <  2:115558/148‑1 (MQ=255)                                    ggACGTTAAGCCCGACGGACAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGAc                                                                                              <  1:404427/149‑1 (MQ=255)                                           aaGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgt                                                                                       <  1:102389/149‑1 (MQ=255)                                               ccGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTCAACGCAATATCAACACGAACCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGGCAGTTCTGACGCAGGgttgtt                                                                                   >  2:195465/1‑149 (MQ=255)                                                 gACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGCGGTCAGTTCTCACGCAGGgttgattg                                                                                 >  2:33657/1‑145 (MQ=255)                                                 gACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACAAAAATTGCCGAGCCGCTTAAACGCAATATCAAAACGATCCGTGAAAATAAATTTAATGTGATGGCGAAATTGGGCGTCAGTTCTTCCGCAGAGatgttga                                                                                 >  2:33315/1‑148 (MQ=255)                                                            aaaTATTGCGCTTTATGTCGCGTGGCTACTCACTGACCCAAATTGCCGAGCAGCTTAACCACAATATCACACCGATCCGTGAAAATAAATATAATGTGATGTCGATGCTGGGCGTCAGTTTTGACTCAGGGTTGTTGGAGGCCGCAGag                                                                      >  2:303989/1‑148 (MQ=255)                                                               gatTGCGCTTTATGTCGCGTGGCTACTCACTGAGACAAATGCCCGAGCAGCTTAAACGCAATATCAAAACGATCCGCGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATAtt                                                                   <  2:413754/148‑1 (MQ=255)                                                                           tGTCGCGTGGCTACTCAATGACACAAATTGCCGTGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatgtatg                                                       <  1:350325/149‑1 (MQ=255)                                                                             tCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCg                                                     <  1:152249/149‑1 (MQ=255)                                                                                        ttCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaaa                                          <  1:347208/148‑1 (MQ=255)                                                                                        cTCAATGACATAAATGGCCGAGCAGCTGAAACGCAAGATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCTCAGATATTCTGTTATGTATGCGGCATTGCGaaa                                          <  2:117754/149‑1 (MQ=255)                                                                                                      ttGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTg                                                                >  1:246070/1‑113 (MQ=255)                                                                                                      ttGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTg                                                                <  2:246070/113‑1 (MQ=255)                                                                                                         gcGAACAGCTTAAACGAAATAATAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCAt                         <  2:441214/148‑1 (MQ=255)                                                                                                                 cTTAAACGCAATAGCAAAACGATCCGTGAACATAAATTTAATGTGCTGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGGATCCCTAtgg                 >  1:162538/1‑147 (MQ=255)                                                                                                                  ttAAACGCAATATCAAAACGATCCGTGAACATAACTTTAATGTGAAGTCGAAACTCGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGGCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTTCATCCCTATTaa                >  1:19694/1‑149 (MQ=255)                                                                                                                     aaCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATcc             <  1:461015/149‑1 (MQ=255)                                                                                                                       cGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCAGTACTAAACAAATAATGTGTTGCATCCATATTAATCCGc           >  1:12872/1‑149 (MQ=255)                                                                                                                         cAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCat         <  2:151038/149‑1 (MQ=255)                                                                                                                                aaaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg  >  2:164768/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                        CCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCG  >  NC_000913_3_vch_tpiA/4604606‑4604881

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