Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 480,295 2 bp→AA coding (13‑14/219 nt) hha ← hemolysin expression‑modulating protein Hha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913480,2950GA93.8% 46.2 / NA 16P5F (CCT→TTT) hhahemolysin expression‑modulating protein Hha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (5/10);  minor base C (1/0);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
*NC_000913480,2960GA100.0% 47.3 / NA 16P5F (CCT→TTT) hhahemolysin expression‑modulating protein Hha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/10);  total (6/10)

AGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTTGTAATTAAGTGCGAA  >  NC_000913/480068‑480537
                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                 
aGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACTCGCTCCAGCGTGTCAATTGTCTGGCAACGACGGAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCt                                                                                                                                                                                                                             <  1:228401/251‑1 (MQ=255)
       aaaGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTTCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTCTTTTTCTCGATAACACGCTCCAGCGTGTCCATTTTCTGGCAACGACGTAAACTCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTTTGGttt                                                                                                                                                                                                                      >  2:97105/1‑250 (MQ=255)
                                        tACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGtt                                                                                                                                                                                                                       <  1:16699/217‑1 (MQ=255)
                                        tACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGtt                                                                                                                                                                                                                       >  2:16699/1‑217 (MQ=255)
                                                                                   tCGGCGAGGGTGTTATTTGCGGCGGAGAAAAAGACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTATCTTTCGTCGAATTGACAAAACAACGGAATAGTCGATTACCt                                                                                                                                           <  2:234927/250‑1 (MQ=255)
                                                                                                 atttGCGGCTGAGTAAAATCCCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTc                                                                                                                            <  2:232434/248‑1 (MQ=255)
                                                                                                            aatAAAATACCGCCAGTTCATTATCTGAGAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTGGTAATAATTCTAACAAGATAAACTCGCAGGATTCTTTTTCGTCGCATTGATAAAACAACGGAAGAGTCGATTTCCTGACTTCCGCCATTTCTTAAGATCCTg                                                                                                                 <  2:155773/249‑1 (MQ=255)
                                                                                                                                     tGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGaac                                                                                        >  2:121525/1‑251 (MQ=255)
                                                                                                                                                atttatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAgg                                                                              <  2:91524/250‑1 (MQ=255)
                                                                                                                                                  ttatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCg                                                                                                                                                  >  2:72/1‑180 (MQ=255)
                                                                                                                                                            cGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGat                                                                  <  2:294690/250‑1 (MQ=255)
                                                                                                                                                              atacgacGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAACATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCCTAACTACTGAACAACATAAAGGTGTCATCCAGATAtt                                                               >  2:30265/6‑251 (MQ=255)
                                                                                                                                                                aaCACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATtc                                                              <  1:362670/250‑1 (MQ=255)
                                                                                                                                                                                 cAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGa                                             <  1:52394/250‑1 (MQ=255)
                                                                                                                                                                                                     aaaCGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACt                        >  2:86730/1‑251 (MQ=255)
                                                                                                                                                                                                                            cGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTTGTAATTAAGTGCGaa  <  1:294891/250‑1 (MQ=255)
                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                 
AGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTTGTAATTAAGTGCGAA  >  NC_000913/480068‑480537

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