Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,872,114 G→T 20.0% V633L (GTG→TTG) ‡ xdhB → xanthine dehydrogenase subunit XdhB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,872,1140GT20.0% 43.5 / 6.6 20V633L (GTG→TTG) ‡xdhBxanthine dehydrogenase subunit XdhB
Reads supporting (aligned to +/- strand):  ref base G (5/11);  new base T (2/2);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

CAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCAT  >  NC_002947/4871991‑4872235
                                                                                                                           |                                                                                                                          
cagGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGAtcg                                                                                                               <  2:169304/136‑1 (MQ=255)
       aCTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCCGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATAc                                                                                                        >  1:210750/1‑136 (MQ=255)
              ccAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGAcc                                                                                                 <  2:292474/136‑1 (MQ=255)
                aGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCgg                                                                                               >  1:166116/1‑136 (MQ=255)
                     tCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGc                                                                                                      <  1:249646/124‑1 (MQ=255)
                     tCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGc                                                                                                      >  2:249646/1‑124 (MQ=255)
                             cAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATg                                                                                  <  2:424773/136‑1 (MQ=255)
                               gCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGtt                                                                                <  2:34102/136‑1 (MQ=255)
                               gCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGtt                                                                                <  1:399389/136‑1 (MQ=255)
                                cACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTg                                                                               >  2:53081/1‑136 (MQ=255)
                                    gggTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTgcgcg                                                                           <  2:78990/136‑1 (MQ=255)
                                              gCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCt                                                                 <  1:114681/136‑1 (MQ=255)
                                                         aTCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTgg                                                      <  1:468863/136‑1 (MQ=255)
                                                                                            gTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGTCCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGTCGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAg                   <  1:161737/136‑1 (MQ=255)
                                                                                            gTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGAGATCGTCGATACGCTGACCGGCGAATACACGATGTTGCGCGCCGACATCCTGCACGACGGGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAg                   <  1:130252/136‑1 (MQ=255)
                                                                                                           ctatGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTc    >  1:392411/4‑136 (MQ=255)
                                                                                                           ctatGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTc    <  1:104055/133‑1 (MQ=255)
                                                                                                           cTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGTTGGTTTc    >  2:52571/1‑136 (MQ=255)
                                                                                                            tatGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTcg   <  1:400930/134‑2 (MQ=255)
                                                                                                             atGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTcgt  >  2:511496/2‑134 (MQ=255)
                                                                                                                           |                                                                                                                          
CAGGCCTACTTCGCCCAGGTGTCGCTGTCCAGCACCGGGTTCTACCGCACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCAT  >  NC_002947/4871991‑4872235

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