Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,025,601 (T)5→4 coding (1251/1401 nt) xanQ → xanthine:H(+) symporter XanQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,025,6010T.90.9% 36.1 / ‑2.8 11coding (1251/1401 nt)xanQxanthine:H(+) symporter XanQ
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (4/6);  total (4/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTG  >  NC_000913/3025366‑3025814
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                     
gcGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATTACGTTGATGTTTTCCATGAGTGCCATCGCATGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGACGCCTACGATCCCGAAATTTTTaaatagttgcagcct                                                                                                                                                                                                        <  2:42633/250‑13 (MQ=255)
  gggCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCctct                                                                                                                                                                                                     <  1:75236/251‑1 (MQ=255)
      gAACCATCGCGGTAATTCTGGTTATCCTCTGCTTATTTCCGATGAGTGGCGGCTTCTTCACGTCCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGTGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACCCTTATTGTCGCCACTTCTTTAGTTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTTCCAGCCTCTAttt                                                                                                                                                                                                 <  2:128773/251‑1 (MQ=255)
                            tATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCCCCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTTAAAACCCTa                                                                                                                                                                            <  1:207547/250‑1 (MQ=255)
                                                                            cTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTACCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGAGACGTGAAACACTTATTGTCGGCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGATGGCGGGTTAACTGCTATTTTTTTAAATATTTTCCTTCCtgg                                                                                                                           >  1:256244/1‑251 (MQ=255)
                                                                                     tCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGACTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTAACTGCGATTTTATTAAATATTAGCCTCCCTGGTGGCTAcc                                                                                                                   <  2:22549/250‑1 (MQ=255)
                                                                                                                                         cGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTAGTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGTTTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCGCAGCGGAAGAGATGGATTAACa                                                               <  1:105158/250‑1 (MQ=255)
                                                                                                                                                 catcaCCAACGGTTTAACGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTTACTGCGAATTTATTAAATATTATCCTCCCTGGGGGCTACCGACAGGGAAACGGTCTGTCTTGTATTAACTCAGCGGAAGAGATgggataaaagaaaagggg                                                      >  1:184959/1‑235 (MQ=255)
                                                                                                                                                                                                     ggtttaggGCTTGGAGTCTCCTACGATCCCGCAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCATATTTGTGCTGGCGGGTTTACTGCGATTTTTTTAAATATTATCCTCCCCGGTGGGTACCAGCAGGAAAACGTTCTGCCTTGTTTTACCTCAGCGGGAGGGATGGATTACCAGTAAAGGAGTCAATTAAGTCAGGAGAACACACGTTAAAACCTGGACGAGGCAGTTTTAttt  >  2:178988/1‑250 (MQ=255)
                                                                                                                                                                                                                                cccGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTATCTGCGATTTTTTTAAAAATTTTCCTCCCTTGGGGCTACCGACAg                                                                                                              >  2:68860/1‑116 (MQ=255)
                                                                                                                                                                                                                                cccGAAATTTT‑AAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAg                                                                                                              <  1:68860/116‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                     
TCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTG  >  NC_000913/3025366‑3025814

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