Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,814 C→T intergenic (‑46/+49) 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,8140CT97.2% 122.8 / ‑3.9 36intergenic (‑46/+49)pyrE/rphorotate phosphoribosyltransferase/truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (22/13);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

ATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGG  >  NC_000913/3815738‑3815888
                                                                            |                                                                          
aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATtcat                                                                        <  1:737968/81‑1 (MQ=255)
   cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAg                                                                     >  1:738002/1‑81 (MQ=255)
     aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTc                                                                   <  1:737969/81‑1 (MQ=255)
        aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGcc                                                                >  2:738001/1‑81 (MQ=255)
        aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGcc                                                                >  1:738020/1‑81 (MQ=255)
            ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTa                                                             >  2:738007/1‑80 (MQ=255)
            ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTa                                                             >  2:738003/1‑80 (MQ=255)
                 tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAAt                                                        >  2:738009/1‑80 (MQ=255)
                    tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGt                                                    >  2:738026/1‑81 (MQ=255)
                       ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTg                                                 >  2:738016/1‑81 (MQ=255)
                        gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGc                                                <  2:737970/81‑1 (MQ=255)
                         tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGc                                                >  2:738027/1‑80 (MQ=255)
                          ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGc                                                >  1:738012/1‑79 (MQ=255)
                           tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAg                                             >  2:738029/1‑81 (MQ=255)
                             aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAATCGCCTTAAAAATCAGTTTGCCAgc                                            <  2:737971/80‑1 (MQ=255)
                               gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGcgcc                                         >  2:737997/1‑81 (MQ=255)
                                 cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGcgccgc                                       <  1:737972/81‑1 (MQ=255)
                                    cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCgccgcc                                      <  1:737974/79‑1 (MQ=255)
                                    cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCgccgcc                                      <  2:737973/79‑1 (MQ=255)
                                     gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTc                                   >  1:738025/1‑81 (MQ=255)
                                     gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTc                                   >  1:738008/1‑81 (MQ=255)
                                      ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCt                                  <  2:737975/81‑1 (MQ=255)
                                         ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGc                                <  1:737976/80‑1 (MQ=255)
                                          tcATCTTACTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGt                              >  2:738010/1‑81 (MQ=255)
                                            aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTc                             <  1:737977/80‑1 (MQ=255)
                                                 aCTTTTCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCt                          <  1:737978/78‑1 (MQ=255)
                                                     ttCTACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAt                     <  1:737979/79‑1 (MQ=255)
                                                        tACAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                >  1:738023/1‑81 (MQ=255)
                                                          cAGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc               >  1:738019/1‑80 (MQ=255)
                                                           aGACAAAAAAAAGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga             >  1:738011/1‑81 (MQ=255)
                                                                   aaaaGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTc     <  1:737980/81‑1 (MQ=255)
                                                                   aaaaGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTc     >  2:738014/1‑81 (MQ=255)
                                                                   aaaaGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTc     >  2:738015/1‑81 (MQ=255)
                                                                     aaGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCgg   <  1:737981/81‑1 (MQ=255)
                                                                     aaGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCgg   >  1:738028/1‑81 (MQ=255)
                                                                      aGGCGATTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCggg  >  1:738018/1‑81 (MQ=255)
                                                                            |                                                                          
ATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGG  >  NC_000913/3815738‑3815888

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