Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,947 Δ1 bp pseudogene (603/669 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,9440C.100.0% 202.7 / NA 49pseudogene (606/669 nt)rphribonuclease PH (defective),enzyme, Degradation of RNA, RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (27/22);  total (27/22)

AAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCA  >  NC_000913/3815807‑3816077
                                                                                                                                         |                                                                                                                                     
aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcc                                                                                                                                      <  2:1184413/139‑1 (MQ=255)
aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcc                                                                                                                                      <  1:711423/139‑1 (MQ=255)
   gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccctg                                                                                                                                   <  2:141130/139‑3 (MQ=255)
   gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccctg                                                                                                                                   <  2:336355/139‑3 (MQ=255)
   gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccctg                                                                                                                                   <  2:458598/139‑3 (MQ=255)
          atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCa                                                                                                                           >  2:189245/1‑139 (MQ=255)
                 gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATgc                                                                                                                    >  1:332505/1‑139 (MQ=255)
                             aGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                                                                                        <  1:1243408/139‑1 (MQ=255)
                              gTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                       <  2:717615/139‑1 (MQ=255)
                              gTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                       >  1:1178297/1‑139 (MQ=255)
                               tttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                          >  1:332961/1‑135 (MQ=255)
                               tttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAGGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                          <  2:332961/135‑1 (MQ=255)
                                ttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                                                                                     <  2:1178297/139‑1 (MQ=255)
                                  gCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatc                                                                                                   <  2:143058/139‑1 (MQ=255)
                                    cAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                 <  2:332505/139‑1 (MQ=255)
                                         ccgccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                       >  1:578741/1‑128 (MQ=255)
                                         ccgccTTCTGCGTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                       <  2:578741/128‑1 (MQ=255)
                                                    gTCGTTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtc                                                                                 >  2:874845/1‑139 (MQ=255)
                                                        ttACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATgcgc                                                                                                                  <  1:132880/101‑1 (MQ=255)
                                                        ttACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATgcgc                                                                                                                  >  2:132880/2‑102 (MQ=255)
                                                          aCAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                           <  2:1288209/139‑1 (MQ=255)
                                                             aTGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                                                     >  2:432207/1‑79 (MQ=255)
                                                             aTGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                                                     <  1:432207/79‑1 (MQ=255)
                                                               ggATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                      <  1:189245/139‑1 (MQ=255)
                                                                  ttcgattcCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                                                   >  1:1287160/1‑139 (MQ=255)
                                                                   tcgattcCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccctg                                                                                                                                   <  2:83467/75‑3 (MQ=255)
                                                                   tcgattcCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccctg                                                                                                                                   >  1:83467/1‑73 (MQ=255)
                                                                                 gCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAg                                                    >  2:707277/1‑139 (MQ=255)
                                                                                    agagCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                                 >  2:325749/1‑139 (MQ=255)
                                                                                    agagCCAACAAGATGAGTAGCGCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                                 >  1:324310/1‑139 (MQ=255)
                                                                                      agCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                                                                                        >  2:25055/1‑82 (MQ=255)
                                                                                      agCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCgg                                                                                                        <  1:25055/82‑1 (MQ=255)
                                                                                      agCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                    >  2:1121291/1‑102 (MQ=255)
                                                                                      agCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                    <  1:1121291/102‑1 (MQ=255)
                                                                                           acaaGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                                                             >  2:179038/1‑88 (MQ=255)
                                                                                           acaaGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                                                             <  1:179038/88‑1 (MQ=255)
                                                                                           acaaGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAcc                                          >  2:1278829/1‑139 (MQ=255)
                                                                                               gATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                      >  2:220855/1‑139 (MQ=255)
                                                                                                         tctTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                            >  1:538315/1‑139 (MQ=255)
                                                                                                                  gTGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                   >  2:964474/1‑139 (MQ=255)
                                                                                                                   tGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                       >  2:1287259/1‑54 (MQ=255)
                                                                                                                   tGAACGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                       <  1:1287259/54‑1 (MQ=255)
                                                                                                                     aaCGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                       >  2:89947/1‑52 (MQ=255)
                                                                                                                     aaCGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                       <  1:89947/52‑1 (MQ=255)
                                                                                                                     aaCGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                >  2:147467/1‑59 (MQ=255)
                                                                                                                     aaCGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                >  2:1234494/1‑59 (MQ=255)
                                                                                                                     aaCGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                <  1:1234494/59‑1 (MQ=255)
                                                                                                                     aaCGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                <  1:147467/59‑1 (MQ=255)
                                                                                                                     aaCGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                      <  2:963669/85‑1 (MQ=255)
                                                                                                                     aaCGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                      >  1:963669/1‑85 (MQ=255)
                                                                                                                       cGGCTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGaaa              >  1:398213/1‑139 (MQ=255)
                                                                                                                          cTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                         <  1:1013100/45‑1 (MQ=255)
                                                                                                                          cTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCg                                                                                                         >  2:1013100/1‑45 (MQ=255)
                                                                                                                          cTCGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTg           >  2:1215356/1‑139 (MQ=255)
                                                                                                                            cGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                    >  2:738197/1‑64 (MQ=255)
                                                                                                                            cGCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                    <  1:738197/64‑1 (MQ=255)
                                                                                                                             gCCTTCTGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg        <  1:707277/139‑1 (MQ=255)
                                                                                                                                ttCTGCAGTCCCTGCAGTTCACTAATGAGCCCGTCTTCGGTCATCACTCCGTTGATTTCGGGCTCTACGGCAGAGTCCTCACCGTACTCGAGATCGCGAACCGCTTCGCCGTTAACAATTCCGCCAGAAACTCCGGCTa     <  2:324310/139‑1 (MQ=255)
                                                                                                                                   tGCCGTCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCa  <  2:538315/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
AAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCA  >  NC_000913/3815807‑3816077

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