Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 122,761 C→T H224Y (CAT→TAT)  pdhR → DNA‑binding transcriptional dual regulator PdhR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913122,7610CT95.2% 61.1 / ‑3.6 21H224Y (CAT→TAT) pdhRDNA‑binding transcriptional dual regulator PdhR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (13/7);  total (14/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

TGGACGCGGAATCAAACGCCGTACTCCAGTATCAGATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGAT  >  NC_000913/122519‑123000
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
tGGACGCGGAATCAAACGCCGTACTCCAGTATCAGATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGcc                                                                                                                                                                                                                                         >  1:203699/1‑251 (MQ=255)
                    gTACTCCAGTATCAGATTGCCGTCACCGAAGCGGCCCACACTGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCt                                                                                                                                                                                                                     >  1:242979/1‑251 (MQ=255)
                                 aGATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGt                                                                                                                                                                                                        <  1:79100/251‑1 (MQ=255)
                                             ccGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGCGcgccgc                                                                                                                                                                                            >  2:191855/1‑251 (MQ=255)
                                                               gggTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGcgt                                                                                                                                                                          <  2:52311/250‑1 (MQ=255)
                                                                       cTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTggg                                                                                                                                                                  >  1:178177/1‑250 (MQ=255)
                                                                          cATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGCAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGc                                                                                                                                                                >  1:140110/1‑250 (MQ=255)
                                                                           aTCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCa                                                                                                                                                               >  1:113965/1‑250 (MQ=255)
                                                                             cTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACg                                                                                                                                                            >  2:122033/1‑251 (MQ=255)
                                                                                                                        aGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACAGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGAtgt                                                                                                                 <  1:111493/251‑1 (MQ=255)
                                                                                                                                            tcgcggcgcgATATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCAtata                                                                                              <  2:190195/250‑1 (MQ=255)
                                                                                                                                                                            ccACCGAACCCGCATATTTGAAGCGATAATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAAAAGTCGTGAAGAGAGCCGCCGAGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTa                                                             <  2:171834/250‑1 (MQ=255)
                                                                                                                                                                                      cGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTc                                                   >  1:47402/1‑251 (MQ=255)
                                                                                                                                                                                                aaGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTCTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGt                                         <  2:203699/251‑1 (MQ=255)
                                                                                                                                                                                                                               gcgcgAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGt          >  2:69120/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGttatt      >  1:137707/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGttatta     <  1:175559/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aaGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGttatt      >  1:136123/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aaGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGttatt      >  1:130715/1‑250 (MQ=255)
                                                                                                                                                                                                                                       catcgcatcgCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGttatt      >  2:141505/1‑247 (MQ=255)
                                                                                                                                                                                                                                          cgcatcgcTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATtagat  >  1:34549/1‑248 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
TGGACGCGGAATCAAACGCCGTACTCCAGTATCAGATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGAT  >  NC_000913/122519‑123000

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: