Predicted mutation
evidence seq id position mutation freq annotation gene description
JC NC_000913_3_bme_pgi 3,815,859 Δ82 bp 67.5% [rph][rph] [rph], [rph]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 381585810 (0.210)21 (0.490) 20/252 NT 67.5% pseudogene (24/48 nt) rph ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
?NC_000913_3_bme_pgi 3815941 = 11 (0.250)pseudogene (609/669 nt) rph ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

AACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGC  >  NC_000913_3_bme_pgi/3815710‑3816003
                                                                                                                                                    |                                                                                                                                                 
aaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGC                                                                                                                                                   <  1:92139‑M1/149‑1 (MQ=255)
               ttgctAAGCGCAAATTCAATAGACTGGCTGTGATATGGTTTCATGCGTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGGGACTCATCAGTCGCCTTAAAAATCAGTTTTCCAGCGCCGCCTTCTGCGTCGCTACAATGgat                                                                                                                                    <  2:344560/149‑1 (MQ=255)
                 gctAAGCCTAAATTCAATAAACTGGCGCTGATATGTTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaa                                                                                                                                  <  2:430567‑M1/149‑18 (MQ=255)
                   taggcgcAAATGCCATGAACTGGCGCTGATATGTTGTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGGCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga                                                                                                                                <  2:110018/146‑1 (MQ=255)
                   tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatg                                                                                                                                >  1:738960‑M1/1‑130 (MQ=255)
                   tAAGCGCAAATTAAATAAACTGGCGCTGATCTGGTGTTTTGCGTGCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga                                                                                                                                <  2:493586/149‑1 (MQ=255)
                    aaGCGCAAATTTATTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatga                                                                                                                               <  1:350297‑M1/149‑21 (MQ=255)
                      gggcaaATTCACTACACTGTCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTTTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCTCCTTCTGCcgtcccctgcacttcaatgatg                                                                                                                             <  2:362881‑M1/147‑23 (MQ=255)
                         cAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgc                                                                                                                          <  2:738960‑M1/149‑26 (MQ=255)
                          aaaTTCAATAAACGGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTATAAAAACAAAAAAAAGGCGACTCATCAGTCGACTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcc                                                                                                                         <  1:367025‑M1/149‑27 (MQ=255)
                                  tAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAACAAAGGCCACTCCTCAGTCGCATTAAAAATCAATTTGCCAGCGCCGCCATCTGCcctcccccacacttccatgatgcccccctctcac                                                                                                                 >  2:154316‑M1/1‑115 (MQ=255)
                                    aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGATAAAAAAAAGGCGACGCATCAGTCGCCTTAAAAATAAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggg                                                                                                               <  2:679901‑M1/149‑37 (MQ=255)
                                            cTGATAGGGTTTCATGCCTTCGGTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                                                                                       <  2:691020‑M1/149‑45 (MQ=255)
                                                  tGGTTTCATGCATTCGCTCCTCATCTTACTGTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcgttcatcactccgttca                                                                                                 >  2:31756‑M1/1‑99 (MQ=255)
                                                   ggTTTCAGGCCTTCGCTCCTGATCTTATTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTATAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTa                                                                                                <  2:132549/149‑1 (MQ=255)
                                                     tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgt                                                                                              >  2:483277‑M1/1‑96 (MQ=255)
                                                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                                                                 >  1:726200‑M1/1‑92 (MQ=255)
                                                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaat                                                                                                                                 <  2:726200‑M1/110‑19 (MQ=255)
                                                             cTTCGCTCCTCATCTTATTTTTCTAAATACAAAAAAAAGCGACTCATCATTCGCCTTAAAAATCAGTTTGCTAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATgg                                                                                     <  1:269133/149‑1 (MQ=255)
                                                                 gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcgg                                                                                  <  2:110289‑M1/149‑66 (MQ=255)
                                                                     cttatctTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTCCCACGCAGTGCGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcaga                                                                              <  1:882902‑M1/146‑70 (MQ=255)
                                                                     ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcctcactacgttcatgtcggtctctgcggcaga                                                                              >  2:365885‑M1/1‑80 (MQ=255)
                                                                        aTCTTACTTTTTTACAGACAAAAAAAAGGCGACTCATCAGTCGCTTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtc                                                                           <  1:374474‑M1/149‑73 (MQ=255)
                                                                                       aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGACAGCGCCGCCTTCTCAGTTGCTACAATGGATACGGTCCCCCTCGGGCCAGAGCCAAAAAAATGAGTATCTCTCCATCGGTAAACGAATCGCCTTCTACCTTcc                                                            >  1:662050/1‑149 (MQ=255)
                                                                                        gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                                           >  1:225886/1‑149 (MQ=255)
                                                                                               aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                                                                                                                >  1:479162‑M1/1‑54 (MQ=255)
                                                                                               aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgt                                                                                                                                                <  2:479162‑M1/57‑4 (MQ=255)
                                                                                                     cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgctt                                              <  1:900629‑M1/149‑102 (MQ=255)
                                                                                                         tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgatctctacgccagagacttcaacgtattccagatcgcaaaccgcttcgcc                                          >  1:697853‑M1/1‑44 (MQ=255)
                                                                                                                    cTTAAAATGCAGTTTGCCAGCGCCCCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcgggctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattc                               <  2:594543‑M1/149‑117 (MQ=255)
                                                                                                                              aGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattacagatcgcaaaccgcttcgccgttcacaattccgacagaaac                     >  1:549649‑M1/1‑23 (MQ=255)
                                                                                                                              aGTTTGCCAGCGCCGCCTTCTGCGTCGCGACAATGGATTCTATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTa                     <  2:708013/149‑1 (MQ=255)
                                                                                                                                tttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                   >  1:33856/1‑149 (MQ=255)
                                                                                                                                      agggccgccTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATg             <  2:432081/146‑1 (MQ=255)
                                                                                                                                       gcgcCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgtattccagatcgcaaaccgcttcgccgttcacaattccgaccaaaactgcggctac            >  2:533539‑M1/1‑14 (MQ=255)
                                                                                                                                           cgccTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctccggcagaatctttactgtaagccagatcgcaaaccgcttcgccgttcacaattccgacagaaactgcggctaccatc        >  1:659310‑M1/1‑10 (MQ=255)
                                                                                                                                                 cTGCGTCGCTACAATGGATTCGATTCCCCTCGTGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc  <  1:60459/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGC  >  NC_000913_3_bme_pgi/3815710‑3816003

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.