Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 1,421,597 IS5 (+) +4 bp 5.9% coding (479‑482/858 nt) ydaU → Rac prophage; conserved protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1421597 =223 (0.890)12 (0.050) 6/278 NT 5.0% coding (479/858 nt) ydaU Rac prophage; conserved protein
?NC_000913 2101749 = NA (NA)noncoding (1195/1195 nt) IS5 repeat region
Rejected: Frequency below/above cutoff threshold.
* ? NC_000913 = 1421600223 (0.890)16 (0.060) 9/278 NT 6.5% coding (482/858 nt) ydaU Rac prophage; conserved protein
?NC_000913 = 2290113 NA (NA)noncoding (1/1195 nt) IS5 repeat region

AAAAAATCAGCCGCAGGAAAAGCATCTGTTCAGGCCAGAAGAAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1421452‑1421600
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGA  <  NC_000913/2290113‑2289981
                                                                                                                                                                                                                                                                                          
AAAAAATCAGCCGCAGGAAAACCATCTGTTCAGTCCAGAAGAAGCAGAAAGGAAGCAGATCTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGAAGTTGAAGAGGTGTTTGAACAAGATGTCAACACAAAGGCAACTAAG                                                                                                                                      >  2:3239872/1‑150
      TCAGCCGCAGGAAAAGCATCTGTTCAGGCCAGAAGAAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGTAAGGT                                                                                                                                >  1:2563792/1‑150
      TCAGCCGCAGGAAAAGCATCTGTTCAGGCCAGAAGAAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGT                                                                                                                                >  1:2150608/1‑150
      TCAGCCGCAGGAAAAGCATCTGTTCAGGCCAGAAGAAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGT                                                                                                                                >  1:2150636/1‑150
                                 GCCAGAAGAAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCAT                                                                                                     >  2:1537568/1‑150
                                                             TTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGG                                                                         <  2:1042735/150‑1
                                                                   CAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACACAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATC                                                                   <  2:3610786/150‑1
                                                                   CAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATC                                                                   <  2:3610647/150‑1
                                                                                        CAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTC                                              <  2:4223343/150‑1
                                                                                          GGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGC                                            <  1:3948596/150‑1
                                                                                          GGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGC                                            <  1:1794707/150‑1
                                                                                          GGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGC                                            <  1:2344115/150‑1
                                                                                                    CAGATGTTGAAGTGGTGTTTGAACTTGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACTAGTCACTGATATGAGTTCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAA                                  <  1:3239872/150‑1
                                                                                                                                  aCAACACAAACCCAACTAACGAACGAGCGAACAAGTCCCAGATAAGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAA    <  2:2563792/149‑1
                                                                                                                                  TCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAA    <  2:2150608/150‑1
                                                                                                                                  TCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAA    <  2:2150636/150‑1
                                                                                                                                    AACACAAAGTCAACTAAGGAAGGTGCGATCAAGTCCCTGATATGAGATCATGTTTGTAATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCTTCAGACCAGAAAAGA  <  2:2247088/150‑1
                                                                                                                                                                                                                                                                                          
AAAAAATCAGCCGCAGGAAAAGCATCTGTTCAGGCCAGAAGAAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1421452‑1421600
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGA  <  NC_000913/2290113‑2289981

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ 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.