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_tpiA4,604,7460CA96.2% 69.3 / NA 27A186E (GCA→GAA) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (15/10);  total (17/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGG  >  NC_000913_3_vch_tpiA/4604603‑4604882
                                                                                                                                               |                                                                                                                                        
cgaccgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGGGAGAGAAATATTGCGCTTTATGTCGCGTGGCTATTCAATGACACACAGTGCCGAGCAGCTTAAATGCAATATCAAAACGATCCGTGaaccaa                                                                                                                                     >  1:161836/1‑146 (MQ=255)
      acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTa                                                                                                                               >  1:300250/1‑149 (MQ=255)
        tGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAAt                                                                                                                             >  1:309676/1‑149 (MQ=255)
        tGAACAATCAGTGGTACATACACCAACGCCACACGTAAACCCCGCAGGTCAAAGAATTATAACGAATCAATTCGCGTGGCTACTCAATGTCACAAATTGCCGAGCAGCTTAACCTCAAAACCAAAACAATCCGTGGACATAAATTTAAt                                                                                                                             >  2:86843/1‑149 (MQ=255)
          aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGATAGAAATATTGGTCATTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACACAATATCAAAAAGATCCGTGAACAAAAATTTAAtgt                                                                                                                           >  1:324342/1‑149 (MQ=255)
                       aCATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAAAATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTg                                                                                                              >  1:138554/1‑149 (MQ=255)
                       aCATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACACAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTg                                                                                                              >  1:232208/1‑149 (MQ=255)
                               cAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGAATCGAGAACATAAATTTAATGAGATACCGAACCTGCGCGTCAg                                                                                                      >  2:71340/1‑149 (MQ=255)
                                            tAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGAAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                                                                                         <  1:83322/149‑1 (MQ=255)
                                                cccGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgt                                                                                     >  1:141519/1‑149 (MQ=255)
                                                cccGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgt                                                                                     >  1:400928/1‑149 (MQ=255)
                                                    aCGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGATACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                 <  1:53346/149‑1 (MQ=255)
                                                        agagagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGcc                                                                             >  1:17468/1‑149 (MQ=255)
                                                           gagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAGCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCa                                                                          <  2:103090/149‑1 (MQ=255)
                                                                       gcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTta                                                              <  2:17468/149‑1 (MQ=255)
                                                                                    gagCTACTCAATGACACAAATTGCCGATCAGCTTAAACGCAATATAAAAACGATCCGTGAACATAAATTTAATTGAATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCAtt                                                 <  1:404258/147‑1 (MQ=255)
                                                                                      acTACTCAATAACACAAATTGGCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAGTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGc                                               <  1:75071/148‑1 (MQ=255)
                                                                                        tACTCAAGGACACAAATTGCCGAGCAGTTTAAACACAATATCAAAACGATCCGTGAACATAAATTTAGTGTGATGTCGAAACTGGGCGCCAGTTCTGACGCAGTGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGa                                             <  2:324342/149‑1 (MQ=255)
                                                                                                 tcacacaTTTCCGCGCAGTTTAAACTCAAGATCAAGACGATCCGTGAATATAAATTTAATGTGATGTCGAAACGGGGCGGCAGTGCTGTCGCAGGGTTGTTTGAGGCTGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTaa                                    <  2:232208/148‑1 (MQ=255)
                                                                                                       aTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                              <  2:300173/149‑1 (MQ=255)
                                                                                                         tGCCGAGCAGCTTAAAGGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGc                            >  1:383593/1‑149 (MQ=255)
                                                                                                             gAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGATGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAGGTAATGTGTTGCATcc                        >  1:272472/1‑149 (MQ=255)
                                                                                                                   cTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTACTGTGTTGCATCCCTATTa                  >  2:120737/1‑149 (MQ=255)
                                                                                                                    ttAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTaa                 >  2:198452/1‑149 (MQ=255)
                                                                                                                          gCAATATCAAATTGATTCGTGAACATAAACTCGATGTTATGACAAAACGGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCa           <  1:80664/149‑1 (MQ=255)
                                                                                                                            aaTATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTATGATATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTCATCCGAatg         >  1:109867/1‑149 (MQ=255)
                                                                                                                                tCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGc     >  2:253341/1‑149 (MQ=255)
                                                                                                                                   aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTAGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATGGCGAAACAAGGAATGTGTTGCATCCCTATTAATCCGCCTTATGCCgg  >  2:156657/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                        
CGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGG  >  NC_000913_3_vch_tpiA/4604603‑4604882

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: