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,2950GA95.2% 65.9 / ‑3.7 21P5F (CCT→TTT) hhahemolysin expression‑modulating protein Hha
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (12/8);  total (13/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 = 9.58e-01
*NC_000913480,2960GA95.2% 66.1 / ‑3.4 21P5F (CCT→TTT) hhahemolysin expression‑modulating protein Hha
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (12/8);  total (13/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 = 5.62e-01

AAACAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTAC‑AGTTT‑ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTTGTAATTAAGTGCGA  >  NC_000913/480054‑480536
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                
aaaCAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTAC‑AGTTT‑ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCTGCTGAGTAAAATACCTCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACTCGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTGTCGTGAAAAATTTTTCg                                                                                                                                                                                                                                           <  2:526709/250‑1 (MQ=255)
    aGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTAC‑AGTTT‑ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGGATTTATTTTTCTCGAAAACACGCTCCAGCGTGTCAATTGTCTGGCAAAGACGTCAACGCATTAAATAATCGGGTTTCGTTAAAAATTTTTCAGACa                                                                                                                                                                                                                                       >  2:422328/1‑250 (MQ=255)
                 tCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTAC‑AGTTT‑ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATg                                                                                                                                                                                                                         >  1:211654/1‑251 (MQ=255)
                        gCGAATTATTTATTAGCGAATAAATTTCCAGACTGAGGAAGGGATCTTGTCGTAC‑AGTTT‑ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCGGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGtaa                                                                                                                                                                                                                  <  1:78421/251‑1 (MQ=255)
                                               aTTTCCATACTGAGGAAGGGATCTTGTCGTAC‑AGTTT‑ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCCACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGc                                                                                                                                                                                           >  1:58788/1‑251 (MQ=255)
                                                 tcccttACTGAGGAAGGGTTTTTGTCG‑ACTAGTTT‑ATTCATGGTCAATTCTGCGTGTCGGTGATCTGCGGCTGATTAAAATTCCGCCAGTTCATTATCTGATAATTCGTATTTATTGTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGTCAACGACGTAAACTCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTAGGGTTGTAATAATTCTAACAAGATAAACTCGCa                                                                                                                                                                                          <  2:463982/245‑1 (MQ=255)
                                                  tcaaTACTGAGGAAGGGATCTTGTCGTA‑‑AGTTTTATTCATGTTCAATTCGGCGAGGGGGTGATCTGCGGTTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAg                                                                                                                                                                                         <  2:527322/247‑1 (MQ=255)
                                                   ccATACGGAGGAAGGGTTTTTGTCGTAC‑ATTTT‑ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAgg                                                                                                                                                                                        <  2:133899/250‑1 (MQ=255)
                                                    cATACTGAGGAAGGGATCTTGTCGTCC‑AGTTT‑CTTCATGGTCAATTCGGCGTGGCGGTGATCTGCGGCTGAGTAAAATACCGCCCGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATACCACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGAt                                                                                                                                                                                      >  2:260274/1‑251 (MQ=255)
                                                           aGGAAGGGATCTTGTCGTAC‑AGTTT‑ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCttt                                                                                                                                                                               >  1:118180/1‑251 (MQ=255)
                                                                            tAC‑AGTTT‑ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAaac                                                                                                                                                              >  1:221936/1‑251 (MQ=255)
                                                                                             gTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTCAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGAt                                                                                                                                               >  1:448343/1‑251 (MQ=255)
                                                                                             gTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGAt                                                                                                                                               >  1:159301/1‑251 (MQ=255)
                                                                                                                                        ccAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATAcc                                                                                                    <  1:17361/251‑1 (MQ=255)
                                                                                                                                            ttCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATa                                                                                                 >  1:150497/1‑250 (MQ=255)
                                                                                                                                                   tCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTg                                                                                          >  2:346636/1‑250 (MQ=255)
                                                                                                                                                      gATAATTCGTATTTAGTTTTCTCGATAACCCGCTCCAGCGGGTCAAGTGTCTGGCAACGACGTAAACGCATTAAATAAGCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACGTCCGCCATTTCGGAAGATCCTGCATATTAATACCATAACTACTGaaca                                                                                      <  2:58788/251‑1 (MQ=255)
                                                                                                                                                                             gATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGatat                                                               >  1:245334/1‑251 (MQ=255)
                                                                                                                                                                               tAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATtc                                                             >  1:339063/1‑251 (MQ=255)
                                                                                                                                                                                        ccAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAAATTTTTCAGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCt                                                     <  1:422328/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cGTTAAAAATTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTTGTAATTAAGTGCGa  >  2:361773/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                
AAACAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTAC‑AGTTT‑ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGATCTGCTCAATGAGCTTATTGTCTTCATTATACTTAATTTTGTAATTAAGTGCGA  >  NC_000913/480054‑480536

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