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

ACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_vch_tpiA/4604605‑4604891                                                                                                                                              |                                                                                                                                                  accgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATaaa                                                                                                                                            >  1:424266/1‑149 (MQ=255)    gaatGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATGAAAACGATCCGTGAACATAAAttt                                                                                                                                         >  2:502846/4‑149 (MQ=255)        gAACAATCAGTGGTACATTAACCAAAGCCGGAAGTTAAGCCACACGGAGAGCACAATATCAGGCTTTATGTCGCGTAATAACTCAATGACGTAAATTGCCGCGCAGCTTAAACGCAGTAACAAAACGATCCGGCAACAAAAATCTAAtg                                                                                                                                     >  2:185358/1‑149 (MQ=255)         aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtgt                                                                                                                                    >  2:259419/1‑149 (MQ=255)            aaTCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGAt                                                                                                                                 >  2:246514/1‑149 (MQ=255)             aTCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGGGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAGCAAAACGATCCGTGAACATAAATTTAATGTGATg                                                                                                                                >  1:204730/1‑149 (MQ=255)              tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGtgttgt                                                                                                                               >  1:39340/1‑149 (MQ=255)              tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGt                                                                                                                               >  2:46630/1‑149 (MQ=255)                  tGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGGGAGAGAAATATTGCGCTTTATGTAGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGaa                                                                                                                           >  2:23548/1‑149 (MQ=255)                  tGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGaa                                                                                                                           <  1:585183/149‑1 (MQ=255)                      agaTTAACCAAGGCGGGACGTTCAAGCCGACTGAAATAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTATACACAATATCACAACGAGCCGTGAACATATATGTAATGTGATGTCGAAACTg                                                                                                                       <  2:369423/147‑1 (MQ=255)                       cATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTgg                                                                                                                      <  2:585153/149‑1 (MQ=255)                        aTTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCTCGTTATGTCGCGTGGCGACTCACTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTggg                                                                                                                     >  1:185537/1‑149 (MQ=255)                             ccAAAGCTTGACGTTAAGCCCTACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCa                                                                                                                <  2:650606/149‑1 (MQ=255)                              cAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAg                                                                                                               >  1:513383/1‑149 (MQ=255)                                     ggAAGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACACAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGAc                                                                                                        <  2:131870/149‑1 (MQ=255)                                         gTTAGGCCCGACGGAGAGAGAAATGGAACGCGTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAACAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAGGTAATGTCGAAAATGAGCGTCAGTTCTGACGCAg                                                                                                    >  2:599039/1‑149 (MQ=255)                                          ttAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAgg                                                                                                   <  2:513383/149‑1 (MQ=255)                                            aaGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgt                                                                                                 <  1:46630/149‑1 (MQ=255)                                                  gACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTgg                                                                                           >  2:475146/1‑149 (MQ=255)                                                   aCGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                          <  1:313853/149‑1 (MQ=255)                                                     ggAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                                                                                        <  2:704669/149‑1 (MQ=255)                                                      gagagagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGGCAGTTCTGACGCGGGGTTGTTGGAGGc                                                                                       >  1:176422/1‑149 (MQ=255)                                                       agagagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                                                                                            <  2:728769/143‑1 (MQ=255)                                                       agagagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                                                                                            >  1:728769/1‑143 (MQ=255)                                                         agagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGc                                                                                    >  2:435224/1‑149 (MQ=255)                                                              aaTATTGCGCTTTATTTCGCGGGGCTATTCAACGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACCTAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGata                                                                               <  1:281754/149‑1 (MQ=255)                                                                tatTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATAtt                                                                             >  2:117979/1‑149 (MQ=255)                                                                      tctctaagTCGCGTGTCTCCTCTAGGACACAAATCGCCGAGCACCTTAAACGAAATATCAAAACGATCCGTGAACATAAATTTAATTTGATGGCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTta                                                                       <  2:713772/142‑1 (MQ=255)                                                                      gcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTta                                                                       <  1:435224/149‑1 (MQ=255)                                                                       ctgtATGGCGCGTGGCTACGCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat                                                                      <  1:48344/146‑1 (MQ=255)                                                                                gcgTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACTCAGGGTTATTGGAGGCCGCAGTTTTTCTGTTATGTATGCGGc                                                             >  2:65142/1‑149 (MQ=255)                                                                                            aaTTACACAAATAGCCAAGCAGCTTAAACGCAAACTCAAAACGATCCGCGAACAAAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGCGGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATGTCGAACCaa                                                 >  2:584065/1‑149 (MQ=255)                                                                                            aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACag                                                 >  2:137428/1‑148 (MQ=255)                                                                                            aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                                 <  1:557651/149‑1 (MQ=255)                                                                                            aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTATGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaaaggg                                                 >  2:588874/1‑146 (MQ=255)                                                                                             aTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAg                                                <  1:246514/149‑1 (MQ=255)                                                                                                    aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGCGATTGTAGGGCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgtg                                         <  1:350524/149‑1 (MQ=255)                                                                                                       ttGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTg                                      <  1:117979/149‑1 (MQ=255)                                                                                                            gAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTg                                                                          <  2:224063/108‑1 (MQ=255)                                                                                                            gAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTg                                                                          >  1:224063/1‑108 (MQ=255)                                                                                                                agcTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTAt                             <  1:259419/149‑1 (MQ=255)                                                                                                                 gcTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAAGTGGGCGGCGATTGTGACGCAGGGTTGTTGGAGGCCGCTGATATTCTGTTATGTGTGCGGCATTGCGCAACAAGTAACGTTTTGGATACatatg                            >  2:438528/1‑148 (MQ=255)                                                                                                                  cTTAAACGCAAGATCAAAACGATCCGTGAACATAAATTTAATGTGCTGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGAGATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTa                           >  1:476905/1‑149 (MQ=255)                                                                                                                      aaCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTATGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTTTTGCATCCCGATTAATcc                       >  2:319309/1‑149 (MQ=255)                                                                                                                          cAATATCAAAACGTTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCat                   <  1:588874/149‑1 (MQ=255)                                                                                                                                  aaaCGATCCGAGAACATAAATTTAATGCGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg           <  1:96758/149‑1 (MQ=255)                                                                                                                                   aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCggg          <  1:85054/149‑1 (MQ=255)                                                                                                                                        tCCGTGAACATAAATCTAATATGATGTCGAAACTGGGCGCCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCCGTCATCTATACGGCATTGCGAAACAAGTAATGTGTTGCAACCCTATTAACCCGCATGATGCCGGGTTTAc     <  2:121335/149‑1 (MQ=255)                                                                                                                                           gTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGTTATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTc  >  1:597617/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                  ACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_vch_tpiA/4604605‑4604891

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