| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | NC_000913 | 2,620,218 | T→A | 21.1% | 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_000913 | 2,620,218 | 0 | T | A | 21.1% | 37.5 / 7.8 | 19 | intergenic (‑58/+28) | uraA/upp | uracil permease/uracil phosphoribosyltransferase | 
| Reads supporting (aligned to +/- strand): ref base T (7/8); new base A (2/2); total (9/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01 | |||||||||||
CAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCG  >  NC_000913/2620073‑2620345                                                                                                                                                 |                                                                                                                               cAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTaaa                                                                                                                              <  1:824487/149‑1 (MQ=255)               aaTGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTaaa                                                                                                                              >  1:1618985/1‑134 (MQ=255)               aaTGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTaaa                                                                                                                              <  2:1618985/134‑1 (MQ=255)                                 ggATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTaaa                                                                                                                              <  4:1569494/116‑1 (MQ=255)                                 ggATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTaaa                                                                                                                              >  3:1569494/1‑116 (MQ=255)                                                         tttCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGc                                                                                                                        >  4:676400/1‑98 (MQ=255)                                                         tttCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTTTAAGAGTCGGc                                                                                                                        <  3:676400/98‑1 (MQ=255)                                                                         aGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCg                                                                                                                          >  1:1858354/1‑80 (MQ=255)                                                                         aGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCg                                                                                                                          <  2:1858354/80‑1 (MQ=255)                                                                                                           tgtgTTATAGGGGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTc                                                                                                                           >  1:1394302/1‑45 (MQ=255)                                                                                                           tgtgTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTc                                                                                                                           <  2:1394302/45‑1 (MQ=255)                                                                                                                   aGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTAttt                                                                                                             >  2:394157/1‑51 (MQ=255)                                                                                                                   aGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTAttt                                                                                                             <  1:394157/51‑1 (MQ=255)                                                                                                                     gcgcTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTc                                                                                               <  3:926420/63‑1 (MQ=255)                                                                                                                     gcgcTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTc                                                                                               >  4:926420/1‑63 (MQ=255)                                                                                                                                   ttaaaGCCGACTTTAAGAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGtt    >  3:1435127/3‑140 (MQ=255)                                                                                                                                   ttaaaGCCGACTTTAAGAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGtt    <  4:1435127/138‑1 (MQ=255)                                                                                                                                           gACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCg  >  3:1284494/1‑134 (MQ=255)                                                                                                                                           gACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCg  <  4:1284494/134‑1 (MQ=255)                                                                                                                                                ttAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCgg             >  1:1030898/1‑118 (MQ=255)                                                                                                                                                ttAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCgg             <  2:1030898/118‑1 (MQ=255)                                                                                                                                                 |                                                                                                                               CAAACATGGCGAACAAATGTTGCAAACTAAGCGGGATTGTCTGTAAAAGTGGCGGTCTTTCACTCACCCCGATAGCACGGCGCGTCATAGTATTATCCTCTGTATTATGTGTTATAGGCGCTTTACTCAAAAAAAAGCCGACTCTTAAAGTCGGCTTTAATTATTTTTATTCTTTATTTCGTACCAAAGATTTTGTCACCGGCATCGCCGAGGCCCGGAATAATGTATCCGTGCTCGTTCAGTCCCTGATCAATCGATGCGGTATACAGTTCG  >  NC_000913/2620073‑2620345 | 
| Alignment Legend | 
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG | 
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑ |