Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,828 T→G intergenic (‑60/+35) pyrE ← / ← rph orotate phosphoribosyltransferase/truncated RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,8280TG100.0% 59.7 / NA 19intergenic (‑60/+35)pyrE/rphorotate phosphoribosyltransferase/truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (6/13);  total (6/13)

TTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAG  >  NC_000913/3815763‑3815892
                                                                 |                                                                
tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGaa                                                                >  1:1153607/1‑68 (MQ=255)
tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGaa                                                                >  1:1305874/1‑68 (MQ=255)
tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGaa                                                                >  1:2375359/1‑68 (MQ=255)
    aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAAt                                                            <  1:423062/68‑1 (MQ=255)
     tGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATc                                                           >  1:2080635/1‑68 (MQ=255)
      gCCTTCGCTCCTCATCTTAGTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCa                                                          <  1:2963357/68‑1 (MQ=255)
      gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCa                                                          >  1:463416/1‑68 (MQ=255)
              tcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCa                                                  >  1:72926/1‑68 (MQ=255)
               cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAg                                                 <  1:536072/68‑1 (MQ=255)
                       tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCt                                         <  1:2816597/68‑1 (MQ=255)
                       tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCt                                         <  1:407095/68‑1 (MQ=255)
                             tCTACAGACAAAAAAAAGCGGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCg                                   <  1:3234275/68‑1 (MQ=255)
                                aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCg                                <  1:1711096/68‑1 (MQ=255)
                                    aCAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAc                            <  1:421933/68‑1 (MQ=255)
                                     cAAAAAAAAGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACa                           <  1:3352195/68‑1 (MQ=255)
                                       aaaaaaaGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAt                         <  1:1558287/68‑1 (MQ=255)
                                         aaaaaGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATgg                       <  1:1129132/68‑1 (MQ=255)
                                         aaaaaGGCGACTCATCAGTCGCCTGAAAAATCAGTGTGCCAGCGCCGCCTTCTGCGTCGCTACAATgg                       <  1:3356410/68‑1 (MQ=255)
                                            aaGGCGACTCATCAGTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                    >  1:3098039/1‑68 (MQ=255)
                                                     catcaCTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTcccc           <  1:909207/68‑1 (MQ=255)
                                                          gTCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCggg      <  1:123226/68‑1 (MQ=255)
                                                           tCGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGc     <  1:2658148/68‑1 (MQ=255)
                                                            cGCCTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGcc    >  1:2249306/1‑68 (MQ=255)
                                                              ccTGAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCag  >  1:1096359/1‑68 (MQ=255)
                                                                 |                                                                
TTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAG  >  NC_000913/3815763‑3815892

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: