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

ACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_pae_tpiA/4604491‑4604777                                                                                                                                              |                                                                                                                                                  accgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGATAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATaaa                                                                                                                                            <  2:984954/149‑1 (MQ=255)     acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGCGACAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTa                                                                                                                                        >  1:463332/1‑149 (MQ=255)       tGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGCGCAGATCAATCGCAATATCAAAACGATCCGTGAACATAAATTTAAt                                                                                                                                      <  1:773711/149‑1 (MQ=255)        gaagaaTCAGAGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATAACACAAATTGCCGAGCAGATTAAACGCAATATCAAAACGATCCGTGAACATAAATTAAAtg                                                                                                                                     >  2:182620/1‑149 (MQ=255)         aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtgt                                                                                                                                    <  1:708004/149‑1 (MQ=255)              tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAAAGATCCGTGAACATAAATTTAATGTGATGt                                                                                                                               >  1:125407/1‑149 (MQ=255)              acAGTGGTAAATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGt                                                                                                                               <  1:755124/148‑1 (MQ=255)                aGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAACCGATCCGTGAACATAAATCTAATGTGATGTCg                                                                                                                             >  2:741710/1‑149 (MQ=255)                                       aCGTAAAGCCCGACGGAGAAAGAAATATTGCGCTTTATGTCGCTTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGc                                                                                                      <  1:92489/149‑1 (MQ=255)                                         gTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGAACTGACGccc                                                                                                    <  1:752680/149‑3 (MQ=255)                                          ttAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATTAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTTACGCAgg                                                                                                   >  1:887646/1‑149 (MQ=255)                                           tAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                                                                                                  <  2:500037/149‑1 (MQ=255)                                                cagACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGTTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                             <  1:17238/147‑1 (MQ=255)                                                   aCGGAGAAAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                          <  1:614204/149‑1 (MQ=255)                                                       agagagaAATATTCCGCTTGTTGTAGCTTGGCTACTCAATAACACAAATTTCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTGGGAGGcc                                                                                      <  2:926431/149‑1 (MQ=255)                                                          gagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGACGTCAGTTCTGAAGCAGGGTTGTTGGAGGCCGCa                                                                                   >  1:66710/1‑149 (MQ=255)                                                           agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAACTGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAg                                                                                  >  1:982712/1‑149 (MQ=255)                                                           agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCCATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAg                                                                                  >  1:412996/1‑149 (MQ=255)                                                              aaTATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGata                                                                               <  1:526028/149‑1 (MQ=255)                                                                   tGCGCTTTTTTTTGGGTGGGTAATCAATGACACAAATTGCCGAGCAGCGTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTg                                                                          <  1:152925/149‑1 (MQ=255)                                                                    gcgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGt                                                                         <  1:953202/149‑1 (MQ=255)                                                                       ccttATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                                                                                         <  1:177215/112‑1 (MQ=255)                                                                       ccttATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                                                                                         >  2:177215/3‑114 (MQ=255)                                                                        tttATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATATCGTAACTGGTCGTCAGTTCTACCGCAGGGTTGTTGGGGGCCGCAGAACTTCTGTtgtg                                                                     >  2:367147/1‑149 (MQ=255)                                                                            tGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTgg                                                                                                                      <  2:283636/96‑1 (MQ=255)                                                                            tGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTgg                                                                                                                      >  1:283636/1‑96 (MQ=255)                                                                                    ggCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGAGCGTCAGTTCTGACGCAAGGTTGTTGGAGGCCGCAGATATTATGTTATGTATGCGTCAttt                                                         >  1:83824/1‑148 (MQ=255)                                                                                    ggCTACTCAATGACACAAATTGACGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTg                                                         <  2:566781/149‑1 (MQ=255)                                                                                            aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATATATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                                 <  1:949321/149‑1 (MQ=255)                                                                                                acacAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGCAACAAGTaa                                             >  2:875743/1‑149 (MQ=255)                                                                                                acacAAAGTGGCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATTTGATGTCGAAACCGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTaa                                             <  1:492042/149‑1 (MQ=255)                                                                                                  acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGAGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                           >  1:897458/1‑149 (MQ=255)                                                                                                    aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgtg                                         >  2:103129/1‑149 (MQ=255)                                                                                                    aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGACGCAGAAATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgtg                                         >  1:81551/1‑149 (MQ=255)                                                                                                      aTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                                       <  1:979357/149‑1 (MQ=255)                                                                                                          ccGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAAATGGGAGTCAATGCTTACTCAGTTTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGGAACAAGTAATGTGTTGCAt                                   >  2:899262/1‑149 (MQ=255)                                                                                                           cGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGTGCGTCAGTTCTGAAGCAGGGTTATTGGAAGCCGAAGATATTCTGTTATGTATGAGGCATTTCGAAAAAAGTAATGTGTTGTAtg                                  >  1:1264/1‑148 (MQ=255)                                                                                                            gAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGACATGTCGAAACAAGTAATGTGTTGCATcc                                 <  2:749816/149‑1 (MQ=255)                                                                                                                       aCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCg                      <  1:28282/149‑1 (MQ=255)                                                                                                                         gCAATATCAATATGTTCCGTGAACATAAGTTTAATGTGATGTCGAAACTGGGCGTCAGTGCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCa                    <  1:182620/149‑1 (MQ=255)                                                                                                                         gCAATATCAAAACGATCCATGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTAACGCAGGGTTGTTGGAGGCCGCAGATATTTTGTTATGTATGCGGCATTGCGAAACAAATAATGTGTTGCATCCCTCTTAATCCGCa                    >  2:179187/1‑149 (MQ=255)                                                                                                                          cAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTTGCGAATTGTCCTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCat                   <  1:81086/149‑1 (MQ=255)                                                                                                                           aaTATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTGAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATCCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatg                  >  1:680701/1‑149 (MQ=255)                                                                                                                                 aaaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGGGGCATTGCGAAACAAGTAATGGGTTGCATCCCGATTAATCCGCATGATGccc            >  2:306949/1‑148 (MQ=255)                                                                                                                                  aaaCGATCCGTGAACATAAATTTAATGTGATGTCGACACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGTCGCAGATATTATGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCACTATTAATCCGCTTGATGCCgc           >  1:180844/1‑148 (MQ=255)                                                                                                                                  aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATTCGGCATTGCGAAACAAGTAATGTGTTGCATCCATATTAATCCGCATGATGCCgg           >  1:35076/1‑149 (MQ=255)                                                                                                                                   aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGTGGCATCGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCggg          >  1:815490/1‑149 (MQ=255)                                                                                                                                       aTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATACCGGGTTTa      >  1:78119/1‑149 (MQ=255)                                                                                                                                           gTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCATTTCTGACGCAGGATTGTTGGAGGCCGTAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTc  >  1:30943/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                  ACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_pae_tpiA/4604491‑4604777

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