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.4% 154.1 / ‑6.7 55A186E (GCA→GAA) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (25/28);  total (26/29)
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

CGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTAC  >  NC_000913_3_vch_tpiA/4604603‑4604888
                                                                                                                                               |                                                                                                                                              
cgaccgacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCGGCTTAAACGCAATATCAAAACGATCCGTGAACATa                                                                                                                                           >  1:464986/1‑149 (MQ=255)
       cTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTaa                                                                                                                                    <  1:564204/149‑1 (MQ=255)
         gAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATGGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtg                                                                                                                                  >  1:40917/1‑149 (MQ=255)
          aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtgt                                                                                                                                 >  1:279759/1‑149 (MQ=255)
          aaCAATCAGCGCTACATTAACCAAAGCCGGACATTAAGCCCAACGGAGAGGGAAATATTGCGCTTTATGGCCCGTGGCTACTCAATGACACAACTTCCCGAGCAGCTTAAACGCCATACCAAACCCAGCCCGGCACATAACTTTacagt                                                                                                                                 >  2:541839/1‑145 (MQ=255)
           aCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtgtg                                                                                                                                <  1:492928/149‑1 (MQ=255)
            cAATCAGTGGTACATTAAACAAAGCCGGACGTTAAGCCCTACGGAGAGAGAAATATTGCGCTTTATGGCGCGTGGCTACGCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGAGCCGTGAACATAAATTTAATGTGa                                                                                                                               <  2:685270/149‑1 (MQ=255)
             aaTCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAAAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGATTAAACGCAATCTCAAAACGATCCGTGAACATAAATTTAATGTGAt                                                                                                                              >  1:118948/1‑149 (MQ=255)
              aTCAGTGGTACATTAACCAAAGCCGGCCGATTATGAATACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTtaatgtaatg                                                                                                                             >  1:824551/1‑149 (MQ=255)
               tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTTGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTCAACATAAATTTAATGTGATGt                                                                                                                            >  2:697129/1‑149 (MQ=255)
               tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGTCATCGCACAACGATCCGTGAACATAAATTTAATGTGATGt                                                                                                                            <  1:409943/149‑1 (MQ=255)
                    ggTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGAAACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGaaa                                                                                                                       <  1:369115/149‑1 (MQ=255)
                        cATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTgg                                                                                                                   >  1:296123/1‑149 (MQ=255)
                         aTTAACCAAAGCCGGATGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATGTCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAACATggg                                                                                                                  >  1:728029/1‑149 (MQ=255)
                          taaaCCAATGCCGTACGTTAAGCCCGACGGAGATAGAAATATTGCGCTTTATGTCGCGTGGTTATTCAATGACAAAAATTGTAAAGCCACTTAAACTCAATATCAAACCCAGCCGTTAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                                 <  2:304833/147‑1 (MQ=255)
                              ccAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCa                                                                                                             >  1:143575/1‑149 (MQ=255)
                                 tagtcGGGCGTTAAGCGCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATGAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                                                                                          <  2:449835/145‑1 (MQ=255)
                                  aGCGGGACGTTGAGCGCGACGGAGAGGGAAATATTGCTCTTTATGTGGGGTGGCTACTCAATGACACCAATTGCCGGGCAGCTTAAACGCAATATCAAAACGTTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTc                                                                                                         <  2:71440/149‑1 (MQ=255)
                                     cGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                                                                                      >  2:827331/1‑149 (MQ=255)
                                      ggACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGCGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGAc                                                                                                     >  2:267447/1‑149 (MQ=255)
                                          gTTAAGCCCGACGGAGAGAGAAATATTCCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAg                                                                                                 <  2:460122/149‑1 (MQ=255)
                                                 ccGAGGGGGGGAGAAAGAGTGCGTTTTATGGCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                          <  2:279759/149‑1 (MQ=255)
                                                 ccGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                          <  2:296123/149‑1 (MQ=255)
                                                 ccGACGGAGAGAGAAATATTGCGCTTTAGGTTGCTGGGCTACTGAATGTCACAACTTGCCGAGCAGCTTAGACGCAATATGAAGACGATGCGTGAACATAAGTTTTATGTGATGTCGAACCTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                          <  2:629751/149‑1 (MQ=255)
                                                   gACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACt                                                                                                                     >  2:54195/1‑120 (MQ=255)
                                                   gACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACt                                                                                                                     <  1:54195/120‑1 (MQ=255)
                                                   gACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTgg                                                                                        <  1:478417/149‑1 (MQ=255)
                                                    gcGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                       <  1:208346/148‑1 (MQ=255)
                                                            agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAg                                                                               >  2:51513/1‑149 (MQ=255)
                                                                  atTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTc                                                                         >  1:215573/1‑149 (MQ=255)
                                                                   ttGCGCTTTATGTCGCGTGTGTACTCAATGACAAAAATTGTCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGACACTGGGCGTCAGTTCTCACGCAGGGTTGTTGGAGGCCGCAGATATCCt                                                                        <  2:533416/149‑1 (MQ=255)
                                                                   ttGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCt                                                                        >  1:551851/1‑149 (MQ=255)
                                                                                        tACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGa                                                   >  2:436407/1‑149 (MQ=255)
                                                                                        tACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGa                                                   >  2:397480/1‑149 (MQ=255)
                                                                                         aCTCAATGACCAAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaa                                                  <  2:506633/149‑1 (MQ=255)
                                                                                          cTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGaaa                                                 <  1:288474/149‑1 (MQ=255)
                                                                                           tCAATGACAAAAATTACCGAGCAAGCTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAAc                                                <  1:674539/149‑1 (MQ=255)
                                                                                             aaTGACCAAAATAGCCGAGGAGCTTAAACGCAATATCAAAACGTTCCGTGAACATAAATTTATTGTGATGGCGAAACGGGGCGTTAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                              <  2:140087/149‑1 (MQ=255)
                                                                                             aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                              >  1:638055/1‑149 (MQ=255)
                                                                                             aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGAGATTCTGTTATGTATGCGGCATTGCGAAACaa                                              >  2:71016/1‑149 (MQ=255)
                                                                                                  cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCATCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAt                                         <  2:551851/149‑1 (MQ=255)
                                                                                                      aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCTTGAACATAAATTTCATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCATGGTTGTTGGAGGCCGCAGATATTCTGTCATGTATGCGGCATTTCGCACCACGAAATgtgt                                     >  2:748986/1‑149 (MQ=255)
                                                                                                       aTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                                    <  1:397480/149‑1 (MQ=255)
                                                                                                        ttGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTg                                   <  2:215573/149‑1 (MQ=255)
                                                                                                               gcagcGTAAACGCAAGATCACAACGATCCGTGCACATAAATTTAATGGGATGTCGAAACTGGTCGTCAGTTCTGACGCAGGGTTGTTGGTGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCt                            <  1:556607/149‑1 (MQ=255)
                                                                                                                cagcTTAAACGCAATATCAAAACGATCCGTGAACATAAATATAACGTGAAGTGGAAAATGGGCATCCGGTATATAACAGGCTCGCTAGTGGCCGCAGATATGATGTTAGGTATGCGGCATCGCGAAACAAGTAATGTCTTGAATACCtt                           >  2:165101/1‑148 (MQ=255)
                                                                                                                      aaaCGAAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCTGATATTCTGGTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATc                     <  2:464986/149‑1 (MQ=255)
                                                                                                                            aaTATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCatg               >  2:657989/1‑149 (MQ=255)
                                                                                                                                  aaaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCg         >  1:669413/1‑149 (MQ=255)
                                                                                                                                  aaaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTCGTTTAAGGCCGCCGAAATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGAATCCCTCTTAATCCGCATGATGCCg         >  2:13308/1‑149 (MQ=255)
                                                                                                                                   aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGCTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg        >  1:716519/1‑149 (MQ=255)
                                                                                                                                    aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCggg       <  1:144294/149‑1 (MQ=255)
                                                                                                                                      ggagccGTGAACTTAAATTTTATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGtt     <  1:375082/145‑1 (MQ=255)
                                                                                                                                       gATCCGTGAAAATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGttt    <  1:291745/149‑1 (MQ=255)
                                                                                                                                         tCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTAc  <  1:827331/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTTAC  >  NC_000913_3_vch_tpiA/4604603‑4604888

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