Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,621,644 G→A intergenic (+42/‑125) deoD → / → yjjJ purine nucleoside phosphorylase/toxin YjjJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,621,6440GA84.6% 24.7 / ‑0.9 14intergenic (+42/‑125)deoD/yjjJpurine nucleoside phosphorylase/toxin YjjJ
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (3/8);  total (5/9)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.67e-01

AAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGT  >  NC_000913/4621397‑4621843
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                       
aaaaaatcgggctcctcGGCGTGGGAAATGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGga                                                                                                                                                                                                      <  1:605187/237‑1 (MQ=255)
 aaaaTACGGCAGCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTGTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGgag                                                                                                                                                                                                     <  1:342470/251‑1 (MQ=255)
             cctcGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCt                                                                                                                                                                                         >  1:22900/2‑251 (MQ=255)
                   ggTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCAGCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAg                                                                                                                                                                                   <  2:71038/250‑1 (MQ=255)
                       gAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATc                                                                                                                                                                               <  1:114074/251‑1 (MQ=255)
                         aaaTGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCt                                                                                                                                                                              >  1:673515/1‑251 (MQ=255)
                          aaTGGAAGCGGCTGGTATCTCCGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTa                                                                                                                                                                             <  2:122353/251‑1 (MQ=255)
                             ggAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTAtt                                                                                                                                                                           <  1:498748/250‑1 (MQ=255)
                                       ggTATCTACGGTGTCGCTTCAGGGTTTGGCGCGAATGCCCGGACCATCTGCTCCTGGGCTGACCACATCCGCCCTCACGCGCAGACCCCTGTCGCTGATCGTCATACCACTTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGATATGTGAACGTCGTGTTAAGCCGGAGCCGGGGACTGCTCCGGCTTTTTAGTAGCTATTcttcaatatc                                                                                                                                                                 <  2:535636/250‑10 (MQ=255)
                                          aTCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagcaac                                                                                                                                                              <  1:444011/250‑11 (MQ=255)
                                              aCGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagcaactga                                                                                                                                                           <  1:679443/249‑14 (MQ=255)
                                                ggCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagcaactgattt                                                                                                                                                        >  1:384982/1‑234 (MQ=255)
                                                                                                      cacaTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGACGATTGCCTTGTGAAGCCGGAGCGGGCGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGACTATTTTCGCTATACTTTTCAGTGAAATGTGTTCATAAATCTATTCAAGTATCt                                                                                                 >  2:331008/1‑251 (MQ=255)
                                                                                                                                                                                                     tAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAAtcagt  >  1:87777/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                       
AAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTTCGCTATACTTTTCAGTGAAATGTGTTAATAAATCTATTCAAGTATCTATTCACGAATCTATTCATTAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGT  >  NC_000913/4621397‑4621843

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