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

GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACC  >  NC_000913_3_bme_pgi/3560307‑3560581
                                                                                                                                                      |                                                                                                                              
gCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACGTTTCGGGCGCTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGGTGGc                                                                                                                                 >  1:173167/1‑149 (MQ=255)
        cGATAAACAGACTGGAGCGAGTGGGGATTTGGTCCGCCACTTTGCGGGCGATGCGCTTTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                         <  2:548393/149‑1 (MQ=255)
                 gCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGc                                                                                                                <  1:33652/149‑1 (MQ=255)
                   gTCGAGCCATTGGGGATTTGCTCGTCCACTTTGCGCGAGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCACCGCACCGCCAGGATGGCGCAGCACCAGGTTGTGGTc                                                                                                              >  2:190245/1‑149 (MQ=255)
                    tCGAGCCATTGGGGATTTGCTCCGCAACTTTGAGGGCGATGCGCTCTTTTTATTCGGAATGGGTGGCCTTGGGATAGTGAAACGGCGTGTTAAACGAACTGGAAGGCACCGCCGCAGCGCCATCATAGGGCAGACGCAGGTTTTgtcca                                                                                                             >  2:318016/1‑145 (MQ=255)
                            ttGGGGATTTGGGACTGCACTTTGCGGGCGAGGCGCTCTTTTTTTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                     <  2:167675/149‑1 (MQ=255)
                            ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                     >  1:107815/1‑149 (MQ=255)
                            ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGAGCTCTTTTTCTTCGGTCTGGGTGGGCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCGGCACCGCCATGATGGCGCAGAATCAGGTTTTGCTTCGCCGGCTc                                                                                                     >  1:87862/1‑149 (MQ=255)
                                ggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCTGTCTGGGTGGCCTTGCGGTCGTGCCACGGCGTGTTAACCGACATGGAATGCAGCGCCGCACCGCCATGATGGCGCAACAATAGGTTTTTCTCGGCCAGCTCAttc                                                                                                 >  2:107589/1‑148 (MQ=255)
                                   tttGCCCCGCCACTTTGCGGGCGATGCGCTCTCTTTCTTCGGCCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAgg                                                                                              <  1:195279/149‑1 (MQ=255)
                                               cTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCTCACCGCGATGATGGCGCAGGTTCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                  <  1:294001/149‑1 (MQ=255)
                                                              gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCCGCACCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGCGGACGCGCCGACTAGTCTACGGGCTGGCg                                                                   >  2:286606/1‑149 (MQ=255)
                                                              gcTCTTTTTCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGCCGCGGCGAATAGTCTGCGGGCTGACg                                                                   >  2:584676/1‑149 (MQ=255)
                                                                 ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                <  1:149281/149‑1 (MQ=255)
                                                                        ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                         <  1:70755/149‑1 (MQ=255)
                                                                        ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                         >  1:218495/1‑149 (MQ=255)
                                                                        ttcGGTCTGGGTGGCCTTGCGAGCGTGCCACGGCGTGTTAACCGAAGTGGAAGGCAGCGCCTCACCGCCATGAGGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                         <  2:657760/149‑1 (MQ=255)
                                                                             tCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGTTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                    <  2:327889/149‑1 (MQ=255)
                                                                                             ctCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCACCACCGCCATGATGACGCAGCACCAGGTTTCACTCCGCCAGCCCGTTGAGATCGCGGCGAATAGTGTGCGGCCGGACGGAGAAATGCTCTACCAGGCCTTCGGtcttgc                                    >  2:356195/2‑144 (MQ=255)
                                                                                             aTCGTGCAACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                    >  2:638378/1‑149 (MQ=255)
                                                                                              tCGTCCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGGCGGAGAAGTGCTGTACTGGCCCTTCGGTAGTGAc                                   >  2:556737/1‑149 (MQ=255)
                                                                                                gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATTATAGCGCAGGATCCGGTAGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGTCTCACGGAGAAATGCCCTACCAGCTCTTCGGTACTGACAt                                 >  2:83528/1‑149 (MQ=255)
                                                                                                gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                 >  1:720571/1‑149 (MQ=255)
                                                                                                   ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCGCCGCCATGATGGCGCAGGAACAGGTTTTGCTCCGCCAGCTCATTTAGGTGGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCGCTAACAGGTCTTCGAGAGTGACATAAc                              >  2:69825/1‑149 (MQ=255)
                                                                                                         cGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                        <  1:350146/149‑1 (MQ=255)
                                                                                                            gtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTCCTGACATAACCCTGCGGtt                     >  1:146460/1‑149 (MQ=255)
                                                                                                             ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAGCCTGCTGttt                    >  1:325001/1‑149 (MQ=255)
                                                                                                             ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                    >  2:164865/1‑149 (MQ=255)
                                                                                                              tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTa                   <  2:107815/149‑1 (MQ=255)
                                                                                                               aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                  <  1:674272/149‑1 (MQ=255)
                                                                                                                aCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                 >  1:386838/1‑149 (MQ=255)
                                                                                                                 ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCGGCCAGATCATTCAGGTCGCGTCGAATAGTCGGCGGGCTGACGGAGAAATGCTCTAGCAGCTCTTCGTTACTGACATAACCCTGCTAGTTAAcc                >  1:152329/1‑149 (MQ=255)
                                                                                                                  cGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCa               <  2:136267/149‑1 (MQ=255)
                                                                                                                           aGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa      <  2:720571/149‑1 (MQ=255)
                                                                                                                               aGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATCATAcc  >  2:132858/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                              
GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACC  >  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: