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,2950GA86.7% 33.4 / ‑2.1 15P5F (CCT→TTT) hhahemolysin expression‑modulating protein Hha
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (5/8);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
*NC_000913480,2960GA86.7% 31.1 / ‑0.8 15P5F (CCT→TTT) hhahemolysin expression‑modulating protein Hha
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (5/8);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

ACAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTT  >  NC_000913/480056‑480522
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                  
aCAGCCGTTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGAc                                                                                                                                                                                                                          <  1:116797/251‑1 (MQ=255)
                        gAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATAACGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAAAACGCTCCAGCGTGTCAATTGTCTGGAAACGTCGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGGTTTTCGGACATACTTCTACCTATGGTTGtaat                                                                                                                                                                                                   >  2:132122/1‑250 (MQ=255)
                                                                                         gTCATTTCGGCGAGGCGGGGATTAGCGGCTATTTAAAATACCGCCAGTTCATTATCTGATAATTCTTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCATGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGAt                                                                                                                                 <  2:250318/251‑1 (MQ=255)
                                                                                                      ggCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCg                                                                                                                    <  1:212188/251‑1 (MQ=255)
                                                                                                               ggtgggctgagaaaaaacacGCCAGTTCATTAGCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGGTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTg                                                                                                           <  2:231013/232‑1 (MQ=255)
                                                                                                                       gAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTAAAGATCCt                                                                                                   >  2:23171/1‑251 (MQ=255)
                                                                                                                                            ttATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAactac                                                                              >  1:207943/1‑251 (MQ=255)
                                                                                                                                                                                        cGTGTCAATTGTCTGGCAACGACGGAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTGTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATGCGTCGATCTGCTCa                                  <  2:230659/251‑1 (MQ=255)
                                                                                                                                                                                                    cTGGCAACGACGTAAACGCATTCAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGTCTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTg                      >  2:37855/1‑251 (MQ=255)
                                                                                                                                                                                                      ggCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTg                      <  1:26649/249‑1 (MQ=255)
                                                                                                                                                                                                         aaCGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTc                 >  1:166094/1‑251 (MQ=255)
                                                                                                                                                                                                            gacgTAAACGAATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAAATCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCATGCATATTTATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTTTCTTCAtt              >  2:157046/1‑251 (MQ=255)
                                                                                                                                                                                                             acgTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGTCATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTTTCTTCATta             <  1:55006/251‑1 (MQ=255)
                                                                                                                                                                                                                         taaataaaCGGTTTTTGTTTAAAATTTTTCGGACATACTTCTTCTGATGGTTGTAATAATTATAACAATATAAACTCGCAGTATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGGGTCATCCAGATATTCGTCGATCTGCGCAATGAGCTTATTGTCTTCATTATACTTAAtttt  <  2:267273/250‑1 (MQ=255)
                                                                                                                                                                                                                         tAAATAATCGGTTGTCTTGAAAAATTTTTCGGACATTCTTCTCCCTTTGTTTTTAATAATTATAACAAGATAAACTCGCATGATTCTCTTTCGTCGCATTGACTAATCAACGGAATAGTCGATTACCTGACTCCCTCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCAGTATACTTAAtttt  <  1:246434/250‑1 (MQ=255)
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                  
ACAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTT  >  NC_000913/480056‑480522

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