Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,621,627 C→T intergenic (+25/‑142) 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,6270CT95.2% 53.2 / NA 21intergenic (+25/‑142)deoD/yjjJpurine nucleoside phosphorylase/toxin YjjJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (6/14);  minor base G (0/1);  total (6/15)
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.81e-01

GGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTT  >  NC_000913/4621395‑4621701
                                                                                                                                                                                                                                        |                                                                          
ggAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTGCGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTTACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCTCCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGct                                                           <  2:84235/250‑2 (MQ=255)
ggAAAAATACGGCATCCTCGGCGTGGAAAGGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGca                                                           <  2:84217/250‑2 (MQ=255)
      agaCGGCATCCTTGGCGTGGAAAGGGAAGCTGCTGGTATCTACGGCGTCGCTTCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACCTCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGATACt                                                    <  2:480627/249‑1 (MQ=255)
        aCGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGc                                                  >  2:205268/1‑251 (MQ=255)
        aCGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGc                                                  <  1:531811/251‑1 (MQ=255)
            catcctcGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCg                                              >  1:558348/5‑251 (MQ=255)
            catcctcGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGGCGTGTTAAGCCGGAGCAGGAGGCTGCTCCg                                              <  1:470496/247‑1 (MQ=255)
                  cGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCtttt                                        >  1:376534/1‑251 (MQ=255)
                  cGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCttt                                         <  1:205268/250‑1 (MQ=255)
                  cGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCtt                                          <  2:321283/249‑1 (MQ=255)
                   ggCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCttttt                                       <  2:291151/251‑1 (MQ=255)
                        ggAAATGGAAGCGGCTGGGATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAACAGGAAACTGCTCCGGCTTTTTAGTAt                                  >  1:533979/1‑251 (MQ=255)
                           aaTGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTa                               <  1:117872/251‑1 (MQ=255)
                                        ggTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagc                   <  1:459970/250‑8 (MQ=255)
                                         gTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagcaa                 <  1:506076/251‑10 (MQ=255)
                                               aCGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagcaactgatt           <  2:188640/251‑16 (MQ=255)
                                                    gTCGCTGCAGAGTTTGGCGCGAAAGCCCTGAACATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCCCTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagcaactgatttat        >  2:396076/1‑231 (MQ=255)
                                                         tGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagcaactgatttata       <  2:176341/245‑20 (MQ=255)
                                                         tGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCAtcaatagcaactgatttata       >  1:176341/1‑226 (MQ=255)
                                                          gCAGAGTTTGGCGCGCAAGCCCTGACCATCTGCACCGTGTCGGACCACATCCGCTCTCACTAGCAGACCACTGCTGCTGAGCGGCATACCACCGTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGTGGAAGTGTGTTTCTCTGAAAGGCGAACGTCGTGTTAAGCCGGAGCAGGAGACTGCTCCGGCGTTTTAGGATCTATTCAtcaatagcaactgatttatggattt  <  2:585523/249‑25 (MQ=255)
                                                                                                                         cAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACGACAGGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGAAAGGCTAACGTCGTGTTAAGCCGTAGCAGGAGACTGCTCCGGCTTTTTAGTATCTATTCagcaatagcaactgatt           <  1:396947/177‑17 (MQ=255)
                                                                                                                                                                                                                                        |                                                                          
GGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACGACATGATCAAAATCGCACTGGAATCCGTTCTGCTGGGCGATAAAGAGTAATTGTGTTTCGCTGCAAGGCGATTGCCTTGTGAAGCCGGAGCGGGAGACTGCTCCGGCTTTTTAGTATCTATTCATTTTTCTCTCCAGCTTGAATATTTT  >  NC_000913/4621395‑4621701

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