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

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 56814152 (0.060)20 (0.580)
+CGGG
18/282 1.1 91.1% intergenic (‑133/+71) PP_4986/PP_4987 hemolysin III family channel protein/putative Chemotaxis protein
?NC_002947 5681416 = 2 (0.060)intergenic (‑134/+70) PP_4986/PP_4987 hemolysin III family channel protein/putative Chemotaxis protein

TAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGC‑‑‑‑G‑‑‑CTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGC  >  NC_002947/5681275‑5681555
                                                                                                                                                  |                                                                                                                                              
gtaatacatgcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggg‑‑‑c                                                                                                                                            <  2:87296‑M2/2‑1 (MQ=255)
     acatgcaaacacagacctcctgagctgggctgggttcctgttgtgatgctgatggcttcgttctttgtcggagcggccttgtgtcccgaaaggggcccgaagccgccccaggctttcatttccactgcatgaaccgccgcgg‑‑‑cggctt                                                                                                                                       >  7:57162‑M2/142‑147 (MQ=255)
         gcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttgatttgcccagcatgaaccgccgggg‑‑‑cggcggcggc                                                                                                                                   >  3:33506‑M2/138‑146 (MQ=255)
         gcaaacacagacctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggG‑‑‑CTGCTGCGCa                                                                                                                                   <  4:44302‑M2/11‑1 (MQ=255)
               acagacctcctgaggtggtctgggttcctgttgtgatgctgatgggttcgttctttgtgggagcggccttgtttcgcgaaaggggcgcgaagaggccccaggctttcattgccactgcatgaaccgccgggG‑‑‑CTGCTGCGCAGCCCTa                                                                                                                             <  2:21348‑M2/17‑1 (MQ=255)
                   agctgctgtggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtggcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggG‑‑‑CTGCTGCGCAGCCCTATcgc                                                                                                                         <  6:16144‑M2/21‑1 (MQ=255)
                   acctcctgaggtgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccacgctttcattgccactgcatgaaccgccgggG‑‑‑CTGCTGCGCAGCCCAATagc                                                                                                                         >  6:83273‑M2/128‑148 (MQ=255)
                              tgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggG‑‑‑CTGCTGCGCAGCCCTATCGCGACACAAGGc                                                                                                               <  8:89448‑M2/31‑1 (MQ=255)
                              tgggctgggttcctgttgtgatgctgatggcttcgttctttgtgggagcggccttgtgtcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggG‑‑‑CTGCTGCGCAGCCCTATCGCGACACAAGGc                                                                                                               <  8:37515‑M2/31‑1 (MQ=255)
                                                                      tgtggtgCGGCATTGGTTCGTGAATGGTGCGAGAAGCGGCCCCTGTCTTTCATTGCCACTGCATGAACCGAC‑‑‑‑GGTTCTGCTTCCCAGCCCTATCTCGAAACAAGGCCGCTCCCACTCCAGACCACCATGCCCTGCGCCTCACATc                                                                        <  4:128975/139‑1 (MQ=255)
                                                                                       gtcgcgaaaggggcgcgacgcggcccctggctttctttgccagtgattgaaccgccgggG‑‑‑CTGCGGGGCCGCCCTATCGCGAAACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGc                                                     <  6:74715‑M2/89‑1 (MQ=255)
                                                                                        tcgcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggG‑‑‑CTGCTGCGCAGCCCTATAGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAACGTCAGGtt                                                       >  8:7031‑M2/59‑145 (MQ=255)
                                                                                          gcgaaaggggcgcgaagcggccccaggctttcattgccactgcatgaaccgccgggG‑‑‑CTGCTGCGCAGCCCTATCGCGCCACAAGGCCGCTCCCACTCCAGACCACCCCGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGa                                                  >  7:134364‑M2/57‑148 (MQ=255)
                                                                                            aaaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGCGG‑‑‑CAGCTGCGCACCCCTATACCCAAACAAGGCCCCTCCCCCTCCAAAACACCACCCCCTGAGCCTCAAAAACAAACCAAGGTCATATTGCacaaca                                               >  2:141667/1‑143 (MQ=255)
                                                                                                          ggggcccctggctttaatgcccactgcatcaaccgcggggG‑‑‑CGGCGGCGCAGCCATAGCGCCACAAATGCCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCa                                  <  5:83273‑M2/108‑1 (MQ=255)
                                                                                                            ggccccaggatttcattgccactgaatgaaccgccgggG‑‑‑CGGCTGCGCACCCCTAGCGCGACACATGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGc                                <  8:186690‑M2/110‑1 (MQ=255)
                                                                                                               cccaggctttcattgccactgcatgaaccgccgggG‑‑‑CTGCTGCGCAGCCCGATCGCGCCACAAGGCCGCGCCCACTCCCGCCCACCAGGCCCGGCGCCTCACATCCTAACCAAGGTCACGTTGCCGCTCCGCGGCACGCCCCCGCCtt                             >  5:147062‑M2/36‑148 (MQ=255)
                                                                                                                                actgcatgaaccgccgggG‑‑‑CTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGCCCACCAGGCCCTGCGCCGCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCCCGCCCCAGCCGTTCCAATGCGCCCAGGTc            >  3:76789‑M2/19‑148 (MQ=255)
                                                                                                                                actgcatgaaccgccgggG‑‑‑CTGCTGCGCAGCCCTATAGCGACACAAGGCCGCTCCCACTACAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTc            >  8:45778‑M2/19‑148 (MQ=255)
                                                                                                                                     atgaaccgccgggG‑‑‑CTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt                            >  8:106723‑M2/14‑127 (MQ=255)
                                                                                                                                      tgaaccgccgggG‑‑‑CTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTc                           <  7:106723‑M2/115‑1 (MQ=255)
                                                                                                                                             ccgggG‑‑‑CTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCcg   >  6:110846‑M2/6‑144 (MQ=255)
                                                                                                                                              cgggG‑‑‑CTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCcgc  <  5:110846‑M2/140‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGC‑‑‑‑G‑‑‑CTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGC  >  NC_002947/5681275‑5681555

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