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.G95.1% 108.4 / ‑1.3 41intergenic (‑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 . (2/0);  new base G (22/17);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGT  >  NC_000913_3_pae_tpiA/3560307‑3560596
                                                                                                                                                       |                                                                                                                                             
gCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGc                                                                                                                                                <  2:572012/149‑1 (MQ=255)
         gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAATTTGTTGGCGATGCGCTGTTTTTCATCGTGCTTTTTTGCTTTGCGATCGTGCCACGGCGTGTTAACGGAACTGGAAGGCAGCGCCGGACCGCCATCA‑T‑CGCGCCGGATCa                                                                                                                                       >  2:460515/1‑149 (MQ=255)
          aTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCCCCGCCCTGA‑T‑GGCGCCGGATCAg                                                                                                                                      >  2:133490/1‑149 (MQ=255)
             aaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGtt                                                                                                                                   >  1:397721/1‑149 (MQ=255)
                  cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCt                                                                                                                              >  1:249836/1‑149 (MQ=255)
                       aGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTTTTTTTTCTTTGTGGCCTTGCGAGCGTGCCACGGCGAATTAACCGAACGGGAAGGCAGCCCCGCCCCGCCCGAACC‑GCC‑CAACATCAGGTTTTGCTCcccc                                                                                                                         >  2:108729/1‑149 (MQ=255)
                               gggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTACGACGGCTTGTTAACCGAACTGGACGGCAGCGCCGCACCGCCATGA‑C‑GGCGCAGGATCAGGTTTTGCTCCGCCCGCTCAtt                                                                                                                 >  2:325113/1‑149 (MQ=255)
                                        tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCACGG‑TAGGC‑CAGTGACGGGTATTTCTCCGCCAGCTCATTGAGGTcgcg                                                                                                        >  1:602417/1‑149 (MQ=255)
                                        tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATAGTGCCACGGCGTGTTAATCGAACTGGAAGGCAGAGACGCACCGCCATGA‑T‑GGCGCAGGATGAGGTATTGCTCCGCCGGTTGATTGAGGtatgg                                                                                                        >  2:33028/1‑145 (MQ=255)
                                          cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                                      >  1:5567/1‑149 (MQ=255)
                                                     gggCGATGCGCTCTTTTTTTTTGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGTCATGA‑T‑GGCGCAGGATCAGGTTTTCCTCCGCCAGCTCATTGAGGCTGCGGCGAAGAGTCAGc                                                                                           >  2:371644/1‑149 (MQ=255)
                                                             cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                                   >  1:520519/1‑149 (MQ=255)
                                                                        ttcGGGCGGGGTGGCCGTGAGATGGTGCCACGGCGTGTTAACCGAACTGGTAGGTTGCTCCGCTCCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                        <  2:426951/149‑1 (MQ=255)
                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                    <  2:397721/149‑1 (MQ=255)
                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGACTGCGTGCTGACGGAGAAATGCTCTAc                                                                    >  2:134103/1‑149 (MQ=255)
                                                                             tcaagggggCCTGGCGATCGTTCCACGGCGTGTAAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑C‑GGCGCAGAAACAGGTTTTGCACCGCCAGCTCATAGAGGACGCGGCGAATAGTCTGCGGGCTGGCGAAGAAAGCCTCTGcc                                                                   >  2:598944/8‑149 (MQ=255)
                                                                             tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                   >  1:225101/1‑149 (MQ=255)
                                                                                      ccTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                          <  2:520519/149‑1 (MQ=255)
                                                                                       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                         >  1:595130/1‑149 (MQ=255)
                                                                                       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                         <  2:225101/149‑1 (MQ=255)
                                                                                       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                         >  2:642209/1‑149 (MQ=255)
                                                                                            gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCTCCATGA‑T‑GGCGCAGGATCAGATTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGATATGCTCTACCAGCTCTTCGGTACTg                                                    >  1:581467/1‑149 (MQ=255)
                                                                                             aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                   >  2:247969/1‑149 (MQ=255)
                                                                                              tCGTGCCACGGCGTGTTAACCGAACTGGAACACAACCCCGCACCGCCATGA‑T‑AGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTATCCGGTACtacc                                                  >  2:463700/1‑146 (MQ=255)
                                                                                                   ccACGGCCTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                             <  1:108813/149‑1 (MQ=255)
                                                                                                    cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTTACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                            <  1:255547/149‑1 (MQ=255)
                                                                                                    cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                            <  2:432093/149‑1 (MQ=255)
                                                                                                        gggtgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgc                                        <  2:314562/147‑1 (MQ=255)
                                                                                                          gtggTAACCGAACGGGAAGGCAGCGCCGCTCCGCCATGA‑T‑GGCGCAGGATGAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                      <  2:2467/149‑1 (MQ=255)
                                                                                                           tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                                     >  1:538278/1‑149 (MQ=255)
                                                                                                                aCCGAACTGGAAGGCAGCGGCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                                <  2:112553/149‑1 (MQ=255)
                                                                                                                      ccggAAGGGAGCGTCCATCCGCCATGA‑T‑GGCGCAGGTGCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAGAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                          <  1:203623/147‑1 (MQ=255)
                                                                                                                      cTGGAAGGCAGCGCTGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAACAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                          <  1:208310/149‑1 (MQ=255)
                                                                                                                       tGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAGCCAGTTCg                         >  1:203569/1‑149 (MQ=255)
                                                                                                                             gCAGCGCCGCACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCGCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                   <  2:572220/149‑1 (MQ=255)
                                                                                                                                 cgccgcACCGCCATGA‑T‑GGCGCAGTATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCACAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACCCTGCTGTTTAACCAGTTCGGTACTACCGt               >  1:219216/1‑149 (MQ=255)
                                                                                                                                  gccgcACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt              <  1:642209/149‑1 (MQ=255)
                                                                                                                                     gcACCGCCATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgtg           >  1:267711/1‑149 (MQ=255)
                                                                                                                                       gtcGCCCTGA‑G‑GGCGCAGGATCGGGGTTTGCTCCGCCAGTTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCGGGCGACACGCTTTTGGAGCTCTTCGGTACTGACATAACCCTGCTGTTTATCCTGTTCGATAATACCGTTGTGAc         <  2:522055/147‑1 (MQ=255)
                                                                                                                                           ccATGA‑T‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg     <  1:239819/149‑1 (MQ=255)
                                                                                                                                             atga‑t‑gGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAccgttgtgccgttgtg   >  2:95219/1‑149 (MQ=255)
                                                                                                                                              tga‑t‑gGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTgtgt  <  2:293354/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGA‑T‑GGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGT  >  NC_000913_3_pae_tpiA/3560307‑3560596

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