Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 410,666 C→T P175S (CCT→TCT)  mak → fructokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913410,6660CT85.7% 24.2 / ‑2.7 14P175S (CCT→TCT) makfructokinase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (3/9);  total (4/10)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCCCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGG  >  NC_000913/410434‑410879
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                     
gNGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTaaa                                                                                                                                                                                                      >  1:278067/1‑250 (MQ=255)
                                               aaGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGTTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGGGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATTGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTAGTTCGGTAAACAAGGTTGTATTGAAACCTTTATGTCAGTCACGGGATTCGCGATGGAt                                                                                                                                                      <  1:293088/251‑1 (MQ=255)
                                                                            cAGACGTTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGTCATTCAATGGTCGGGCGCATATCGGCGGCAATGGCACGTCAGGTGAGTGGGGCCACAATCCGCTACCGTGGATGGACGAAGACGAACTGCTTTATCGCGAGCAAGTCTCTTGTTATTTCTGGAAACAAGGTTGGATTGAAACCTTTTTTTCGGGCACGGGATTCGCGATTGATTATCGTCGTTTGAGc                                                                                                                                       <  1:320599/237‑1 (MQ=255)
                                                                              gaaggTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAg                                                                                                                       <  1:353082/248‑1 (MQ=255)
                                                                                                          ggATGCGGCGCGGGCGTGGCATTCAATGTGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTAGTGCGGTAAACAAGGTTGTATTGAAGCCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGttt                                                                                                                                           <  1:338222/203‑1 (MQ=255)
                                                                                                                    ggggCGTGGCATTCAGTGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGGTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCCCTTGGTATGGCGGTAAACAAGGGTGTATTGAAACCTTTATTTCGGGGACGTGAGGCGCGATGGATTATCATCTTGTGAGCGGACATGCGCTGAAAGGCAGTGAAAGTATCCGCCTGGTTGAAGAAAGCGATCCg                                                                                 <  1:237390/250‑1 (MQ=255)
                                                                                                                                     gggcgggcGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGTGCACGGTATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTTGTTTAAGAAAGCGAGCCGGTAGCTGAACTGGCAtt                                                                <  1:127904/251‑1 (MQ=255)
                                                                                                                                        cgggcgCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAGCCTTGATGTCGGTTACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAg                                                                                                                       <  1:313900/193‑1 (MQ=255)
                                                                                                                                            cgcATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACCCAATCCGCCACCGTGCAAGGCAGCACAACAACCGCGTTATCGCGAGGAAGTACATTACTATTGCGGTAAACAAGGTTCTATTGTCACCTTTATTTTGGGCACGGGATTCCCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGGAGCGGAACTGGCATTGCggcgg                                                         >  1:188334/1‑247 (MQ=255)
                                                                                                                                                                 cGGCAGGGGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGGTATTGCGGTAAACAAGGGTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTg                                                                                                                                          <  1:163518/149‑1 (MQ=255)
                                                                                                                                                                               gnACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGc                                                                   >  1:277077/1‑206 (MQ=255)
                                                                                                                                                                                                   natggaCGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCgg  >  1:64404/2‑251 (MQ=255)
                                                                                                                                                                                                                      cGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTAGTTCGGGCACGGGATTCGCTATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGc                                                                            <  1:229061/158‑1 (MQ=255)
                                                                                                                                                                                                                       gTTATCGCGAGGAAGGCTCTTGTTAGTGCGGTAAACAAGGTTTTATTGAAACCGTTATGTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCt                                                                                                                            <  1:306077/109‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                     
GCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCCCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGG  >  NC_000913/410434‑410879

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