Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,863 C→T pseudogene (19/48 nt) rph ← 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,8630CT100.0% 135.8 / NA 43pseudogene (19/48 nt)rphribonuclease PH (defective),enzyme, Degradation of RNA, RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/24);  total (19/24)

GCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGT  >  NC_000913/3815734‑3815976
                                                                                                                                  |                                                                                                                 
gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGga                                                                                                          <  2:614484/139‑1 (MQ=255)
         aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTcc                                                                                                 >  2:1107265/1‑139 (MQ=255)
                ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCggg                                                                                          >  1:638628/1‑139 (MQ=255)
                 gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                         <  1:760873/139‑1 (MQ=255)
                     tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCaga                                                                                     <  2:371816/139‑1 (MQ=255)
                       atatGGTTTCATGCCTTCGCTCCTCATCTCACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGc                                                                                   <  1:752063/139‑1 (MQ=255)
                           ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaac                                                                               <  2:638628/139‑1 (MQ=255)
                              ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAg                                                                            <  2:517627/139‑1 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGact                                                                         >  1:603931/1‑137 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                                                         >  1:952371/1‑139 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                                                         >  2:856529/1‑139 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                                                         <  2:952371/139‑1 (MQ=255)
                                      ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                    >  2:307517/1‑139 (MQ=255)
                                           tcctcATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCaga                                                                                     <  2:181081/117‑1 (MQ=255)
                                           tcctcATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCaga                                                                                     >  1:181081/1‑117 (MQ=255)
                                           tcctcATCTTACTTTTCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                                                               <  1:686946/139‑1 (MQ=255)
                                                         ttCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGcc                                                                                        >  1:429901/1‑100 (MQ=255)
                                                         ttCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGcc                                                                                        <  2:429901/100‑1 (MQ=255)
                                                          tCTACAGACAAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGc                                                 >  2:305946/1‑139 (MQ=255)
                                                          tCTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                >  1:972649/1‑139 (MQ=255)
                                                           cTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                              >  1:718188/1‑92 (MQ=255)
                                                           cTACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                              <  2:718188/92‑1 (MQ=255)
                                                            tACAGACAAAAAAAAGGTGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg                                              >  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)
                                                                           aggtgACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                <  1:809469/135‑1 (MQ=255)
                                                                              tgACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                             <  2:998792/138‑1 (MQ=255)
                                                                                aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcga                                                                                                     <  2:231012/65‑1 (MQ=255)
                                                                                aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcga                                                                                                     >  1:231012/1‑65 (MQ=255)
                                                                                            ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgc               >  1:427766/1‑139 (MQ=255)
                                                                                            ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGACTCCCCTCGGGCCAGAGCCAACAAGATGCGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCATCCCCTGCACTGCAATGATgc               >  2:656911/1‑139 (MQ=255)
                                                                                             cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                    <  1:195864/85‑1 (MQ=255)
                                                                                             cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                    >  2:195864/1‑85 (MQ=255)
                                                                                              ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACa                                                                                                               <  2:490475/41‑1 (MQ=255)
                                                                                              ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACa                                                                                                               >  1:490475/1‑41 (MQ=255)
                                                                                              ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcgc             <  1:1107265/139‑1 (MQ=255)
                                                                                                  aaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGt         <  1:305946/139‑1 (MQ=255)
                                                                                                     tCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCtt      <  1:856529/139‑1 (MQ=255)
                                                                                                      cAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc     <  1:885270/139‑1 (MQ=255)
                                                                                                         tttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt  <  2:970793/139‑1 (MQ=255)
                                                                                                                      ccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                       >  1:404835/1‑89 (MQ=255)
                                                                                                                      ccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                       <  2:404835/89‑1 (MQ=255)
                                                                                                                                  |                                                                                                                 
GCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGT  >  NC_000913/3815734‑3815976

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