Predicted mutation
evidence seq id position mutation freq annotation gene description
JC NC_002947 5,555,608 +GCC 100% intergenic (‑84/+87) PP_4887 ← / ← PP_4888 hypothetical protein/methyl‑accepting chemotaxis transducer

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 55556081 (0.020)28 (0.570)
+GCC
25/252 NT 96.6% intergenic (‑84/+87) PP_4887/PP_4888 hypothetical protein/methyl‑accepting chemotaxis transducer
?NC_002947 5555609 = 1 (0.020)intergenic (‑85/+86) PP_4887/PP_4888 hypothetical protein/methyl‑accepting chemotaxis transducer

CGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCC  >  NC_002947/5555461‑5555746
                                                                                                                                                   |                                                                                                                                             
cGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAgcg                                                                                                                                               <  7:189682/148‑1 (MQ=255)
cGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAgcg                                                                                                                                               <  6:268510/148‑1 (MQ=255)
cGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAgcg                                                                                                                                               <  5:170927/148‑1 (MQ=255)
  cTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGTCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgcc                                                                                                                                            <  4:51153‑M1/149‑4 (MQ=255)
       aCTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccgg                                                                                                                                        <  5:180549‑M1/148‑8 (MQ=255)
           tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatt                                                                                                                                   >  7:179442‑M1/1‑137 (MQ=255)
             cAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCGCATAGCGgccccgggcaatta                                                                                                                                 >  6:97319‑M1/1‑135 (MQ=255)
              aGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttcc                                                                                                                                >  5:272892‑M1/1‑134 (MQ=255)
                        ggTCGGCCTTGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgc                                                                                                                      <  6:10252‑M1/149‑26 (MQ=255)
                          tCGGCCTGGGCAGCTATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaa                                                                                                                    <  6:187529‑M1/149‑28 (MQ=255)
                                             gTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctat                                                                                                  <  6:272892‑M1/148‑46 (MQ=255)
                                              tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgc                                                                                                <  8:179442‑M1/149‑48 (MQ=255)
                                                gggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgca                                                                                               <  1:203772‑M1/148‑49 (MQ=255)
                                                                   cTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaagg                                                                            <  5:97319‑M1/148‑68 (MQ=255)
                                                                    ttttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagc                                                                                             >  1:284789‑M1/1‑80 (MQ=255)
                                                                    ttttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagc                                                                                             <  2:284789‑M1/130‑51 (MQ=255)
                                                                              tgccgccGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaagga                                                                 <  7:293195‑M1/145‑79 (MQ=255)
                                                                                  gCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggg                                                                                                          >  1:156677‑M1/1‑66 (MQ=255)
                                                                                        aGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgc                                                                                                >  1:64187‑M1/1‑60 (MQ=255)
                                                                                        aGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgc                                                                                                <  2:64187‑M1/107‑48 (MQ=255)
                                                                                                aaCACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccac                                                                       <  5:107687‑M1/124‑73 (MQ=255)
                                                                                                              gtgGGAGCGGCCTTGTGTACCGAAAAGGCCCCAAAGCGGCCCCCGGAAATTACCCGGGTTTGCAACAAACCAGGGGCGGCTATGCAACCCTTTAGCGACCAAAAGCCCACAAACGACCGCGCCCCGATAAGGCCCCGAACTGGccaccg                                >  4:245940/1‑145 (MQ=255)
                                                                                                                     cGGCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgc                         <  6:233467‑M1/149‑119 (MQ=255)
                                                                                                                       gCCTTGTGTCGCGAAAGGGCCGCATAGCGgccccggaaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgc                       <  3:51504‑M1/149‑121 (MQ=255)
                                                                                                                        ccTTGTGTCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaac                                         >  5:41130‑M1/1‑28 (MQ=255)
                                                                                                                               tCGCGAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagc               >  5:33531‑M1/1‑21 (MQ=255)
                                                                                                                                 gcgAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccccaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagatgagccg             >  5:333635‑M1/1‑19 (MQ=255)
                                                                                                                                  cgAAAGGGCCGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgc            >  4:256896‑M1/1‑18 (MQ=255)
                                                                                                                                           cGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatc   <  6:33531‑M1/149‑141 (MQ=255)
                                                                                                                                           cGCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatc   >  8:135056‑M1/1‑9 (MQ=255)
                                                                                                                                            gCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatcc  <  3:175515‑M1/149‑142 (MQ=255)
                                                                                                                                            gCATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatcc  <  1:221879‑M1/149‑142 (MQ=255)
                                                                                                                                             cATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcagccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatc   >  5:158241‑M1/1‑7 (MQ=255)
                                                                                                                                             cATAGCGgccccggcaatttccgcgtgtatgcaacaaacctggggctgctatgcaaccctttcgcgacacaaggcccacaaaggaccgcgcaggcttcaggcccggaactggcccaggctggcgcgcagctgagccgccagatcatc   >  6:179798‑M1/1‑7 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCG‑‑‑CCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCC  >  NC_002947/5555461‑5555746

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.