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

AGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATC  >  NC_000913_3_bme_pgi/3815797‑3816079
                                                                                                                                                |                                                                                                                                          
aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGACAGCGCCGCCTTCTCAGTTGCTACAATGGATACGGTCCCCCTCGGGCCAGAGCCAAAAAAATGAGTATCTCTCCATCGGTAAACGAATCGCCTTCTACCTTcc                                                                                                                                        >  1:662050/1‑149 (MQ=255)
 gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                                                       >  1:225886/1‑149 (MQ=255)
            gctaagcctaaattcaataaactggcgctgatatgttttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCaa                                                                                                                            <  2:430567‑M2/17‑1 (MQ=255)
              taagcgcaaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatg                                                                                                                          >  1:738960‑M2/131‑149 (MQ=255)
               aagcgcaaatttattaaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAatga                                                                                                                         <  1:350297‑M2/20‑1 (MQ=255)
                 gggcaaattcactacactgtcgctgatatggtttcatgccttcgctcctcatcttacttttttacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgcctccttctgcCGTCCCCTGCACTTCAatgatg                                                                                                                       <  2:362881‑M2/22‑1 (MQ=255)
                    caaattcaataaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATgcgc                                                                                                                    <  2:738960‑M2/25‑1 (MQ=255)
                     aaattcaataaacgggcgctgatatggtttcatgccttcgctcctcatcttacttttataaaaacaaaaaaaaggcgactcatcagtcgacttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGcc                                                                                                                   <  1:367025‑M2/26‑1 (MQ=255)
                             taaactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagacaaaacaaaggccactcctcagtcgcattaaaaatcaatttgccagcgccgccatctgcCCTCCCCCACACTTCCATGATGCCcccctctcac                                                                                                           >  2:154316‑M2/116‑146 (MQ=255)
                               aactggcgctgatatggtttcatgccttcgctcctcatcttacttttctacagataaaaaaaaggcgacgcatcagtcgccttaaaaataagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTcggg                                                                                                         <  2:679901‑M2/36‑2 (MQ=255)
                                       ctgatagggtttcatgccttcggtcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTa                                                                                                 <  2:691020‑M2/44‑1 (MQ=255)
                                       aGTTTGCCAGCGCCGCCTTCTGCGTCGCGACAATGGATTCTATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTa                                                                                                 <  2:708013/149‑1 (MQ=255)
                                         tttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                                                               >  1:33856/1‑149 (MQ=255)
                                             tggtttcatgcattcgctcctcatcttactgttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGTTCATCACTCCGTTCa                                                                                           >  2:31756‑M2/100‑149 (MQ=255)
                                               agggccgccTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATg                                                                                         <  2:432081/146‑1 (MQ=255)
                                                tttcatgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGt                                                                                        >  2:483277‑M2/97‑149 (MQ=255)
                                                    atgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                                                           <  2:726200‑M2/18‑1 (MQ=255)
                                                    atgccttcgctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAat                                                                                                                           >  1:726200‑M2/93‑110 (MQ=255)
                                                          cTGCGTCGCTACAATGGATTCGATTCCCCTCGTGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                              <  1:60459/149‑1 (MQ=255)
                                                            gctcctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCgg                                                                            <  2:110289‑M2/65‑1 (MQ=255)
                                                                cttatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttcccacgcagtgcgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                        <  1:882902‑M2/69‑1 (MQ=255)
                                                                ctcatcttacttttctacagacaaaaaaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                        >  2:365885‑M2/81‑149 (MQ=255)
                                                                cGCTACAATGGATTCGATTCCCCTCTGGCCAGATCTAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga                                                                        <  1:97207/149‑1 (MQ=255)
                                                                   tACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                                                     >  1:756751/1‑149 (MQ=255)
                                                                   atcttacttttttacagacaaaaaaaaggcgactcatcagtcgctttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                                                     <  1:374474‑M2/72‑1 (MQ=255)
                                                                            ttcgattcCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTTTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTa                                                            <  1:379307/149‑1 (MQ=255)
                                                                                     ccTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc                                                   >  1:652412/1‑149 (MQ=255)
                                                                                          aaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGt                                                                                                                                          <  2:479162‑M2/3‑1 (MQ=255)
                                                                                          aaaaggcgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGt                                                                                                                                          >  1:479162‑M2/55‑57 (MQ=255)
                                                                                                cgactcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCtt                                        <  1:900629‑M2/101‑1 (MQ=255)
                                                                                                    tcatcagtcgccttaaaaatcagtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGATCTCTACGCCAGAGACTTCAACGTATTCCAGATCGCAAACCGCTTCGcc                                    >  1:697853‑M2/45‑149 (MQ=255)
                                                                                                               tAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCGGACGCTTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTTCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                         >  1:213208/1‑149 (MQ=255)
                                                                                                               cttaaaatgcagtttgccagcgcccccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                         <  2:594543‑M2/116‑1 (MQ=255)
                                                                                                                         agtttgccagcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTACAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAAc               >  1:549649‑M2/24‑149 (MQ=255)
                                                                                                                           ggTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTg             >  1:695954/1‑149 (MQ=255)
                                                                                                                                  gcgccgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACCAAAACTGCGGCTAc      >  2:533539‑M2/15‑149 (MQ=255)
                                                                                                                                    cTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCa    <  1:716412/149‑1 (MQ=255)
                                                                                                                                      cgccttctgcCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTCCGGCAGAATCTTTACTGTAAGCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATc  >  1:659310‑M2/11‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                          
AGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATC  >  NC_000913_3_bme_pgi/3815797‑3816079

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.