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_tpiA4,604,6320CA97.4% 120.7 / ‑4.1 39coding (557/678 nt)bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (18/20);  minor base . (1/0);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGT  >  NC_000913_3_pae_tpiA/4604493‑4604770
                                                                                                                                           |                                                                                                                                           
cgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGAAGGATAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAttt                                                                                                                                   >  1:322129/1‑150 (MQ=255)
  acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGATCATAAATTTaa                                                                                                                                 <  1:375788/150‑1 (MQ=255)
    tGAACTATAAGTGGTACATTAACCAAAGGCGGACGTTAAGCCCGACGGAGAGTGAAATATTGCGCTTTATGTCGGGTGGCTACTGAATGACACAAATTGCCGAGCAGGTTAAACGCAGTATCAAAACGATCCGTGAACATAAATTTAAtg                                                                                                                               <  2:562858/150‑1 (MQ=255)
    tGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAGTTAAtg                                                                                                                               >  1:50180/1‑150 (MQ=255)
       aCAATCAGTTGTACATTAACCAAAGCCGGAGGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACAAAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGa                                                                                                                            <  2:531544/150‑1 (MQ=255)
       aCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtgtgg                                                                                                                            >  1:120681/1‑149 (MQ=255)
       aCAAGCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCTCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGa                                                                                                                            <  2:529538/150‑1 (MQ=255)
               tGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTAGGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGaaa                                                                                                                   >  1:560617/1‑150 (MQ=255)
                     aTTAACCAAAGCGGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGAATGGCAATTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                             <  2:365640/150‑1 (MQ=255)
                     aTTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                             <  2:140506/150‑1 (MQ=255)
                         aCCAAAGCGGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGATACTGGGCGTCa                                                                                                         <  1:376842/150‑1 (MQ=255)
                          ccAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATCTCAACACGATCCGTGAACATAAAGTTGAGTTGATGTCGAAACTCGGCGTGAg                                                                                                        >  2:73060/1‑150 (MQ=255)
                                  ggACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCTAAACTGGGCGTACGTTCTGACg                                                                                                >  2:338867/1‑150 (MQ=255)
                                         aaGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCCCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgtt                                                                                         <  2:424454/150‑1 (MQ=255)
                                                    agagagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATATCGAAACTGGACGTCAGTTCTGACGCAGGGTTGTTGGAGGccc                                                                              >  2:526729/1‑149 (MQ=255)
                                                    agagagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCACAGCGGTTTTAGCGCAATATCAAAACGATCGGTGAGCTCAAGGGTATTTCAAGTGCTAAACTTGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCg                                                                              >  2:425620/1‑150 (MQ=255)
                                                      agagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGGGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCa                                                                            >  1:581317/1‑150 (MQ=255)
                                                        agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgt                                                                                          >  2:2775/1‑134 (MQ=255)
                                                        agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgt                                                                                          <  1:2775/134‑1 (MQ=255)
                                                                 gcgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAGTGTGATGTCGAAACTGGGCGTCAGTCCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGtt                                                                 <  1:345613/150‑1 (MQ=255)
                                                                             gcgTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCa                                                     <  2:323626/150‑1 (MQ=255)
                                                                                    tACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaa                                              >  1:55455/1‑150 (MQ=255)
                                                                                     aCTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaaa                                             <  1:425620/150‑1 (MQ=255)
                                                                                       tCAATGACACAAATTGCCGAGCAGCTTAAACGCAAAAGCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACa                                           >  2:484130/1‑150 (MQ=255)
                                                                                       tCAATGACACAAATTGCCGAGCAGCTTAAACACAATATCAAAACGATCCGTGAACAAAAATTTAATGTGATGTCCAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGGCATCCTGTTAAGTATGCAGAATTGCGAAACa                                           >  2:240039/1‑150 (MQ=255)
                                                                                       tCAATGACAAAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACa                                           <  2:230522/150‑1 (MQ=255)
                                                                                            gACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTaa                                      <  2:261062/150‑1 (MQ=255)
                                                                                               acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgt                                   >  2:526126/1‑150 (MQ=255)
                                                                                               aacaaTTGCTGCGCAGGTTAAACGAAATATCAAAACGATCTGTGAACATAAATTTAATGTGACGTCGAACCTGGGCGTCAGTTCTGACGAAGGGTTGTTGGAGGCCGCAGATATTCTTTTATGTACGCGGCATTGCGAAACAAGTAAtgt                                   <  2:176857/147‑1 (MQ=255)
                                                                                                  aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                                >  2:331300/1‑150 (MQ=255)
                                                                                                  aaTTGCCCAGCAGCTTAAACGCAATATCAAAGCGTTCCGGGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                                <  1:498919/150‑1 (MQ=255)
                                                                                                      gCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCAt                            <  1:425244/150‑1 (MQ=255)
                                                                                                           gcatcTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTa                       <  1:250902/150‑1 (MQ=255)
                                                                                                            catcTTAAACGCAATATTAAAACGATCCTTGAAGATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCGGACGCAGGGTTGGTGGATGCCGCAGATATTCGGTTATGTATGCGGCATTGCGAAACAAGTAATTTGGTGCTTCCCTAt                      <  2:373255/147‑1 (MQ=255)
                                                                                                                ttAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAAt                  >  1:110033/1‑150 (MQ=255)
                                                                                                                           atCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGc       >  1:224286/1‑150 (MQ=255)
                                                                                                                              aacaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACACGTAATGTATTGCATACCTCTTAATCCGCATGATGCCgg    >  2:557082/4‑150 (MQ=255)
                                                                                                                               aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCATCAGTTCTGACGCAGGGTTGTTGTAGGGCGCAGATATTCTGTTATCTATGCCGCATTGCGAAACAAGTACTGTGTTGCATCCCTATTAATCCGCATGTTGCCggg   >  2:219952/1‑150 (MQ=255)
                                                                                                                                aaCGAGCGGTGAACATAAATTTAATGGGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGt  <  2:89086/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
CGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGT  >  NC_000913_3_pae_tpiA/4604493‑4604770

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: