Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,517,532 C→T 18.4% F48F (TTC→TTT ortT → orphan toxin OrtT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,517,5320CT18.4% 47.0 / 3.9 21F48F (TTC→TTTortTorphan toxin OrtT
Reads supporting (aligned to +/- strand):  ref base C (8/9);  new base T (2/2);  total (10/11)
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.44e-01

TTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAA  >  NC_000913/1517414‑1517729
                                                                                                                      |                                                                                                                                                                                                     
ttttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCcgc                                           >  1:133787/1‑275 (MQ=255)
    gtATGCGGTTAGGGCAGCAATCGCATTTCTTATCACCTGGTGTCTTTCTCACGATAAGAAACGCAGCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCa                                                                                                                                                               <  1:68770/154‑1 (MQ=255)
            tttgggccaataccGCATTTTTTTTCCCCGTGTTTATTTCTCACGTTAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAagagag                               <  2:133787/262‑1 (MQ=255)
                    gcaATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGcac                                                                                                                                                >  1:80678/1‑154 (MQ=255)
                     caATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTa                                                                                                                                                                   >  2:324807/1‑134 (MQ=255)
                        tCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACg                                                                                                                                            <  2:80678/154‑1 (MQ=255)
                         cGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCa                                                                                                                                                               <  1:324807/134‑1 (MQ=255)
                                                  tctcACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGAc                                                                                                                                                     >  2:71276/1‑119 (MQ=255)
                                                      aCGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGca                                                                                                                                                 <  1:71276/119‑1 (MQ=255)
                                                              aaaCGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGcac                                                                                                                                                >  1:23713/1‑112 (MQ=255)
                                                                  gCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACg                                                                                                                                            <  2:23713/112‑1 (MQ=255)
                                                                                                      cATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGAc                                                                                                                                       >  2:77930/1‑81 (MQ=255)
                                                                                                          gCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTtat                                                                                                                                   <  1:77930/81‑1 (MQ=255)
                                                                                                           ccGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGt                                                                                    >  1:233355/1‑127 (MQ=255)
                                                                                                           ccGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGaa                                    >  1:288963/1‑175 (MQ=255)
                                                                                                            cGATGAGTTTTCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAatta      >  2:52886/1‑204 (MQ=255)
                                                                                                            cGATGAGTTTTCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTaa  <  1:52881/208‑1 (MQ=255)
                                                                                                            cGATGAGTTTTCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTaa  >  2:52881/1‑208 (MQ=255)
                                                                                                               tGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGt                                                                                <  2:233355/127‑1 (MQ=255)
                                                                                                               tGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAagaga                                <  2:288963/175‑1 (MQ=255)
                                                                                                                gAGTTTTCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTaa  <  1:52886/204‑1 (MQ=255)
                                                                                                                      |                                                                                                                                                                                                     
TTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAA  >  NC_000913/1517414‑1517729

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: