Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_pgi 3,815,812 G→T intergenic (‑44/+22) pyrE ← / ← rph orotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_pgi3,815,8120GT93.3% 38.3 / ‑6.6 15intergenic (‑44/+22)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (5/9);  minor base C (0/1);  total (5/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 = 1.00e+00

AGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAA  >  NC_000913_3_vch_pgi/3815666‑3815957
                                                                                                                                                    |                                                                                                                                                 
agaaATAGGTGCTTTCGCGCTCGGAGTTAAGCTTAAACCTCACAAACTTTAACACCTGTATGCTAAGCGCAAATACAATAAACTGGCGCTGATATGGTTTCATTCCTTCGCTCCTGATCTTACTTTTATACAGACAAAAAAAAGGCTAc                                                                                                                                                 <  2:88843/149‑1 (MQ=255)
                                                 tAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATTAGTTTGCCAGCGCCGCCTTCTGCGTccc                                                                                                >  2:180627/1‑147 (MQ=255)
                                                                         ttCAAGAAACTGGCGCTAATCTGGTTTCATGACTTCCCTGCCGATCTTACATGTCTACAGAAAAAAGAAATGCTACTCATCAGTCTTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCgg                                                                        <  2:156718/149‑1 (MQ=255)
                                                                           cAATAATCTGGCTCTGGTAGAGGTTTCGGCCTTCGTCCCTCCTCTTACTTTTACACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGc                                                                      <  2:52477/149‑1 (MQ=255)
                                                                                aaCTGGCGCGGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCTACTCATCAGCCGCCTTAAAAATCAGGTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCagag                                                                 <  1:222988/149‑1 (MQ=255)
                                                                                     gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaa                                                             <  1:259047/149‑1 (MQ=255)
                                                                                     gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCTACTCGTCAGTCGCCTTAAAAATCTGTTTGCCAGCGCCGCCTTCTGCGTCGGTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaac                                                            >  2:201178/1‑149 (MQ=255)
                                                                                     gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCCTCGCCGCCTTCTACGTCGCTACAATGGATTCGCTTCCCCTCGGGCCAGAGCcaac                                                            >  2:201207/1‑149 (MQ=255)
                                                                                                        gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCa                                                                                                                           <  1:3375/68‑1 (MQ=255)
                                                                                                        gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCa                                                                                                                           >  2:3375/1‑68 (MQ=255)
                                                                                                          cTTCGCTCCTCATCTTACTTTTCTACATACAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATg                                        <  1:70471/149‑1 (MQ=255)
                                                                                                             cGCTCCTGATCTTACCTTTCTACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGt                                     <  1:135203/149‑1 (MQ=255)
                                                                                                                 cctcATCTTACTTTTCGACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGTCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGCGCCAACAAGATGAGTAGCTCTTCATGGTTGAAc                                 <  2:224644/149‑1 (MQ=255)
                                                                                                                      tCTTACTTTTCTACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAAGAGTAGATCTTCATGGGTGAACGGCtg                            >  2:78234/1‑148 (MQ=255)
                                                                                                                                                 ggCCACTCATCAGCCGCCTCAAAAATCCGTTTGCCAGCGCCCCCTGCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGACCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCCTCCCCTGCACTTCaa  <  2:156081/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAA  >  NC_000913_3_vch_pgi/3815666‑3815957

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