Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 1,421,597 IS5 (+) +4 bp 6.4% 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 =131 (1.120)4 (0.030) 4/278 NT 2.9% 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 = 1421600131 (1.120)14 (0.120) 13/278 NT 9.4% coding (482/858 nt) ydaU Rac prophage; conserved protein
?NC_000913 = 2290113 NA (NA)noncoding (1/1195 nt) IS5 repeat region

AAGCATCTGTTCAGGCCAGAAGAAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1421471‑1421600
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTT  <  NC_000913/2290113‑2289972
                                                                                                                                                                                                                                                                                
AAGCATCTGTTCAGGCCAGAAGAAGCATAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCC                                                                                                                            <  1:881119/150‑1
                CAGAAGAAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGT                                                                                                            >  1:1220866/1‑150
                       AGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCAT                                                                                                     <  1:1041899/150‑1
                                         GTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGG                                                                                   >  1:805465/1‑150
                                           TCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTAGAAGTGGTGTTTGAACATGATGTCAACACAAAGCCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTTGAGCCATAGAACAGGct                                                                                 >  1:1188994/1‑148
                                                CAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCAGCTGGAGCCATAGAACAGGGTTCATC                                                                            <  1:375859/150‑1
                                                            GAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGATTCATCAA                                                                >  1:67799/1‑150
                                                                                 CAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAA                                           >  2:1808125/1‑150
                                                                                 CAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAA                                           >  2:731053/1‑150
                                                                                       TTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATCTGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGC                                     <  2:1220865/150‑1
                                                                                        TGAAGTGGTGTCTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCA                                    <  1:302765/150‑1
                                                                                           AGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTA                                 <  1:104082/150‑1
                                                                                                                 AACACAAAGGCAAATAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTACATCATGAGACATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGAAGTCAGACCAGAAAAta           >  1:81105/1‑148
                                                                                                                          GCAACTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTT  <  1:731054/150‑1
                                                                                                                                                                                                                                                                                
AAGCATCTGTTCAGGCCAGAAGAAGCAGAAAGGAAGCAGATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1421471‑1421600
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTT  <  NC_000913/2290113‑2289972

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