Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,604,045 (A)7→8 coding (569/726 nt) yjjQ → DNA‑binding transcriptional repressor YjjQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,0381.A95.7% 72.8 / ‑3.9 24Q188K (CAA→AAA) yjjQDNA‑binding transcriptional repressor YjjQ
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base A (13/9);  minor base G (0/1);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

GCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCG  >  NC_000913/4603962‑4604109
                                                                              |                                                                       
gCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaa                                                                      >  1:931943/1‑80 (MQ=255)
   gCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaaa                                                                  <  1:931897/81‑1 (MQ=255)
      aTTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCaaaaaaaa                                                                 >  2:931933/1‑78 (MQ=255)
         aGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTg                                                             >  2:931930/1‑81 (MQ=255)
          ggTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGc                                                            <  2:931898/81‑1 (MQ=255)
              aaGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATg                                                        <  2:931899/81‑1 (MQ=255)
                      ggCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACt                                                <  2:931900/81‑1 (MQ=255)
                            aaCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAAt                                          <  1:931901/81‑1 (MQ=255)
                            aaCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAAt                                          <  2:931902/81‑1 (MQ=255)
                             aCCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCTCAAAAAAAACAGTGCGATGCGTCGACTGGCAAtt                                         >  2:931952/1‑81 (MQ=255)
                              ccAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTc                                        <  2:931903/81‑1 (MQ=255)
                               cAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTca                                       <  2:931904/81‑1 (MQ=255)
                               cAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTca                                       >  2:931931/1‑81 (MQ=255)
                                aGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCGAAAAAAACAGTGCGATGCGTCGACTGGCAATTcac                                      <  2:931905/81‑1 (MQ=255)
                                     gCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCaa                                 >  1:931932/1‑81 (MQ=255)
                                      cTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAAc                                <  2:931906/81‑1 (MQ=255)
                                            tGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGa                           >  2:931941/1‑80 (MQ=255)
                                                 aaGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTa                     >  2:931937/1‑81 (MQ=255)
                                                  aGAAAAGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTa                     <  1:931907/80‑1 (MQ=255)
                                                    aaaaGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATg                   >  1:931951/1‑80 (MQ=255)
                                                      aaGTAATAAAACTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCAt                >  2:931940/1‑81 (MQ=255)
                                                             aaaaCTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaaa         >  1:931945/1‑81 (MQ=255)
                                                              aaaCTATCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATaaa         >  2:931948/1‑80 (MQ=255)
                                                                   aTCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgc   >  1:931950/1‑81 (MQ=255)
                                                                    tCAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgcg  >  2:931957/1‑81 (MQ=255)
                                                                     cAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgcg  >  1:931946/1‑80 (MQ=255)
                                                                     cAGCGCCCAAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAgcgcg  >  1:931947/1‑80 (MQ=255)
                                                                              |                                                                       
GCGGCTATTAGGTAAAGGATGGGGCATCAACCAGATAGCTTCATTGCTTAAGAAAAGTAATAAAACTATCAGCGCCCAAAAAAACAGTGCGATGCGTCGACTGGCAATTCACAGCAACGCTGAAATGTATGCATGGATAAATAGCGCG  >  NC_000913/4603962‑4604109

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