Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 1,633,913 IS5 (–) +4 bp coding (292‑295/591 nt) pinQ → putative recombinase PinQ

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1633913 =0 (0.000)17 (0.890) 14/276 0.2 100% coding (292/591 nt) pinQ putative recombinase PinQ
?NC_000913 = 2290113 NA (NA)noncoding (1/1195 nt) IS5 repeat region
* ? NC_000913 = 1633916NA (NA)20 (1.030) 16/278 0.1 NA coding (295/591 nt) pinQ putative recombinase PinQ
?NC_000913 2101749 = NA (NA)noncoding (1195/1195 nt) IS5 repeat region

TCGAAAGTGCAGGTTTTAAAATCAAACCTCAGCAAATAATCGAAGAACACATTAGCGGCTCAGCAGCAACCAGTGAGCGTCCTGGTTTTAACCGGTTGCTTGCTCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAAT  >  NC_000913/1633692‑1634080
                                                                                                                                                                                                                                |                                                                                                                                                                    
gccattaccagccttggagtaaaacggctcggtgacttccaccatgttttgaaatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctggcggcgcttaatgctgaattcactgacggcgaaggtaagttcatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGCATTGGGGGGCATTGACCTGACCAGt                                                                                                                                            <  1:151605‑M1/30‑1 (MQ=255)
   nttaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctgggctgatggcgcttactgctgaattcactgtcggcgaaggtaagttgctgactcatgatgaaccctgttctatggctccagatgacaaacatgatcccatatcagggacttggtcgcaccttccTTCGCAGTGGGGGTCATtgac                                                                                                                                                     <  1:300401‑M1/21‑1 (MQ=255)
   attaccagccttggggtaaaacggctcgatgtcttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGCATTGGGGGGCATTGACCTGACCAGt                                                                                                                                            >  1:288211‑M1/219‑248 (MQ=255)
        cagccgtggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctctgttctggtctgacggcgcgtactgctgaagtcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttcgatggctccagatgacaaacatgatcgcagatcatggacttgttcgcaccttccGTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAgg                                                                                                                                    >  1:219373‑M1/214‑251 (MQ=255)
                                      ccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGCATTGGGGGGCATTGACCTGACCAg                                                                                                                                             <  1:324245‑M1/29‑1 (MQ=255)
                                                        ncagaatctgatccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACg                                                                                       <  1:146832‑M1/83‑1 (MQ=255)
                                                                                        nnnnnnnnnntcnnnnngcgnttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGCATTGGGGGGCATtgac                                                                                                                                                     >  2:300401‑M1/134‑154 (MQ=255)
                                                                                                                ctgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGaa                                                                                         <  1:410196‑M1/81‑1 (MQ=255)
                                                                                                                          cactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgccaaacatgatctcatatcagggaattgttcgcaacttccTTAGCATTGGGGGGCATTGACCTGCCCAGTCCAACCGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGACTTTGAACgag                                                                                     <  1:434994‑M1/85‑1 (MQ=255)
                                                                                                                                     aaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcacctgccTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACctgtgn       >  1:58986‑M1/89‑250 (MQ=255)
                                                                                                                                       ggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGCATTgggggg                                                                                                                                                            >  1:442481‑M1/87‑100 (MQ=255)
                                                                                                                                          aagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAAt  <  1:172142‑M1/168‑1 (MQ=255)
                                                                                                                                                      catgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGCATTGGGGTGCATTGACCTGACCAGGCCAACAGGAAAAATGATGATTCAATTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTc                                                                >  1:323038‑M1/72‑177 (MQ=255)
                                                                                                                                                                                                    tatcagggacttgttcgcaccttccTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGTATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCg                  >  1:225628‑M1/26‑177 (MQ=255)
                                                                                                                                                                                                     atcagggacttgttcgcaccttccTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCg                                                          >  1:201275‑M1/25‑136 (MQ=255)
                                                                                                                                                                                                         tggactttttcgcgccttccTTAGCATTGGGGGTCATTGACCTGTCCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTca                                                                                                           >  1:205397‑M1/21‑83 (MQ=255)
                                                                                                                                                                                                                   cgcaccttccTTAGCAGTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATGTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCggg                                                        >  1:88628‑M1/11‑124 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                    
TCGAAAGTGCAGGTTTTAAAATCAAACCTCAGCAAATAATCGAAGAACACATTAGCGGCTCAGCAGCAACCAGTGAGCGTCCTGGTTTTAACCGGTTGCTTGCTCGCCTGAAATGTGGTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCTGTGTTAAAT  >  NC_000913/1633692‑1634080

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ 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.