Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,233,640 G→A 11.8% R17R (CGC→CGT)  glpR ← DNA‑binding transcriptional repressor GlpR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,233,6400GA11.8% 84.5 / 3.9 34R17R (CGC→CGT) glpRDNA‑binding transcriptional repressor GlpR
Reads supporting (aligned to +/- strand):  ref base G (17/13);  new base A (2/2);  total (19/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 = 7.54e-01

TCGATGCTCGAGTCGTAAGCTGCGCCGCCGTGGTAACGGCGCAGCAGGTTGAGCTCGGCAAGCTGATTGATGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAGGGGATGACGGCCAGTAGACGATTGGG  >  NC_002947/1233503‑1233767
                                                                                                                                         |                                                                                                                               
tcGATGCTCGAGTCGTAAGCTGCGCCGCCGTGGTAACGGCGCAGCAGGTTGAGCTCGGCAAGCTGATTGATGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTAcc                                                                                                                     >  1:689292/1‑150 (MQ=255)
         gAGTCGTAAGCTGCGCCGCCGTGGTAACGGCGCAGCAGGTTGAGCTCGGCAAGCTGATTGATGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAgg                                                                                                            >  1:334531/1‑150 (MQ=255)
                    tGCGCCGCCGTGGTAACGGCGCAGCAGGTTGAGCTCGGCAAGCTGATTGATGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATttgttgttg                                                                                                 >  1:650088/1‑150 (MQ=255)
                          gccgTGGTAACGGCGCAGCAGGTTGAGCTCGGCAAGCTGATTGATGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCgggg                                                                                           >  2:157478/1‑150 (MQ=255)
                                     ggCGTAGCAGGGGGTGCTCAGCAAGCGGATTGATGTGGCGTCGGATGGTCTGGGGATTGACAACGAACAGCTGCGCCATTTCCTGGATACGGACATGCCCGCGTTCGCGTACCAGTGCGAGGGATTGTTTTTGTCGGGGGGGCCGAtgta                                                                                >  2:365762/1‑147 (MQ=255)
                                         cagcagGTTGAGCTCGGCAAGCTGATTGATGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATggg                                                                            <  2:650088/150‑1 (MQ=255)
                                           gcagGTTGAGCTCGGCAAGCTGATTGATGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCg                                                                          <  1:133930/150‑1 (MQ=255)
                                             agGTTGAGCTCGGCAAGCTGATTGATGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCgga                                                                        >  1:798564/1‑149 (MQ=255)
                                                              cTGATTGATGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGa                                                       <  2:677533/150‑1 (MQ=255)
                                                               tgattgatGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGAc                                                      >  1:76648/1‑150 (MQ=255)
                                                                      tGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCt                                               <  2:774960/150‑1 (MQ=255)
                                                                      tGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCt                                               <  1:157478/150‑1 (MQ=255)
                                                                          gcgTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCggggg                                                                                          <  2:119043/103‑1 (MQ=255)
                                                                          gcgTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCggggg                                                                                          >  1:119043/1‑103 (MQ=255)
                                                                                gATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCgg                                                                         <  2:582735/114‑1 (MQ=255)
                                                                                gATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCgg                                                                         >  1:582735/1‑114 (MQ=255)
                                                                                         gggAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAg                            <  2:798564/150‑1 (MQ=255)
                                                                                            ggtcacacCGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCTGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAgggg                         <  2:656000/146‑1 (MQ=255)
                                                                                              tGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAGGGGAt                       >  1:1222/1‑150 (MQ=255)
                                                                                               gACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAGGGGATg                      >  2:737399/1‑150 (MQ=255)
                                                                                                  aaCGAACAGCTGCGCCATTTCTTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTCAGATGATGCCGTAGGAAGGGGATGACg                   >  1:72/1‑150 (MQ=255)
                                                                                                  aaCGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAGGGGATGACg                   <  2:1222/150‑1 (MQ=255)
                                                                                                  aaCGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCACGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATggg                                                                            >  1:230782/1‑93 (MQ=255)
                                                                                                  aaCGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCACGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATggg                                                                            <  2:230782/93‑1 (MQ=255)
                                                                                                  aaCGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCACGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATggg                                                                            <  2:230889/93‑1 (MQ=255)
                                                                                                  aaCGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCACGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATggg                                                                            >  1:230889/1‑93 (MQ=255)
                                                                                                      aaCAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTa                                 <  2:244674/132‑1 (MQ=255)
                                                                                                      aaCAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTa                                 >  1:244674/1‑132 (MQ=255)
                                                                                                            tGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAGGGGATGACGGCCAGTAGAc         >  1:43076/1‑150 (MQ=255)
                                                                                                                 cATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAGGGGATGACGGCCAGt             >  2:15978/1‑141 (MQ=255)
                                                                                                                 cATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAGGGGATGACGGCCAGt             <  1:15978/141‑1 (MQ=255)
                                                                                                                   tttCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGCGGCGTTAGGAAGGGGATGACGGCCAGTAGACGATTggg  >  1:243331/1‑150 (MQ=255)
                                                                                                                               tGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAGGGGATGACGGCCAGTAGACGATTgg   <  1:240903/137‑1 (MQ=255)
                                                                                                                               tGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAGGGGATGACGGCCAGTAGACGATTgg   >  2:240903/1‑137 (MQ=255)
                                                                                                                                         |                                                                                                                               
TCGATGCTCGAGTCGTAAGCTGCGCCGCCGTGGTAACGGCGCAGCAGGTTGAGCTCGGCAAGCTGATTGATGTCGCGTCGGATGGTCTGGGGAGTGACAACGAACAGCTGCGCCATTTCCTCGATACTGACATAGCCGCGTTCGCGTACCAGTTCGAGGATTTGTTGTTGGCGGGGGGGCAGATTCATGGGCGGTCCTTTGGGGCTGCCGGACAAATTCTGCAATGATGCCGTAGGAAGGGGATGACGGCCAGTAGACGATTGGG  >  NC_002947/1233503‑1233767

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