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

GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGT  >  NC_000913_3_pae_tpiA/3560307‑3560596
                                                                                                                                                       |                                                                                                                                             
gCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGc                                                                                                                                                <  2:572012/149‑1 (MQ=255)
         gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCAATTTGTTGGCGATGCGCTGTTTTTCATCGTGCTTTTTTGCTTTGCGATCGTGCCACGGCGTGTTAACGGAACTGGAAGGCAGCGCCGGACCGCCATCATCGCGCCGGATCa                                                                                                                                       >  2:460515/1‑149 (MQ=255)
          aTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCCCCGCCCTGATGGCGCCGGATCAg                                                                                                                                      >  2:133490/1‑149 (MQ=255)
             aaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGtt                                                                                                                                   >  1:397721/1‑149 (MQ=255)
                  cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCt                                                                                                                              >  1:249836/1‑149 (MQ=255)
                       aGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTTTTTTTTCTTTGTGGCCTTGCGAGCGTGCCACGGCGAATTAACCGAACGGGAAGGCAGCCCCGCCCCGCCCGAACCGCCCAACATCAGGTTTTGCTCcccc                                                                                                                         >  2:108729/1‑149 (MQ=255)
                               gggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTACGACGGCTTGTTAACCGAACTGGACGGCAGCGCCGCACCGCCATGACGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCAtt                                                                                                                 >  2:325113/1‑149 (MQ=255)
                                        tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCACGGTAGGCCAGTGACGGGTATTTCTCCGCCAGCTCATTGAGGTcgcg                                                                                                        >  1:602417/1‑149 (MQ=255)
                                        tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATAGTGCCACGGCGTGTTAATCGAACTGGAAGGCAGAGACGCACCGCCATGATGGCGCAGGATGAGGTATTGCTCCGCCGGTTGATTGAGGtatgg                                                                                                        >  2:33028/1‑145 (MQ=255)
                                          cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                                      >  1:5567/1‑149 (MQ=255)
                                                     gggCGATGCGCTCTTTTTTTTTGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGTCATGATGGCGCAGGATCAGGTTTTCCTCCGCCAGCTCATTGAGGCTGCGGCGAAGAGTCAGc                                                                                           >  2:371644/1‑149 (MQ=255)
                                                             cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                                   >  1:520519/1‑149 (MQ=255)
                                                                        ttcGGGCGGGGTGGCCGTGAGATGGTGCCACGGCGTGTTAACCGAACTGGTAGGTTGCTCCGCTCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                        <  2:426951/149‑1 (MQ=255)
                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                    <  2:397721/149‑1 (MQ=255)
                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGACTGCGTGCTGACGGAGAAATGCTCTAc                                                                    >  2:134103/1‑149 (MQ=255)
                                                                             tcaagggggCCTGGCGATCGTTCCACGGCGTGTAAACCGAACTGGAAGGCAGCGCCGCACCGCCATGACGGCGCAGAAACAGGTTTTGCACCGCCAGCTCATAGAGGACGCGGCGAATAGTCTGCGGGCTGGCGAAGAAAGCCTCTGcc                                                                   >  2:598944/8‑149 (MQ=255)
                                                                             tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                   >  1:225101/1‑149 (MQ=255)
                                                                                      ccTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                          <  2:520519/149‑1 (MQ=255)
                                                                                       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                         >  1:595130/1‑149 (MQ=255)
                                                                                       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                         <  2:225101/149‑1 (MQ=255)
                                                                                       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                         >  2:642209/1‑149 (MQ=255)
                                                                                            gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCTCCATGATGGCGCAGGATCAGATTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGATATGCTCTACCAGCTCTTCGGTACTg                                                    >  1:581467/1‑149 (MQ=255)
                                                                                             aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                                   >  2:247969/1‑149 (MQ=255)
                                                                                              tCGTGCCACGGCGTGTTAACCGAACTGGAACACAACCCCGCACCGCCATGATAGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTATCCGGTACtacc                                                  >  2:463700/1‑146 (MQ=255)
                                                                                                   ccACGGCCTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                             <  1:108813/149‑1 (MQ=255)
                                                                                                    cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTTACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                            <  1:255547/149‑1 (MQ=255)
                                                                                                    cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                            <  2:432093/149‑1 (MQ=255)
                                                                                                        gggtgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgc                                        <  2:314562/147‑1 (MQ=255)
                                                                                                          gtggTAACCGAACGGGAAGGCAGCGCCGCTCCGCCATGATGGCGCAGGATGAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                      <  2:2467/149‑1 (MQ=255)
                                                                                                           tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                                     >  1:538278/1‑149 (MQ=255)
                                                                                                                aCCGAACTGGAAGGCAGCGGCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                                <  2:112553/149‑1 (MQ=255)
                                                                                                                      ccggAAGGGAGCGTCCATCCGCCATGATGGCGCAGGTGCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAGAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                          <  1:203623/147‑1 (MQ=255)
                                                                                                                      cTGGAAGGCAGCGCTGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAACAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                          <  1:208310/149‑1 (MQ=255)
                                                                                                                       tGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAGCCAGTTCg                         >  1:203569/1‑149 (MQ=255)
                                                                                                                             gCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCGCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                   <  2:572220/149‑1 (MQ=255)
                                                                                                                                 cgccgcACCGCCATGATGGCGCAGTATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCACAGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACCCTGCTGTTTAACCAGTTCGGTACTACCGt               >  1:219216/1‑149 (MQ=255)
                                                                                                                                  gccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt              <  1:642209/149‑1 (MQ=255)
                                                                                                                                     gcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgtg           >  1:267711/1‑149 (MQ=255)
                                                                                                                                       gtcGCCCTGAGGGCGCAGGATCGGGGTTTGCTCCGCCAGTTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCGGGCGACACGCTTTTGGAGCTCTTCGGTACTGACATAACCCTGCTGTTTATCCTGTTCGATAATACCGTTGTGAc         <  2:522055/147‑1 (MQ=255)
                                                                                                                                           ccATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg     <  1:239819/149‑1 (MQ=255)
                                                                                                                                             atgatgGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAccgttgtgccgttgtg   >  2:95219/1‑149 (MQ=255)
                                                                                                                                              tgatgGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTgtgt  <  2:293354/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGT  >  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: