Predicted mutation
evidence seq id position mutation freq annotation gene description
JC AE015451 5,681,415 +CGGG 100% intergenic (‑133/+71) PP_4986 ← / ← PP_4987 Channel protein, hemolysin III family/putative Chemotaxis protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 = 56814150 (0.000)18 (0.720)
+CGGG
16/310 NT 100% intergenic (‑133/+71) PP_4986/PP_4987 Channel protein, hemolysin III family/putative Chemotaxis protein
?AE015451 5681416 = 0 (0.000)intergenic (‑134/+70) PP_4986/PP_4987 Channel protein, hemolysin III family/putative Chemotaxis protein

GCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGCGCC  >  AE015451/5681197‑5681606
                                                                                                                                                                                                                          |                                                                                                                                                                                               
gCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcag                                                                                                                                                                                 <  1:164921‑M1/235‑17 (MQ=255)
                              gcGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCGGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgc                                                                                                                                                           <  2:219075‑M1/227‑39 (MQ=255)
                                                     ggTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaag                                                                                                                                                                <  1:173759‑M1/199‑34 (MQ=255)
                                                     ggTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaag                                                                                                                                                                >  2:173755‑M1/1‑166 (MQ=255)
                                                                                                                         gttgtGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctc                                                                                                                                                         <  2:415535‑M1/138‑41 (MQ=255)
                                                                                                                               aTGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgc                                                                                                                               >  2:238303‑M1/1‑92 (MQ=255)
                                                                                                                                 gCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtca                                                                                        >  1:193444‑M1/1‑90 (MQ=255)
                                                                                                                                  cTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgt                                                                                          >  1:195148‑M1/1‑89 (MQ=255)
                                                                                                                                     aTGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccc                                                                                                                                    <  2:141895‑M1/147‑62 (MQ=255)
                                                                                                                                       ggCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggg                                                                                                                                                                                             >  2:201253‑M1/1‑84 (MQ=255)
                                                                                                                                       ggCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggg                                                                                                                                                                                             <  1:201257‑M1/88‑5 (MQ=255)
                                                                                                                                                           cGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgc                                        <  1:432289‑M1/217‑154 (MQ=255)
                                                                                                                                                                   tgtCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtca                                                                                        >  2:272805‑M1/1‑56 (MQ=255)
                                                                                                                                                                                      aGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgctgccagctccaacctgcccag  <  2:281223‑M1/228‑192 (MQ=255)
                                                                                                                                                                                      aGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctaccgcgacacaacgccgctcccactccagaccaccaggccctgcgcctcacctcctcacccagctcacgttgcag                                                                                                  >  2:2853‑M1/1‑37 (MQ=255)
                                                                                                                                                                                                        ttGCCACTGCATGAACCGCcggggctgctgcgcagccctatcg                                                                                                                                                                         >  2:505305‑M1/1‑19 (MQ=255)
                                                                                                                                                                                                                  aTGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgc                                                                                            >  1:499742‑M1/1‑9 (MQ=255)
                                                                                                                                                                                                                  aTGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgc                                                                                            <  2:499738‑M1/110‑102 (MQ=255)
                                                                                                                                                                                                                  aTGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggt                                                                                                       >  2:252739‑M1/1‑9 (MQ=255)
                                                                                                                                                                                                                  aTGAACCGCcggggctgctgcgcagccctatcgcgacacaaggc                                                                                                                                                              >  2:195372‑M1/1‑9 (MQ=255)
                                                                                                                                                                                                                     acccgccggggctgcggcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagccttt                                                                           <  2:20728‑M1/122‑119 (MQ=255)
                                                                                                                                                                                                                        cgccggggctgctgcgcagccctatcgcgccacaaggccgctcccact                                                                                                                                                    <  2:74824‑M1/48‑46 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                               
GCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGCGCC  >  AE015451/5681197‑5681606

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