Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,383,225 G→A 20.0% P35L (CCG→CTG)  pdhR ← DNA‑binding transcriptional regulator LldR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,383,2250GA20.0% 34.5 / 4.8 20P35L (CCG→CTG) pdhRDNA‑binding transcriptional regulator LldR
Reads supporting (aligned to +/- strand):  ref base G (8/8);  new base A (2/2);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

TGAAGGTGGACCCCAACGACTCCGCCACAAAATTGCCCCCGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGACAATGTCGTCGGACAGGCGCCGTTGGCGTACCTGATCAAAAACCATCACATGCTCTCCACAAACCCCCGGCAAATTTCAGGGGG  >  NC_002947/5383082‑5383366
                                                                                                                                               |                                                                                                                                             
tGAAGGTGAAACCCAACGCCTCCGCAACAAAATTTCCCCCGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGcc                                                                                                                                         <  2:247188/150‑1 (MQ=255)
          ccccAACGACTCCGCCACAAAATTGCCCCCGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGa                                                                                                                               <  1:191981/150‑1 (MQ=255)
                 gACTCCGCCACAAAATTGCCCCCGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGc                                                                                                                        <  1:532214/150‑1 (MQ=255)
                         cacaAAATTGCCCCCGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCAGCAGCCGCTGCCCCGAGGTCAGc                                                                                                                        <  1:33175/142‑1 (MQ=255)
                         cacaAAATTGCCCCCGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCAGCAGCCGCTGCCCCGAGGTCAGc                                                                                                                        <  1:33221/142‑1 (MQ=255)
                         cacaAAATTGCCCCCGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCAGCAGCCGCTGCCCCGAGGTCAGc                                                                                                                        >  2:33221/1‑142 (MQ=255)
                         cacaAAATTGCCCCCGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCAGCAGCCGCTGCCCCGAGGTCAGc                                                                                                                        >  2:33175/1‑142 (MQ=255)
                           caAAATTGCCCCCGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCa                                                                                                              >  1:393127/1‑150 (MQ=255)
                                 tGCCCCCGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAatca                                                                                                        >  2:267137/1‑150 (MQ=255)
                                    ccccGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATcc                                                                                                     >  1:65413/1‑150 (MQ=255)
                                           cTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAg                                                                                              >  2:274128/1‑150 (MQ=255)
                                               cGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGa                                                                                          <  2:393127/150‑1 (MQ=255)
                                                    gACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGAc                                                                                     <  1:267137/150‑1 (MQ=255)
                                                    gACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGAc                                                                                     <  1:274128/150‑1 (MQ=255)
                                                          cagcCCCTTGGCCACCAGATTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGACAATgtc                                                                               >  2:226853/1‑150 (MQ=255)
                                                                       accaGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGACAATGTCGTCGGACAGgcgc                                                                  >  1:356953/1‑150 (MQ=255)
                                                                                                      ggcgggACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGACAATGTCGTCGGACAGGCGCCGTTGGCGTACCTGATCAAAAACCATCACAt                                   <  2:65413/150‑1 (MQ=255)
                                                                                                                acCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGACAATGTCGTCGGACAGGCGCCGTTGGCGTACCTGATCAAAAACCATCACATGCTCTCcaca                         >  1:377390/1‑150 (MQ=255)
                                                                                                                   aaaCTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTCCCAGAATCATCCCTTCCAGCCGATCGACAATGTCGTCGGACAGGCGCCGTTGGCGTACCTGATCAAAAACCATCACATGCTCTCCACAAAc                      <  2:377390/150‑1 (MQ=255)
                                                                                                                                       gcTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGACAATGTCGTCGGACAGGCGCCGTTGGCGTACCTGATCAAAAACCATCACATGCTCTCCACAAACCCCCGGCAAATTTCAggggg  >  2:595081/1‑150 (MQ=255)
                                                                                                                                              cggcAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGACAATGTCGTCGGACAGGCGCCGTTGGCGTACCTGATCAAAAACCATCACATGCTCTCCACAAACCCCCGGCaaa            <  1:391895/133‑1 (MQ=255)
                                                                                                                                              cggcAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGACAATGTCGTCGGACAGGCGCCGTTGGCGTACCTGATCAAAAACCATCACATGCTCTCCACAAACCCCCGGCaaa            >  2:391895/1‑133 (MQ=255)
                                                                                                                                               ggcAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGACAATGTCGTCGGACAGGCGCCGTTGGCGTACCTGATCAAAAACCATCACATGCTCTCCACAAACCCCCGGCAAATTTCa       >  1:382841/1‑137 (MQ=255)
                                                                                                                                               ggcAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGACAATGTCGTCGGACAGGCGCCGTTGGCGTACCTGATCAAAAACCATCACATGCTCTCCACAAACCCCCGGCAAATTTCa       <  2:382841/137‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
TGAAGGTGGACCCCAACGACTCCGCCACAAAATTGCCCCCGCCCTGACGGCTGACCAGCAGCCCCTTGGCCACCAGCTTCTGGATCGCCTCACGCAGTGACGGGCGGGACACACCAAACTGCTCGGCCAGGACGCGCTCGGCCGGCAGCCGCTGCCCCGAGGTCAGCGTGCCTTCCAGAATCATCCCTTCCAGCCGATCGACAATGTCGTCGGACAGGCGCCGTTGGCGTACCTGATCAAAAACCATCACATGCTCTCCACAAACCCCCGGCAAATTTCAGGGGG  >  NC_002947/5383082‑5383366

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