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% 173.0 / NA 53pseudogene (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 (26/27);  total (26/27)

GCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAC  >  NC_000913/3815734‑3815982
                                                                                                                                 |                                                                                                                       
gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGga                                                                                                                >  1:1232506/1‑139 (MQ=255)
  aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgatt                                                                                                              <  2:703247/139‑1 (MQ=255)
         aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTcc                                                                                                       >  1:421284/1‑139 (MQ=255)
                 gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                               <  2:157157/139‑1 (MQ=255)
                  cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                    >  1:510874/1‑133 (MQ=255)
                  cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                    <  2:510874/133‑1 (MQ=255)
                        tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGcc                                                                                        <  1:1396227/139‑1 (MQ=255)
                             tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTCCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaacaa                                                                                   >  1:836264/1‑139 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                    >  2:1359176/1‑118 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                    <  1:1359176/118‑1 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                               >  2:266448/1‑123 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                               <  1:266448/123‑1 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                               >  1:1332282/1‑123 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGc                                                                                               <  2:1332282/123‑1 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaaca                                                                                    >  2:6425/1‑134 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaaca                                                                                    <  1:6425/134‑1 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                                                               >  1:72677/1‑139 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                                                               >  2:473523/1‑139 (MQ=255)
                                 aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                                                               >  2:879385/1‑139 (MQ=255)
                                   gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAg                                                                             <  2:1232506/139‑1 (MQ=255)
                                      ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                    <  1:1044477/113‑1 (MQ=255)
                                      ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCt                                                                                                    >  2:1044477/1‑113 (MQ=255)
                                      ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                          <  1:626809/139‑1 (MQ=255)
                                      ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                          <  1:879385/139‑1 (MQ=255)
                                        cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGTGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGct                                                                        <  1:618185/139‑1 (MQ=255)
                                        cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgatt                                                                                                         <  2:1580487/106‑1 (MQ=255)
                                        cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgatt                                                                                                         >  1:1580487/1‑106 (MQ=255)
                                                         ttCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATg                                                                                                                  <  2:1268602/80‑1 (MQ=255)
                                                         ttCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATg                                                                                                                  >  1:1268602/1‑80 (MQ=255)
                                                         ttCTACAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGc                                                       >  1:891681/1‑139 (MQ=255)
                                                             aCAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGc                                                   <  2:891681/139‑1 (MQ=255)
                                                              cAGACAAAAAAAAGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                                  <  2:836264/139‑1 (MQ=255)
                                                                   aaaaaaaaGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                             >  1:847442/1‑139 (MQ=255)
                                                                    aaaaaaaGGCGGCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGc                                            >  1:1228097/1‑139 (MQ=255)
                                                                           ggcggcTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgat                                                                                                          >  2:531180/1‑70 (MQ=255)
                                                                           ggcggcTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgat                                                                                                          <  1:531180/70‑1 (MQ=255)
                                                                           ggcggcTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGcc                                                                                              >  2:1597408/1‑82 (MQ=255)
                                                                           ggcggcTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGcc                                                                                              <  1:1597408/82‑1 (MQ=255)
                                                                           ggcggcTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc                                     >  1:349533/1‑139 (MQ=255)
                                                                                      aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAat                          >  2:1140939/1‑139 (MQ=255)
                                                                                       gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatg                         <  2:847442/139‑1 (MQ=255)
                                                                                          gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgatt                                                                                                         <  1:1115126/56‑1 (MQ=255)
                                                                                          gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgatt                                                                                                         >  2:1115126/1‑56 (MQ=255)
                                                                                          gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                                                            <  2:569983/101‑1 (MQ=255)
                                                                                          gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                                                            >  1:569983/1‑101 (MQ=255)
                                                                                             ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAg                                                                             >  2:1262275/1‑81 (MQ=255)
                                                                                             ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAg                                                                             <  1:1262275/81‑1 (MQ=255)
                                                                                             ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                   >  2:400690/1‑123 (MQ=255)
                                                                                             ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                   <  1:400690/123‑1 (MQ=255)
                                                                                                      aGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgat                                                                                                          >  1:996764/1‑43 (MQ=255)
                                                                                                      aGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGgattcgat                                                                                                          <  2:996764/43‑1 (MQ=255)
                                                                                                       gTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCgg         <  2:1228097/139‑1 (MQ=255)
                                                                                                              aGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc  <  1:413513/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                       
GCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAC  >  NC_000913/3815734‑3815982

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