Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,810 G→T intergenic (‑42/+24) 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,8100GT96.3% 87.9 / ‑3.5 27intergenic (‑42/+24)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (14/12);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCA  >  NC_000913/3815664‑3815951
                                                                                                                                                  |                                                                                                                                             
gaagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCg                                                                                                                                             >  5:7886/1‑149 (MQ=255)
  agaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGAc                                                                                                                                           >  4:143941/1‑149 (MQ=255)
  agaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAACAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGAc                                                                                                                                           <  3:139487/149‑1 (MQ=255)
          gggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGt                                                                                                                                   >  5:103189/1‑149 (MQ=255)
                                       cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGcgccgc                                                                                                      <  6:103189/149‑1 (MQ=255)
                                                       acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACa                                                                                      <  3:143941/149‑1 (MQ=255)
                                                                 tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga                                                                            >  6:141063/1‑149 (MQ=255)
                                                                 tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                              >  6:162403/1‑147 (MQ=255)
                                                                 tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                              <  5:162403/147‑1 (MQ=255)
                                                                            tCAATAAACTGGCGCTTATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACaa                                                                                     <  1:92443/129‑1 (MQ=255)
                                                                            tCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACaa                                                                                     >  2:92443/1‑129 (MQ=255)
                                                                               aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAt                                                                                    >  7:153327/1‑127 (MQ=255)
                                                                               aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATAAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAt                                                                                    <  8:153327/127‑1 (MQ=255)
                                                                                         gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCt                                                                                                 >  2:57069/1‑104 (MQ=255)
                                                                                          cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAg                                                   <  6:7886/149‑1 (MQ=255)
                                                                                                  gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                           >  5:53919/1‑149 (MQ=255)
                                                                                                   tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTAAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAAGGATTCCATTCCCCTCGCGCCAGAGCCAACAAGATaagaag                                           >  6:114655/1‑148 (MQ=255)
                                                                                                        tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                     <  2:172137/149‑1 (MQ=255)
                                                                                                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGctc                                        >  7:27857/1‑143 (MQ=255)
                                                                                                           cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGctc                                        <  8:27857/143‑1 (MQ=255)
                                                                                                                   ctcATCTTACTTTTCTACAGACAAAAAAAAGTGGACTCATCAGTCTCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAAc                           <  5:114655/148‑1 (MQ=255)
                                                                                                                   ctcATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAGTGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACg                          >  4:21852/1‑149 (MQ=255)
                                                                                                                           aCTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGc                    <  5:141063/147‑1 (MQ=255)
                                                                                                                               ttCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg              >  1:192705/1‑149 (MQ=255)
                                                                                                                                           ttttaaGGCGACTGATGAGTCGACTTAAGAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa  >  5:57573/5‑149 (MQ=255)
                                                                                                                                             aaaaGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa  >  5:54785/1‑147 (MQ=255)
                                                                                                                                             aaaaGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa  <  6:54785/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCA  >  NC_000913/3815664‑3815951

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