Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 3,560,455 +G intergenic (‑2/+1) glpR ← / ← glpR pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA3,560,4551.G98.0% 160.9 / NA 49intergenic (‑2/+1)glpR/glpRpseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (19/29);  minor base A (1/0);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGT  >  NC_000913_3_pae_tpiA/3560309‑3560596
                                                                                                                                                   |                                                                                                                                             
cGATATCGTGAAACAGCGGTGAGCGAGTGGGTATTTGCTCCGCCACTTTGCGGGCGTTGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgc                                                                                                                                              <  2:572812/149‑3 (MQ=255)
       gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCa                                                                                                                                       <  2:491730/149‑1 (MQ=255)
       gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCCCACCGCCATCATCGCGCAGGAACa                                                                                                                                       >  2:262527/1‑149 (MQ=255)
       gATAAACAGCGTCGACCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCa                                                                                                                                       <  1:269897/149‑1 (MQ=255)
         tAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                                                                                                                     <  1:482873/149‑1 (MQ=255)
              aGCGGCGGGCCATTGGGGATTTGCTCCGCCAGTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTg                                                                                                                                <  2:915801/149‑1 (MQ=255)
                cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCt                                                                                                                              >  1:507506/1‑149 (MQ=255)
                  tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTcc                                                                                                                            <  1:539317/149‑1 (MQ=255)
                    gAGCCATTGGGGCTTTGATCCGCCACTTTGCGGGCGATGCGCTCTTTGTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGc                                                                                                                          >  2:321764/1‑149 (MQ=255)
                    gAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGc                                                                                                                          <  1:851584/149‑1 (MQ=255)
                     aGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTCGCTCCGcc                                                                                                                         >  1:90112/1‑149 (MQ=255)
                             gggATTTGCTCCGCCACTTTGCGGGCTATGCGCTTTTTTTCTTCGGTCTGGTTGGCCTGTCGACCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                                 <  2:792164/149‑1 (MQ=255)
                             gggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                                 >  1:361780/1‑149 (MQ=255)
                              ggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCGGCTCATTg                                                                                                                >  1:277700/1‑149 (MQ=255)
                               gATTTGCGCCGCCACTTTGCGGGCGATGCGGTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGa                                                                                                               <  2:590569/149‑1 (MQ=255)
                                gtttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                              <  2:657985/148‑1 (MQ=255)
                                 tttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCGCATTGAgg                                                                                                             >  2:516534/1‑149 (MQ=255)
                                     atCCGCCGCTTTGCGGGCGATGCGCCCGTGTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCAGCACCGTCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                         <  2:477405/148‑1 (MQ=255)
                                           cACTTTGCGGGCGATGCGCGCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                                   <  2:38536/149‑1 (MQ=255)
                                             cTTTGCGGTCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGCTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                                 <  1:61225/149‑1 (MQ=255)
                                                tGCGGGCGATCCGCTCTTTTCCTTCGGTCTGGGTGGCCGTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc                                                                                              <  1:643476/149‑1 (MQ=255)
                                                    ggtgATGCGCTCTTTTTGTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCg                                                                                          <  2:51907/146‑1 (MQ=255)
                                                     gCGAGGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                         <  2:367512/149‑1 (MQ=255)
                                                         tGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGAGCGTGCCACGGCGTGGTAACCGAACTGGAAGGCAGCGCCGCACCGCCATCATGGCGCGGGATCAGGTGTTGCTCCGCCAGCTCCTCGGGGGCGCGCGGAATAATCCGCGGACta                                                                                     >  2:651272/1‑148 (MQ=255)
                                                          gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGa                                                                                    <  1:297822/149‑1 (MQ=255)
                                                               tttttcCTTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGGAGCGCCGCACCGCCATGATGGCGCAGGATCAGAATATGCTCCGTCAGCTCATTGAGGTCGCGGCGAATAGGCTGCGGGCTGACGgag                                                                               <  2:4549/148‑1 (MQ=255)
                                                                  tttCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCGGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                                            <  2:108213/149‑1 (MQ=255)
                                                                      ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                        >  1:723603/1‑149 (MQ=255)
                                                                          gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                    >  1:624041/1‑149 (MQ=255)
                                                                          gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGAAGCGATGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                    >  1:883755/1‑149 (MQ=255)
                                                                                   gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                           <  1:640785/149‑1 (MQ=255)
                                                                                    ccTTGCGATCGTCGCACGGCGTGGTAACCCAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                          <  1:434865/149‑1 (MQ=255)
                                                                                     cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                         >  1:610774/1‑149 (MQ=255)
                                                                                     cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCGGACGGAGAAATGCTCTCCCAGCTCTTCgg                                                         >  2:328412/1‑149 (MQ=255)
                                                                                          gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTGGCGGCGAATAGTCTGCGGGCTGACGGAGGAATGCTCTACCAGCTCTTCGGTACTg                                                    >  2:206022/1‑149 (MQ=255)
                                                                                           aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                   <  1:676679/149‑1 (MQ=255)
                                                                                                 ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                             <  1:91333/149‑1 (MQ=255)
                                                                                                         tgtTAACCGAACTGGAAGGCAGAGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                                     <  2:132701/149‑1 (MQ=255)
                                                                                                              aCCGAACTGGAAGGCAGCCCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                                <  1:164945/149‑1 (MQ=255)
                                                                                                                 gAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                             >  1:619120/1‑149 (MQ=255)
                                                                                                                     tGGAAGGCAGCGCCGCACCGCCATGATGGCACAGGATCAGGGTTTGGTCCGCCACCTCAGTGAGTTCGCGCCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGCTACTCACATAACCCTCCTGTTTAACCAGTTCg                         >  2:496921/1‑149 (MQ=255)
                                                                                                                          ggCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTGACCAGTTCGataat                    >  2:922195/1‑149 (MQ=255)
                                                                                                                           gCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                   <  1:639665/149‑1 (MQ=255)
                                                                                                                             aGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                         <  1:851534/61‑1 (MQ=255)
                                                                                                                             aGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                         >  2:851534/1‑61 (MQ=255)
                                                                                                                             aGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                 >  1:873659/1‑149 (MQ=255)
                                                                                                                                  cgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgt            <  1:329020/149‑1 (MQ=255)
                                                                                                                                  cgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgt            <  1:321764/149‑1 (MQ=255)
                                                                                                                                   gcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgtg           >  2:858868/1‑149 (MQ=255)
                                                                                                                                            tgatgGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTgtgt  <  2:877512/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGT  >  NC_000913_3_pae_tpiA/3560309‑3560596

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: