Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,811 C→T intergenic (‑43/+23) 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_0009133,815,8110CT97.1% 106.9 / ‑3.4 35intergenic (‑43/+23)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective),enzyme, Degradation of RNA, RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/17);  minor base . (0/1);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCC  >  NC_000913/3815675‑3815941
                                                                                                                                         |                                                                                                                                  
ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGa                                                                                                                                  >  1:614484/1‑139 (MQ=255)
       gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACtcatca                                                                                                                           <  2:18580/139‑1 (MQ=255)
                                   aaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGG‑GACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGcgcc                                                                                              <  1:1106071/139‑1 (MQ=255)
                                    aCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGcgcc                                                                                              >  1:970793/1‑139 (MQ=255)
                                                gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGTGACTCATCAGt                                                                                                                         >  1:292013/1‑99 (MQ=255)
                                                gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGTGACTCATCAGt                                                                                                                         <  2:292013/99‑1 (MQ=255)
                                                           gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGga                                                                       <  2:614484/139‑1 (MQ=255)
                                                                    aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTcc                                                              >  2:1107265/1‑139 (MQ=255)
                                                                           ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCggg                                                       >  1:638628/1‑139 (MQ=255)
                                                                            gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                      <  1:760873/139‑1 (MQ=255)
                                                                                tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCaga                                                  <  2:371816/139‑1 (MQ=255)
                                                                                  atatGGTTTCATGCCTTCGCTCCTCATCTCACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGc                                                <  1:752063/139‑1 (MQ=255)
                                                                                      ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaac                                            <  2:638628/139‑1 (MQ=255)
                                                                                         ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAg                                         <  2:517627/139‑1 (MQ=255)
                                                                                            aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGact                                      >  1:603931/1‑137 (MQ=255)
                                                                                            aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                      >  2:856529/1‑139 (MQ=255)
                                                                                            aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                      <  2:952371/139‑1 (MQ=255)
                                                                                            aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                      >  1:952371/1‑139 (MQ=255)
                                                                                                 ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                 >  2:307517/1‑139 (MQ=255)
                                                                                                      tcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCaga                                                  >  1:181081/1‑117 (MQ=255)
                                                                                                      tcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCaga                                                  <  2:181081/117‑1 (MQ=255)
                                                                                                      tcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                            <  1:686946/139‑1 (MQ=255)
                                                                                                           aTCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTg                                                                                      <  2:486541/76‑1 (MQ=255)
                                                                                                           aTCTTACTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTg                                                                                      >  1:486541/1‑76 (MQ=255)
                                                                                                                aCTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCg                                                                                    <  2:498451/73‑1 (MQ=255)
                                                                                                                aCTTTTCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCg                                                                                    >  1:498451/1‑73 (MQ=255)
                                                                                                                    ttCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGcc                                                     <  2:429901/100‑1 (MQ=255)
                                                                                                                    ttCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGcc                                                     >  1:429901/1‑100 (MQ=255)
                                                                                                                     tCTACAGACAAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGc              >  2:305946/1‑139 (MQ=255)
                                                                                                                     tCTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt             >  1:972649/1‑139 (MQ=255)
                                                                                                                      cTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                           >  1:718188/1‑92 (MQ=255)
                                                                                                                      cTACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                           <  2:718188/92‑1 (MQ=255)
                                                                                                                       tACAGAC‑AAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg           >  2:809469/1‑139 (MQ=255)
                                                                                                                               aaaaaaaaGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgat                                                                      >  1:806969/1‑73 (MQ=255)
                                                                                                                               aaaaaaaaGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgat                                                                      <  2:806969/73‑1 (MQ=255)
                                                                                                                                 aaaaaaGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc  <  1:307517/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
GGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCC  >  NC_000913/3815675‑3815941

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑