Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 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_bme_pgi3,560,4551.G97.1% 102.0 / ‑3.8 34intergenic (‑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 . (1/0);  new base G (19/14);  total (20/14)
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

GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACC  >  NC_000913_3_bme_pgi/3560307‑3560581
                                                                                                                                                      |                                                                                                                              
gCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACGTTTCGGGCGCTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGGT‑GGc                                                                                                                                 >  1:173167/1‑149 (MQ=255)
        cGATAAACAGACTGGAGCGAGTGGGGATTTGGTCCGCCACTTTGCGGGCGATGCGCTTTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATc                                                                                                                         <  2:548393/149‑1 (MQ=255)
                 gCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGc                                                                                                                <  1:33652/149‑1 (MQ=255)
                   gTCGAGCCATTGGGGATTTGCTCGTCCACTTTGCGCGAGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCACCGCACCGCCAGGAT‑GGCGCAGCACCAGGTTGTGGTc                                                                                                              >  2:190245/1‑149 (MQ=255)
                    tCGAGCCATTGGGGATTTGCTCCGCAACTTTGAGGGCGATGCGCTCTTTTTATTCGGAATGGGTGGCCTTGGGATAGTGAAACGGCGTGTTAAACGAACTGGAAGGCACCGCCGCAGCGCCATCATAGGG‑CAGACGCAGGTTTTgtcca                                                                                                             >  2:318016/1‑145 (MQ=255)
                            ttGGGGATTTGGGACTGCACTTTGCGGGCGAGGCGCTCTTTTTTTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                     <  2:167675/149‑1 (MQ=255)
                            ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                     >  1:107815/1‑149 (MQ=255)
                            ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGAGCTCTTTTTCTTCGGTCTGGGTGGGCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCGGCACCGCCATGAT‑GGCGCAGAATCAGGTTTTGCTTCGCCGGCTc                                                                                                     >  1:87862/1‑149 (MQ=255)
                                ggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCTGTCTGGGTGGCCTTGCGGTCGTGCCACGGCGTGTTAACCGACATGGAATGCAGCGCCGCACCGCCATGAT‑GGCGCAACAATAGGTTTTTCTCGGCCAGCTCAttc                                                                                                 >  2:107589/1‑148 (MQ=255)
                                   tttGCCCCGCCACTTTGCGGGCGATGCGCTCTCTTTCTTCGGCCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAgg                                                                                              <  1:195279/149‑1 (MQ=255)
                                               cTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCTCACCGCGATGAT‑GGCGCAGGTTCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                  <  1:294001/149‑1 (MQ=255)
                                                              gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCCGCACCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGCGGACGCGCCGACTAGTCTACGGGCTGGCg                                                                   >  2:286606/1‑149 (MQ=255)
                                                              gcTCTTTTTCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGCCGCGGCGAATAGTCTGCGGGCTGACg                                                                   >  2:584676/1‑149 (MQ=255)
                                                                 ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                <  1:149281/149‑1 (MQ=255)
                                                                        ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                         <  1:70755/149‑1 (MQ=255)
                                                                        ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                         >  1:218495/1‑149 (MQ=255)
                                                                        ttcGGTCTGGGTGGCCTTGCGAGCGTGCCACGGCGTGTTAACCGAAGTGGAAGGCAGCGCCTCACCGCCATGAG‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                         <  2:657760/149‑1 (MQ=255)
                                                                             tCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGTTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                    <  2:327889/149‑1 (MQ=255)
                                                                                             ctCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCACCACCGCCATGAT‑GACGCAGCACCAGGTTTCACTCCGCCAGCCCGTTGAGATCGCGGCGAATAGTGTGCGGCCGGACGGAGAAATGCTCTACCAGGCCTTCGGtcttgc                                    >  2:356195/2‑144 (MQ=255)
                                                                                             aTCGTGCAACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                    >  2:638378/1‑149 (MQ=255)
                                                                                              tCGTCCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGGCGGAGAAGTGCTGTACTGGCCCTTCGGTAGTGAc                                   >  2:556737/1‑149 (MQ=255)
                                                                                                gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATTAT‑AGCGCAGGATCCGGTAGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGTCTCACGGAGAAATGCCCTACCAGCTCTTCGGTACTGACAt                                 >  2:83528/1‑149 (MQ=255)
                                                                                                gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                 >  1:720571/1‑149 (MQ=255)
                                                                                                   ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCGCCGCCATGAT‑GGCGCAGGAACAGGTTTTGCTCCGCCAGCTCATTTAGGTGGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCGCTAACAGGTCTTCGAGAGTGACATAAc                              >  2:69825/1‑149 (MQ=255)
                                                                                                         cGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                        <  1:350146/149‑1 (MQ=255)
                                                                                                            gtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTCCTGACATAACCCTGCGGtt                     >  1:146460/1‑149 (MQ=255)
                                                                                                             ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAGCCTGCTGttt                    >  1:325001/1‑149 (MQ=255)
                                                                                                             ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                    >  2:164865/1‑149 (MQ=255)
                                                                                                              tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTa                   <  2:107815/149‑1 (MQ=255)
                                                                                                               aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                  <  1:674272/149‑1 (MQ=255)
                                                                                                                aCCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                 >  1:386838/1‑149 (MQ=255)
                                                                                                                 ccGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCGGCCAGATCATTCAGGTCGCGTCGAATAGTCGGCGGGCTGACGGAGAAATGCTCTAGCAGCTCTTCGTTACTGACATAACCCTGCTAGTTAAcc                >  1:152329/1‑149 (MQ=255)
                                                                                                                  cGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCa               <  2:136267/149‑1 (MQ=255)
                                                                                                                           aGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa      <  2:720571/149‑1 (MQ=255)
                                                                                                                               aGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATCATAcc  >  2:132858/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                              
GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACC  >  NC_000913_3_bme_pgi/3560307‑3560581

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