Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 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_0009133,560,4551.G96.5% 86.0 / ‑3.1 29intergenic (‑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/1);  new base G (14/14);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTG  >  NC_000913/3560316‑3560587
                                                                                                                                             |                                                                                                                                    
gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCCGATGGcgca                                                                                                                                    >  1:69393/1‑140 (MQ=255)
gATAAACAGCGTCGAGCCATTGCGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCATCGCCATGATGGcgca                                                                                                                                    <  2:602905/143‑4 (MQ=255)
   aaaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGa                                                                                                                                 >  2:643741/1‑143 (MQ=255)
        tcctCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                                                                                                            <  1:438951/140‑1 (MQ=255)
            cGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGACCAGGcttt                                                                                                                        >  2:535172/1‑143 (MQ=255)
                 cATTGGGGATGTGCGCCGGCACTTTGCGGGCGATGCGCTCTGTTGCTCCGGTCTGGGTGGCCTTGCGATCGTGCCGCGGCGTGTTAACCGAACTGGAAGGCAGCGCCACACCGCCATGATGGCGCAGGATCAGGTTTTGCTcc                                                                                                                   <  1:201559/143‑1 (MQ=255)
                  aTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGAGTCAGGTGTCGCGCCg                                                                                                                  >  2:417748/1‑143 (MQ=255)
                   ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGc                                                                                                                 >  1:152511/1‑143 (MQ=255)
                   ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGCGTGGCCTTGCGAGCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGCCGCAGGAGCAGGTTGTGCTCCGc                                                                                                                 >  2:455108/1‑143 (MQ=255)
                          tttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCCGCGCCGCACCGCCCTGATGGCGCACAATCAGGTTCGGCTCCGCCGGcccc                                                                                                          >  2:424063/1‑140 (MQ=255)
                             gCTCCGCCTCTTTGCGGGCGATGCGCTCTTTTGCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                       <  2:140274/143‑1 (MQ=255)
                             gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGCGCAGGATCAGGGTTTGCTCCGCCAGCTCATTg                                                                                                       >  1:439077/1‑143 (MQ=255)
                                ccGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAgg                                                                                                    >  2:635250/1‑143 (MQ=255)
                                          gCGGGGGATGCGCTCGTTTTCTTCGGTCTGGGTGGCGTTGCGAGCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                          <  2:171200/143‑1 (MQ=255)
                                          gCGGGCGATGCGCTCTTTTGCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                          <  2:178123/143‑1 (MQ=255)
                                                  tGCTCTCTTTTTCTTCGGTCTGGGTGCCCTCGCGATCGTGCCACGGCGTGCTAACCGAACTGGAAGGCAGCGCCGCACCGCCAGGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                  <  1:642921/143‑1 (MQ=255)
                                                       tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                             >  2:35732/1‑143 (MQ=255)
                                                           tttCTGCGGTCGGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCGCCGCCATGATGGCGCAGGCTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                         <  1:657723/143‑1 (MQ=255)
                                                              cttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCACCGCACCGCCATGACGGCGCAGGATCAGGTTTTGCTCCGCCCGCCCAGTGAGGGCGCGGCGACTAGTCGGCGGGCTGACGgag                                                                      >  2:440302/1‑143 (MQ=255)
                                                               ttggggcTGGGTGGCCTTGCGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTCTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgaga                                                                     <  2:468146/137‑1 (MQ=255)
                                                                      tGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGGGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                              <  1:643741/143‑1 (MQ=255)
                                                                        ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                                                            <  1:246128/143‑1 (MQ=255)
                                                                              cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAGCAGGTTTTGCTCCGCCAGCTCATCGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                      >  2:410477/1‑143 (MQ=255)
                                                                                   gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                 <  2:439077/143‑1 (MQ=255)
                                                                                          ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGTTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                          >  1:196248/1‑143 (MQ=255)
                                                                                                    tgagccGGACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                <  1:550972/139‑1 (MQ=255)
                                                                                                         cGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATGACCctgctggt                           >  1:578021/1‑141 (MQ=255)
                                                                                                               ggAAGGCGGTGCCGCACCGCCATGATGGGGCAGGGTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCa                     <  2:167314/143‑1 (MQ=255)
                                                                                                                gAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTGGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCCGGCTGACGGAGAAACGCTCTACCAGCCCTTCGCTACAGCCATGACCCTGCTGTTCAAcccg                    >  2:97313/1‑141 (MQ=255)
                                                                                                                    gCAGCGCCGCACCTCCATGATGGCGCAGGATCAGGTTTTGCTCTGCCATCTTATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCg                <  1:125136/143‑1 (MQ=255)
                                                                                                                                  ccATGAGGGCGCGGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCGGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgtg  <  1:242689/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                    
GATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTG  >  NC_000913/3560316‑3560587

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: