Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,253,002 C→T intergenic (+17/+83) dhaR → / ← ycgV DNA‑binding transcriptional dual regulator DhaR/putative autotransporter adhesin YcgV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,253,0020CT95.7% 59.5 / ‑4.3 23intergenic (+17/+83)dhaR/ycgVDNA‑binding transcriptional dual regulator DhaR/putative autotransporter adhesin YcgV
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (6/16);  minor base . (1/0);  total (7/16)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTCTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATATCCTCCGTTGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGG  >  NC_000913/1252776‑1253245
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
gaccAGGCGACATATGATGTCAGCTCATACCGCCTTTCCACCAGTCTGGCATTTGCGGAAGTTGAAAAATAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCAGGGCAGTTTAAGCGCCGGGTATTAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTAtt                                                                                                                                                                                                                              <  2:153778/248‑1 (MQ=255)
                              cGCCTTTCCCCCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTTGTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCTGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATTCATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGc                                                                                                                                                                                               <  2:373477/251‑1 (MQ=255)
                                     ccaccaggcgttCATGTGCGGAAGTTGAAAAAGAGGTAATTAATAACGCAGGCCAGGTCACAGGCTGTCGCATGCAGGAAATGTCGGCTTTACTTGTGATCGGCCTAACTATGCTGTGGCTGAAAATGAAGTAACATGGCATTGATGCATGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATGTTTGATACATCAGAGGGATCTGTATTT‑CATTCCGGCGCACGCTAACAAGTTTCAGCATCGttt                                                                                                                                                                                        <  2:453721/240‑1 (MQ=255)
                                       accaGTCCGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCCGCTTTAATTGGGATCGGCCGCCCTACGCTGTGGCGGGAAATGAAGCAACAAGGCATTGATGCAGGGCAGGTTAAGAGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGAACTGTATTT‑CATTCCGGCGCACGCTAAAAATTTTCAGCATCGtttt                                                                                                                                                                                       >  1:331781/1‑249 (MQ=255)
                                                 cATTTGCGTAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACATGCGGTCGCATGCAGGAAATGTCGGCTTTACTTGTGATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATGGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTTTATTT‑CATTCCGGCGCACTCTATCAATTTTCAGCAGCGTTTACGGGCTTGTCt                                                                                                                                                                            <  2:84082/251‑1 (MQ=255)
                                                    ttGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCCTTCAGGAAATTTCGGCTTTACTTGGGGTCGGGCGCACTACGCCGTGGCGGGAAAAGAAGCAACCTGGCATTTTTGCAGGGCAGTTTAAACCCCCGGGATGAAAAACAGAAACGCATTT‑TTTTACATCAGAGAGATCTGTTTTTTCTTTCCG‑CGCCCGCTAAAAAATTTTAACCTCGTTTAAGGGCTTTTCtttt                                                                                                                                                                         >  2:309835/1‑248 (MQ=255)
                                                              tGAAAAGGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTa                                                                                                                                                                <  1:217199/250‑1 (MQ=255)
                                                                     gagGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGa                                                                                                                                                        <  1:429736/251‑1 (MQ=255)
                                                                                  aCGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATAtcctcc                                                                                                                                           <  1:111065/251‑1 (MQ=255)
                                                                                        cccAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATATCCTCCGTTGAc                                                                                                                                     <  1:350563/251‑1 (MQ=255)
                                                                                                   ggCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATATCCTCCGTTGACTTGCTTCTGAt                                                                                                                          <  1:299709/251‑1 (MQ=255)
                                                                                                   ggCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTGCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATATCCTCCGTTGACTTGCTTCTGAt                                                                                                                          <  2:504049/251‑1 (MQ=255)
                                                                                                     cGGTCGCATTCAGGAAATGTCTGCTTTACTTGGGATCGGCAGCACTCCGCTGTGTCGTAAAATGAAGCAACATGTCATTGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTAGATTCATCAGAGTTATCTGTATTT‑CATTCCGGTGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAATCGATATCCTCCGTTGACTTTCTTCTGATCa                                                                                                                        <  1:754/251‑1 (MQ=255)
                                                                                                                   aaaTGTCGGCTTTACTTGGGAGCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATATCCTCCGTTGACTTGCTTCTGATCAAAGAGGTTACCCTg                                                                                                          <  2:222929/251‑1 (MQ=255)
                                                                                                                           gctgtacgtGGGATCGTGCGCACTTCTCTGTGGCGGAAAATGAAGCATCCTTGCATTGATGCAGGGCAGGTTAAGCGCCGGGTATGAAAGACATAAACG‑ATTTTTGAGTCACCAGAGTGATCTGTATTT‑CATGCCCGCGTCCGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTAGCCCACACTTAAAAGCTGAAGCGTTATCCTCCTTTGACTTGGTTCTGATCAAAGAGTTTTCCCTGCGTGTAAt                                                                                                  <  2:296280/243‑1 (MQ=255)
                                                                                                                                    gggATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCTCACTTAAAAGCTGAAGCGATATCCTCCGTTGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCg                                                                                          <  1:36584/250‑1 (MQ=255)
                                                                                                                                                                         aaCATGGCATTGATGCAGGGCAGTTTAAGCGACGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATATCCTCCGTTGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGAc                                                    >  1:133223/1‑251 (MQ=255)
                                                                                                                                                                          acattgcatttctGCAGGGCAGGTTAAGCGCCTGGTATGAAAGACAGAATCG‑ATTTTTGATTCATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATATCCTCCGTGGACTTGCTGCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACt                                                   <  2:170252/239‑1 (MQ=255)
                                                                                                                                                                                     aTGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATAtcctcc                                                                                                                                           <  1:178927/152‑1 (MQ=255)
                                                                                                                                                                                     aTGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATAtcctcc                                                                                                                                           >  2:178927/1‑152 (MQ=255)
                                                                                                                                                                                                                 caaGACAGAAACG‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTCTCCCGCCCTTACAAGCTGAAGCGATATCCTCCGTTGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGGGATAGAATGAGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTAc            >  2:418747/2‑251 (MQ=255)
                                                                                                                                                                                                                      cAGAACCG‑ATTTTTGATAAATCAGAGTGATCTGTATTT‑CATTCCGGCACAAGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGAGATATCCTCCGTTGACTTGCTTATGATCAAAGAGGTTAACCTGCGTGAAATCCGCTTCGAGATAGAATGTGTGCTATTTGTTATACTGTGAACTGACTCCCACACCGTTATTCCAGCCATTACCTGTGAAACTGTAATtctc       >  2:52791/1‑251 (MQ=255)
                                                                                                                                                                                                                             g‑ATTTTTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATATCCTCCGTTGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCggg  >  1:138689/1‑249 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
GAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTATTAACGCAGCCCAGGTCACAGGCGGTCGCATTCAGGAAATGTCGGCTTTACTTGGGATCGGCCGCACTACGCTGTGGCGGAAAATGAAGCAACATGGCATTGATGCAGGGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACG‑ATTTCTGATACATCAGAGTGATCTGTATTT‑CATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCGATATCCTCCGTTGACTTGCTTCTGATCAAAGAGGTTACCCTGCGTGTAATCCGCTTCGAGATAGAATGTGTGCTGTTTGTTATACTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGG  >  NC_000913/1252776‑1253245

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