Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,620,218 T→A 36.4% intergenic (‑58/+28) 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,2180TA36.4% 5.6 / 8.2 11intergenic (‑58/+28)uraA/uppuracil permease/uracil phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base A (2/2);  total (6/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

ACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCTTTTTC  >  NC_000913/2620076‑2620365
                                                                                                                                              |                                                                                                                                                    
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)
                                                                                                                                             ttAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTcc                                                  <  1:532241/101‑1 (MQ=255)
                                                                                                                                             ttAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTcc                                                  >  2:532241/1‑101 (MQ=255)
                                                                                                                                             ttAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCTTTTTc  >  3:346159/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
ACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCGACGTCCGGGTGCGCTTTTTC  >  NC_000913/2620076‑2620365

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