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.2% 177.7 / NA 57A186E (GCA→GAA) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (26/29);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTACTTC  >  NC_000913_3_pae_tpiA/4604491‑4604777
                                                                                                                                             |                                                                                                                                                 
agcgacCGAAAAATCAGGGGTACATAAACCAAAGCCGTACGTGAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATaaa                                                                                                                                            <  2:635933/147‑1 (MQ=255)
 ccgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGCTCCGTGAACATAAAt                                                                                                                                           <  2:8111/149‑1 (MQ=255)
 ccgacTGAACAATAAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAt                                                                                                                                           <  2:694158/149‑1 (MQ=255)
   gacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTAGTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAttt                                                                                                                                         >  1:729551/1‑149 (MQ=255)
   gacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAttt                                                                                                                                         >  1:775072/1‑149 (MQ=255)
       gAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtg                                                                                                                                     >  2:333021/1‑149 (MQ=255)
        aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACACCAAAACGGTGATTGACGAGCCGCTTAAACGGAATATCAAAACGATCGGTGCAAAAAAATATAAtgt                                                                                                                                    >  1:396769/1‑149 (MQ=255)
           aaTCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGAAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGCGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATATGAt                                                                                                                                 >  2:774760/1‑149 (MQ=255)
             tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGGGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATAAAAACGATCCGTGAACATAAATTTAATGTGATGt                                                                                                                               >  1:173052/1‑149 (MQ=255)
                        ttAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTTGGc                                                                                                                    >  2:476783/1‑149 (MQ=255)
                            ccAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCCATATCAAAACGATCCGGGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCa                                                                                                                >  2:514164/1‑149 (MQ=255)
                             cAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTAGGCGCCAg                                                                                                               >  2:473720/1‑149 (MQ=255)
                                 gCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATCAATGTCATCTGATTTCGAAACTGCGCGCCAGTTCt                                                                                                           >  2:475671/1‑149 (MQ=255)
                                   cGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGGACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                                                                                         >  1:92673/1‑149 (MQ=255)
                                       cGTTAAGCCCGACGGAGAGAGAAATATTGCGTTTTATGTTGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCa                                                                                                     <  2:170077/149‑1 (MQ=255)
                                                 gACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTgg                                                                                           <  1:141660/149‑1 (MQ=255)
                                                   cGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAACCCAACCGTGAACATAAATTTAATGTGATGTCGACACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAg                                                                                         <  1:23986/149‑1 (MQ=255)
                                                   cGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAAGTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAg                                                                                         <  1:363926/149‑1 (MQ=255)
                                                     gagagagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACACAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGc                                                                                       <  1:233355/149‑1 (MQ=255)
                                                                  tGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTg                                                                          <  1:807471/149‑1 (MQ=255)
                                                                      cTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat                                                                      <  1:744732/149‑1 (MQ=255)
                                                                             tCGCGTGGCTACTCAATGACACAAATTGCCGCGAAGCTTATACGTATTATCAAAACGTTCCGTGAACATAAATTTAATGTGATGTCGCAACTGGGCGTCAGTTCTGACGGAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCg                                                               <  2:276824/149‑1 (MQ=255)
                                                                               gcgTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGTTATTCTGTTATGTATGCGGc                                                             >  1:83773/1‑149 (MQ=255)
                                                                               gcgTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGc                                                             >  1:788614/1‑149 (MQ=255)
                                                                                cgTGGCTACTCAATGACACAAATTGACGAGCAGCTTAAACGCAATATCAAAACCATCCGTGAACCTAAATTTAATGTGATGTTGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCTCAGATATTCTGTTATGTATGCGGCa                                                            <  1:740124/149‑1 (MQ=255)
                                                                                 gTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCAt                                                           >  1:363546/1‑149 (MQ=255)
                                                                                  tGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCAtt                                                          <  2:523432/149‑1 (MQ=255)
                                                                                    gCTACTCAATGACAAAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGc                                                        <  2:601308/149‑1 (MQ=255)
                                                                                     cTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTggg                                                                                                                     <  2:463717/87‑1 (MQ=255)
                                                                                     cTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTggg                                                                                                                     >  1:463717/1‑87 (MQ=255)
                                                                                      tACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGCTGTCGAAACAAGGCGTCAGTTATGACGCCGGGTTGTTGGAGGCCGCCGATATTCTGTACTGTATGCGGCATTACGa                                                      >  1:374828/1‑149 (MQ=255)
                                                                                        cacAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTAAGTCCTTTTGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaaa                                                    <  2:342073/147‑1 (MQ=255)
                                                                                           aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGt                                                                                                                               <  1:796695/71‑1 (MQ=255)
                                                                                           aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGt                                                                                                                               >  2:796695/1‑71 (MQ=255)
                                                                                           aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                                    <  1:489740/82‑1 (MQ=255)
                                                                                           aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                                    >  2:489740/1‑82 (MQ=255)
                                                                                           aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGACGGGTTGTTGGAGACCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                                 >  2:482090/1‑149 (MQ=255)
                                                                                              gACACAAATTGCCGAGCAGCTTAAACGAAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTa                                              <  1:573834/149‑1 (MQ=255)
                                                                                                cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAt                                            <  1:333021/149‑1 (MQ=255)
                                                                                                 acaAATTGCCGAGCAGCTTAAACGCAATATGAAAACGATCCGGGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATGCTGCTATGTATGCGGCAGTGCGAAACAAGTaagg                                           >  1:210136/1‑147 (MQ=255)
                                                                                                 acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                           >  1:108081/1‑149 (MQ=255)
                                                                                                 acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                                           >  1:356156/1‑149 (MQ=255)
                                                                                                   aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAgaa                                                                                                             <  1:282755/81‑3 (MQ=255)
                                                                                                   aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAgaa                                                                                                             >  2:282755/1‑79 (MQ=255)
                                                                                                                     aaCGCAATAGCAAAACGATCCGGGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGGTGGATGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATcc                       <  2:497348/149‑1 (MQ=255)
                                                                                                                             atCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatgatg               <  2:301954/149‑1 (MQ=255)
                                                                                                                                aaaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg            >  1:356465/1‑149 (MQ=255)
                                                                                                                                aaaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCAGTTATGTATGCGGCATTGCGAAACACGTAATGTGTTGCATCCCTATGAATCCGCATGATGCCg            >  2:755732/1‑149 (MQ=255)
                                                                                                                                 aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg           <  2:775072/149‑1 (MQ=255)
                                                                                                                                  aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATAtt                                                                             <  1:713311/82‑1 (MQ=255)
                                                                                                                                  aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATAtt                                                                             >  2:713311/1‑82 (MQ=255)
                                                                                                                                  aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTGATCCGCATGATGCCggg          <  1:792000/149‑1 (MQ=255)
                                                                                                                                  aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCggg          <  2:635902/149‑1 (MQ=255)
                                                                                                                                  aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCggg          <  2:795182/149‑1 (MQ=255)
                                                                                                                                  aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCggg          <  2:173052/149‑1 (MQ=255)
                                                                                                                                      aTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTa      >  1:347036/1‑149 (MQ=255)
                                                                                                                                          gTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGACCGCAGATATTATGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTCATCCGCATGATGCCGGGTTTACTTc  >  1:803928/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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑