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

AGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATA‑TGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAA  >  NC_000913_3_vch_pgi/3815666‑3815957
                                                                                                                                                    |                                                                                                                                                 
agaaATAGGTGCTTTCGCGCTCGGAGTTAAGCTTAAACCTCACAAACTTTAACACCTGTATGCTAAGCGCAAATACAATAAACTGGCGCTGATA‑TGGTTTCATTCCTTCGCTCCTGATCTTACTTTTATACAGAC‑AAAAAAAAGGCTAc                                                                                                                                                 <  2:88843/149‑1 (MQ=255)
                                                 tAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATA‑TGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATTAGTTTGCCAGCGCCGCCTTCTGCGTccc                                                                                                >  2:180627/1‑147 (MQ=255)
                                                                         ttCAAGAAACTGGCGCTAATC‑TGGTTTCATGACTTCCCTGCCGATCTTACATGTCTACAGAA‑AAAAGAAATGCTACTCATCAGTCTTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCgg                                                                        <  2:156718/149‑1 (MQ=255)
                                                                           cAATAATCTGGCTCTGGTAGAGGTTTC‑GGCCTTCGTCCCTCCTCTTACTTTTACACAGAC‑AAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGc                                                                      <  2:52477/149‑1 (MQ=255)
                                                                                aaCTGGCGCGGATA‑TGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCTACTCATCAGCCGCCTTAAAAATCAGGTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCagag                                                                 <  1:222988/149‑1 (MQ=255)
                                                                                     gcgcTGATA‑TGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaa                                                             <  1:259047/149‑1 (MQ=255)
                                                                                     gcgcTGATA‑TGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCTACTCGTCAGTCGCCTTAAAAATCTGTTTGCCAGCGCCGCCTTCTGCGTCGGTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaac                                                            >  2:201178/1‑149 (MQ=255)
                                                                                     gcgcTGATA‑TGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCCTCGCCGCCTTCTACGTCGCTACAATGGATTCGCTTCCCCTCGGGCCAGAGCcaac                                                            >  2:201207/1‑149 (MQ=255)
                                                                                                        gCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCa                                                                                                                           <  1:3375/68‑1 (MQ=255)
                                                                                                        gCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCa                                                                                                                           >  2:3375/1‑68 (MQ=255)
                                                                                                          cTTCGCTCCTCATCTTACTTTTCTACATAC‑AAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATg                                        <  1:70471/149‑1 (MQ=255)
                                                                                                             cGCTCCTGATCTTACCTTTCTACAGAC‑AAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGt                                     <  1:135203/149‑1 (MQ=255)
                                                                                                                 cctcATCTTACTTTTCGACAGAC‑AAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGTCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGCGCCAACAAGATGAGTAGCTCTTCATGGTTGAAc                                 <  2:224644/149‑1 (MQ=255)
                                                                                                                      tCTTACTTTTCTACAGAC‑AAAAAAAAGGCTACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAAGAGTAGATCTTCATGGGTGAACGGCtg                            >  2:78234/1‑148 (MQ=255)
                                                                                                                                                 ggCCACTCATCAGCCGCCTCAAAAATCCGTTTGCCAGCGCCCCCTGCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGACCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCCTCCCCTGCACTTCaa  <  2:156081/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATA‑TGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAA  >  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: ‑