Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,620,216 C→T 51.0% intergenic (‑56/+30) uraA ← / ← upp uracil permease/uracil phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,620,2160CT51.0% ‑6.2 / 7.4 9intergenic (‑56/+30)uraA/uppuracil permease/uracil phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (2/2);  total (5/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

ACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAAT  >  NC_000913/2620076‑2620326
                                                                                                                                            |                                                                                                               
aCATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTc                                                                                                        >  2:916073/1‑149 (MQ=255)
                                      tGTAAAGGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCTTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGAATTTAAGAGTCGGCt                                                                                                    <  2:73318/115‑1 (MQ=255)
                                      tGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGCt                                                                                                    >  1:73318/1‑115 (MQ=255)
                                                      tttCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTAttt                                                                                          <  3:1142210/109‑1 (MQ=255)
                                                      tttCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTAttt                                                                                          >  4:1142210/1‑109 (MQ=255)
                                                            tcacCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGCttt                                                                                                  <  1:739763/95‑1 (MQ=255)
                                                            tcacCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGCttt                                                                                                  >  2:739763/1‑95 (MQ=39)
                                                                                                            ttATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGAttt                                                               >  3:144334/1‑82 (MQ=255)
                                                                                                            ttATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGAttt                                                               <  4:144334/82‑1 (MQ=255)
                                                                                                                                           tctTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAAt  >  3:1253614/1‑112 (MQ=255)
                                                                                                                                           tctTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAAt  <  4:1253614/112‑1 (MQ=255)
                                                                                                                                            |                                                                                                               
ACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAAT  >  NC_000913/2620076‑2620326

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