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,6320CA98.1% 163.7 / NA 54A186E (GCA→GAA) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (30/22);  total (32/22)
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

GACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGT  >  NC_000913_3_pae_tpiA/4604486‑4604770
                                                                                                                                                  |                                                                                                                                          
gaccgaccgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGaac                                                                                                                                          >  1:285358/1‑149 (MQ=255)
    gaccgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATAGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCAAGCAGCTTAAACGCAATATCAACACGCTCCGTGCACAtca                                                                                                                                      >  2:595267/1‑147 (MQ=255)
     accgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATaaa                                                                                                                                     <  1:624681/149‑1 (MQ=255)
         acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCCATATCAAAACGATCCGTGAACCTAAATTTa                                                                                                                                 >  2:861465/1‑149 (MQ=255)
         acTGAACAATCAGTGGTACATTAACCAAAGCAGATCCCCGAGGCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTa                                                                                                                                 >  1:497720/1‑149 (MQ=255)
            gAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCAACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGACGAGCCGCTTAAACGCAATATCAAAAAGACTCGAGAACATAAATTTAAtg                                                                                                                              >  1:183160/1‑149 (MQ=255)
                gcactgTGGGAGATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGAt                                                                                                                          <  2:462087/144‑1 (MQ=255)
                aaTCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGCGAGAAATATTGCGCTTTATGTCGCGTGGGTACTCAAGGACACAAATTGCCGAGCAGCCTAAACGCAATATCAAAACGATCCGGGAACATAAATTTAAGGTGAt                                                                                                                          >  1:603290/1‑149 (MQ=255)
                          aCATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGACCCGTGAACATAAATTTAATGTGATGTCGAAACTg                                                                                                                <  1:419774/149‑1 (MQ=255)
                             ttAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                             >  1:654714/1‑149 (MQ=255)
                                 ccAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTggg                                                                                                              >  2:312091/1‑144 (MQ=255)
                                 ccAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTggg                                                                                                              <  1:312091/144‑1 (MQ=255)
                                    aaGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                                                                                      >  1:712969/1‑149 (MQ=255)
                                     aGCCGAACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTc                                                                                                     <  1:155414/149‑1 (MQ=255)
                                         ggACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGAc                                                                                                 >  1:132935/1‑149 (MQ=255)
                                                    ccGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                      <  1:316752/149‑1 (MQ=255)
                                                         ggAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                                                                                 <  1:613790/149‑1 (MQ=255)
                                                             agagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGCCACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGc                                                                             <  2:572442/149‑1 (MQ=255)
                                                                gaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGCAACTGGGCGTCAGTCCTGACGCAGGGATGTTGGAGGCCGCAGa                                                                          <  2:497720/149‑1 (MQ=255)
                                                                gaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGa                                                                          <  1:297593/149‑1 (MQ=255)
                                                                         cgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATCTTCTGtt                                                                 >  2:232754/1‑149 (MQ=255)
                                                                         cgcTTTATGTCCGGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGtt                                                                 <  1:431585/149‑1 (MQ=255)
                                                                           cTTTATGTCGGGTGGCTGCTCAATGACACAAATTGCCGAGGAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat                                                               <  2:105595/149‑1 (MQ=255)
                                                                           cTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGCCACTTCTGACGCAGGGGTGTTGGAGGCCGCAGATATTCTGTtat                                                               >  2:425037/1‑149 (MQ=255)
                                                                            tttATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTGCTGACGCAAGGGTGTTGGAAGCCGCAGATATTCTGTtatg                                                              >  2:598065/1‑149 (MQ=255)
                                                                            tttATGTCGCGTGGCTACTCAATGACACAAATTACCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATATGATGTCGAAACTGAGCATCAGTTATGACACAGGGTTGTTGGAGACCGAAGATATTCTGTtatg                                                              >  2:805627/1‑149 (MQ=255)
                                                                            tttATGTCGCGTGGCTACTCAATGAAACAAATTGGCGAGCAGCTTAAACGAAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatg                                                              <  2:235396/149‑1 (MQ=255)
                                                                                  tCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCg                                                        >  1:400229/1‑149 (MQ=255)
                                                                                   cgcgTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTCTGCgg                                                       >  2:102364/1‑149 (MQ=255)
                                                                                           tACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGa                                               <  2:297258/149‑1 (MQ=255)
                                                                                             cTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATGTTCTGTTATGTATGCGGCATTGCGaaa                                             >  1:433283/1‑149 (MQ=255)
                                                                                              tCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATGTAATGTGACGTTGAAACTGGGCGTCAGTTCCGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAAc                                            >  2:164127/1‑149 (MQ=255)
                                                                                                aaTGACACAAATTGCCGAGGAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                          <  2:69810/149‑1 (MQ=255)
                                                                                                aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                          >  1:236026/1‑149 (MQ=255)
                                                                                                      acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                    >  2:630076/1‑149 (MQ=255)
                                                                                                      acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                    >  1:15881/1‑149 (MQ=255)
                                                                                                      acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCAGTGTACAGAAATTTAATATGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                    >  1:747543/1‑149 (MQ=255)
                                                                                                        aaaTTGCCGAGCAGCTTAAACGTAATATCAAAACGATCAGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCTCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgtg                                  <  2:794535/149‑1 (MQ=255)
                                                                                                         aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATgtgt                                 <  1:245431/149‑1 (MQ=255)
                                                                                                          aTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                                <  2:806469/149‑1 (MQ=255)
                                                                                                                gAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAAATGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTGAAATATTTATGTGGCATGGCGAAACAAGTAATGTGTTGCATcc                          >  1:278492/1‑149 (MQ=255)
                                                                                                                       ttAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTaa                   >  1:675345/1‑149 (MQ=255)
                                                                                                                       ttAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTATGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTaa                   >  1:717284/1‑149 (MQ=255)
                                                                                                                           aCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCg               >  2:941013/1‑149 (MQ=255)
                                                                                                                           aCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCg               >  2:156538/1‑149 (MQ=255)
                                                                                                                             gCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTGTTAATCCGCa             >  1:724924/1‑149 (MQ=255)
                                                                                                                             gCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCCCa             >  2:609703/1‑149 (MQ=255)
                                                                                                                             gCAATATCAAAACGAGCCGTGAACATAAATTTAATGTGATGTCGAAACTGGCCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATGCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGct             >  2:258010/1‑148 (MQ=255)
                                                                                                                             gCAAGACCAAAACGATCCGTGAACATAAAGTTAATGTGATGTCGAAACTGGGCGTCAGTCCTAACGCAGGGTTGTTGGAGGCCGCAGATATACTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCa             <  1:834424/149‑1 (MQ=255)
                                                                                                                               aaTATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatg           >  2:46563/1‑149 (MQ=255)
                                                                                                                                    cAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGcc      <  2:712969/149‑1 (MQ=255)
                                                                                                                                      aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCCTCTCTATTCATCCGCATGATGGCgg    >  2:440432/1‑149 (MQ=255)
                                                                                                                                      aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg    <  1:148827/149‑1 (MQ=255)
                                                                                                                                       aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCgggg   >  2:159646/1‑149 (MQ=255)
                                                                                                                                        aCGATGCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCGGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGt  <  2:243627/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
GACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGT  >  NC_000913_3_pae_tpiA/4604486‑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: ‑