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,2950GA94.4% 58.0 / NA 18P5F (CCT→TTT) hhahemolysin expression‑modulating protein Hha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/8);  minor base C (1/0);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913480,2960GA100.0% 59.6 / NA 18P5F (CCT→TTT) hhahemolysin expression‑modulating protein Hha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/8);  total (10/8)

ATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATA  >  NC_000913/480066‑480513
                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                         
ataGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCGCCAGCGTGTCAATTTTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTAc                                                                                                                                                                                                       <  1:202618/251‑1 (MQ=255)
                                                ggAAGGGATCTTGTCGTACATTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTtc                                                                                                                                                       <  1:9286/251‑1 (MQ=255)
                                                     ggATCTTGTCGTACAGTTTATTAATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGc                                                                                                                                                  >  1:239167/1‑251 (MQ=255)
                                                       aTCTTGTCGTACAGTTTATTCATGGACAATTCGGCGAGGCGGTGATCTGCGGCTGCGTAAAATACCGCCAGTTCATTATCTGATAATTCGCATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTTGCAACGACGTAAACGCATTAAATAATCGGATTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTTTAACAAAATAAAATAGCAAGATTCTCTTTCGTCGcct                                                                                                                                                >  2:259989/1‑249 (MQ=255)
                                                               gtacattgtATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATCCCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTGCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACaaa                                                                                                                                         <  1:200235/242‑1 (MQ=255)
                                                                    gTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGt                                                                                                                                                                                                <  2:54622/190‑1 (MQ=255)
                                                                    gTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGt                                                                                                                                                                                                >  1:54622/1‑190 (MQ=255)
                                                                      ttATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTATCGATAACCAGCTCCAGCGTGTCAATTGTCTTGCAAAGAAGTAAGCACATTAAATAATCGGATTTCGTTAAACATTTTTCGGACATACTTCTACCTATTGTTGTAATAATTCTATCAAGATAAAATCGAAGGGTTCTCTTTCGTCGCATTGACAAaacaaagg                                                                                                                                  >  2:85567/1‑247 (MQ=255)
                                                                                                                 aaaaTACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATACTCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATATTTCTAACAAGATAAACTCGCAGGCTTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGCCTTCCGCCATTTCTGAAGATCCTGCat                                                                                      >  2:57249/1‑251 (MQ=255)
                                                                                                                                    aTCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAATCGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTATCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTg                                                                   >  1:261790/1‑251 (MQ=255)
                                                                                                                                                   tttatttTTCTCGATAACACGCTAAAGCGTGTCAATTGTCAGGCAACGACGTAAACGCATTAAATACTCGGTTTTCGTTAAAAATTTTTCGGCCATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGgtg                                                    >  2:30482/1‑251 (MQ=255)
                                                                                                                                                               gATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGatat                                        >  1:27143/1‑251 (MQ=255)
                                                                                                                                                                              cGTGTCAATTGTCTGGCAACGACGTAATCGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTGTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCGGACTTCCGCCATTTCTGAAGAGCCTGCATATTAATACCATAACTACTGAACAACATAAATGTGTCATCCAGATAGTCGTCGATCTGCTCa                         <  1:119677/251‑1 (MQ=255)
                                                                                                                                                                                     acgatgctctttcgatctAAACGCATTAAATAATCGGTTTTCTTTAAAAATTGTTCGGACATTCTTCTTCCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTTGTTCGTCGCCTTGACAAAACAACGGAATAGTCGATTACCTGACGTCCGCCATTTCTGAAGATCCTTCATATTAATACCATAACTACTGAACAACATAAAGGTGTCGTCCAGATATTCGTCGATCTGCTCAATGAGCt                  <  2:234647/234‑1 (MQ=255)
                                                                                                                                                                                                     gTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATAGTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATtata  <  1:116528/251‑1 (MQ=255)
                                                                                                                                                                                                       aaaCGCATTAAATACTCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATACTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCt                  >  1:234647/1‑233 (MQ=255)
                                                                                                                                                                                                                   taatCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAactac                                                                     <  1:25723/170‑1 (MQ=255)
                                                                                                                                                                                                                   taatCGGATTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAAATCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAactac                                                                     >  2:25723/1‑170 (MQ=255)
                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                         
ATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATA  >  NC_000913/480066‑480513

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