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.87.5% 43.6 / 4.6 16coding (1251/1401 nt)xanQxanthine:H(+) symporter XanQ
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (2/12);  total (3/13)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

ACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTAC  >  NC_000913/3025358‑3025827
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                  
aCGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTAAAATAtt                                                                                                                                                                                                                            <  1:95210/251‑1 (MQ=255)
    tATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTAAAATATTGcc                                                                                                                                                                                                                         >  1:265628/1‑250 (MQ=255)
     aTGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTAAAATATTGCCAg                                                                                                                                                                                                                       <  1:233726/251‑1 (MQ=255)
                aCCATCGCGGTAATGCTGGTTATCCTCTGCTTATTTCCGATGAGTGGCGGCTTCGTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTAAAATATTGCCAGCCTCTATTTAt                                                                                                                                                                                                            <  2:265628/251‑1 (MQ=255)
                   aTCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGCGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGCTTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGCAATTTTTa                                                                                                                                                                                                                                   >  1:27449/1‑226 (MQ=255)
                                 ggTTATCCCTGGCTTATTTCCGAGGAGTGGCTGCTTCTTCACGGCCATGCCCTCGGCTGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGTGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACTCTTATTGTCGCAACTTCGTTATGTTTAAGGCTTGGCGTCTCCTACGATCCCCAAATTTTAAAATATTGCCAGCCTCTATTTATGTATTCGTTGaataca                                                                                                                                                                                            <  2:262188/250‑5 (MQ=255)
                                                    ccGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTa                                                                                                                                                                        <  1:288059/251‑1 (MQ=255)
                                                            gtgcgttttcttcACGTGCATTCCCTCGGGATTTCTGGGAGGCGCAAGGACGTTGATGGTTTCCGTGATTGCAACCGCAGGGATTCGCTTCATCATCACAAACTGTTTATAGCGTCGTGAAACAGTTATTGTTGCTACTTCTTTAGTTTTAGGGCTTGTCTTCTCCTACGATCCCGAAATTGTAAATTATTGCCAGCCTCTATTTATGTATTAGTTGAAAGCCCTATTTTTGCTGGCGGGTCAATTGCGAt                                                                                                                                                                <  2:172424/244‑1 (MQ=255)
                                                                               aTTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTAACTGCGATTTTATTAAATATTAtcctc                                                                                                                                             <  1:366199/251‑1 (MQ=255)
                                                                                                                  gATGTTTTCCATGATTGCCATCGCAGGGAGTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTgg                                                                                                          <  1:215906/251‑1 (MQ=255)
                                                                                                                                     tcgcgcttccgatctatcatcttcTCCAACTGTTTTAAGCGACGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATATTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTAACTGCGATTTTATTAAATATTATCCTCCCGGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAaga                                                                                        <  2:94736/235‑1 (MQ=255)
                                                                                                                                                    atcatcatcACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTTACTGCGGTTTTATTTAATATTATTCCCCCTGGGGGCTACCGACAGGCAAACGTTCTGCCTGGTATTACCTCAGCGGgaga                                                                                        >  1:94736/1‑235 (MQ=255)
                                                                                                                                                                    gTTTAAAGCGCCGGGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCGATTTGTGCTGGCGGGTTAACTGCGATTTTATTAAATATTAtcctc                                                                                                                                             <  1:303018/166‑1 (MQ=255)
                                                                                                                                                                                   aatcacTTATTGTGGCCACTTTTTTAGGTGTAGGGCCTGGCGTCTCCTACGATCCCGAAATTTTAAAATATTGCCAGCCTCTATTTATGCATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTATCTGCGATTTTATTAAATATTATCCTCCCGGGTGGCTACCGACAGGAAATCGTTCTGTCTGGTATTACCTTAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGtt                                         <  1:369693/248‑1 (MQ=255)
                                                                                                                                                                                                           taggtttaggGCTTGGCGTCTCCTACGATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGtttt                   <  1:167477/250‑1 (MQ=255)
                                                                                                                                                                                                                          gcgtcactgaCGATCCAGAAATTTAAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTTCTGGCGGGTTAACTGCGATTGTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATTGATTAACAGTAAAGTAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTTATGTAACCCGTaa  <  2:243586/242‑2 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                  
ACGTTATGTCGGGCGAACCATCGCGGTAATGCTGGTTATCCTCGGCTTATTTCCGATGATTGGCGGCTTCTTCACGACCATTCCCTCGGCAGTTCTGGGAGGCGCAATGACGTTGATGTTTTCCATGATTGCCATCGCAGGGATTCGCATCATCATCACCAACGGTTTAAAGCGCCGTGAAACACTTATTGTCGCCACTTCTTTAGGTTTAGGGCTTGGCGTCTCCTACGATCCCGAAATTTTTAAAATATTGCCAGCCTCTATTTATGTATTAGTTGAAAACCCTATTTGTGCTGGCGGGTTAACTGCGATTTTATTAAATATTATCCTCCCTGGTGGCTACCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTAC  >  NC_000913/3025358‑3025827

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