New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 1251009 =18 (0.780)3 (0.140) 3/234 NT 16.7% coding (268/549 nt) PP_1093 electron transport complex protein
?NC_002947 1251060 = 13 (0.600)coding (217/549 nt) PP_1093 electron transport complex protein

GTGCCTTGCGAGCAGGCGCCAGCGGCGTTGGTCAGGTAGCAGGTGGTCACTGGCCAGGCATTGCAGGCACAACAATGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGTTGGTGCAGGCACCCCATCAGTACTGCGCACAGGCCGATGGCCGCGCCTTCAGCCAGCGTTCGGGT  >  NC_002947/1250911‑1251201
                                                                                                                                                     |                                                                                                                                             
gTGCCTTGCGAGCAGGCGCCAGCGGCGTTGGTCAGGTAGCAGGTGGTCACTGGCCAGGCATTGCAGGCACAACAATGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCACCTGCAGGCAACTGGCCaa                                                                                                                                               >  2:266796/1‑150 (MQ=255)
gTGCCTTGCGAGCAGGCGCCAGCGGCGTTGGTCAGGTAGCAGGTGGTCACTGGCCAGGCATTGCAGGCACAACAATGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCaa                                                                                                                                               <  1:590087/150‑1 (MQ=255)
gTGCCTTCCGAGCAGGCGCCAGCGGCGTTGGTCAGGTAGCAGGCGGTCTGATAGTAGTGATTGCAGGCACAACAATGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCaa                                                                                                                                               <  2:324584/150‑1 (MQ=255)
actgatggggtgcctgcaccaactgttgctggcgccactgcgtaaacacatggccagctgggcctatgtactggctagcctgttgctgctggcagcctTGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGTCaa                                                                                                                                               >  1:150200‑M1/99‑150 (MQ=255)
      tGCGAGCAGGCGCCAGCGGCGTTGGTCAGGTAGCAGGTGGTCACTGGCCAGGCATTGCAGGCACAACAATGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGc                                                                                                                                         <  2:231483/150‑1 (MQ=255)
               gcgcCAGCGGCGTTGGTCAGGTAGCAGGTGGTCACTGGCCAGGCATTGCAGGCACAACAATGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCaa                                                                                                                                <  2:201128/150‑1 (MQ=255)
                                     aGCAGGTGGTCACTGGCCAGGCATTGCAGGCACAACAATGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGcc                                                                                                          <  1:21190/150‑1 (MQ=255)
                                          gtggtCACTGGCCAGGCATTGCAGGCACAACAATGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCt                                                                                                     >  2:221004/1‑150 (MQ=255)
                                           tggtCACTGGCCAGGCATTGCAGGCACAACAATGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTg                                                                                                    >  1:301382/1‑150 (MQ=255)
                                                                      aacaaTGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTg                                                                                 <  1:504992/142‑1 (MQ=255)
                                                                      aacaaTGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTg                                                                                 >  2:504992/1‑142 (MQ=255)
                                                                                         ccAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGTTGGTGCAGGCAcc                                                      <  1:521131/150‑1 (MQ=255)
                                                                                             ggCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGTTGGTGCAGGCACCCCAt                                                  <  1:148428/150‑1 (MQ=255)
                                                                                                aaTGGCCAGTGGTAGCAGACAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGTTGGTGCAGGCACCCCATCAg                                               >  2:257681/1‑150 (MQ=255)
                                                                                                 ttgaccagttgcctgcagttggccctgcgggcctggctgctaccactggccaAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGTTGGTGCAGGCACCCCATCAGt                                              >  1:150200‑M2/53‑150 (MQ=255)
                                                                                                    GCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCAtgt                                                                                             <  1:439565‑M1/100‑1 (MQ=255)
                                                                                                    GCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCAtgt                                                                                             >  2:439565‑M1/1‑100 (MQ=255)
                                                                                                      cagttgcctgcagttggccctgcgggcctggctgctaccactggccaAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGcc                                                                                                 <  2:605443‑M2/47‑1 (MQ=255)
                                                                                                      cagttgcctgcagttggccctgcgggcctggctgctaccactggccaAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGcc                                                                                                 >  1:605443‑M2/48‑94 (MQ=255)
                                                                                                         TGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGTTGGTGCAGGCACCCCATCAGTACTGCGCa                                      >  2:408192‑M1/1‑150 (MQ=255)
                                                                                                                 gcCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGTTGGTGCAGGCACCCCATCAGTACTGCGCACAGGCCGa                              <  2:301382/150‑1 (MQ=255)
                                                                                                                    aGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGTTGGTGCAGGCACCCCATCAGTACTGCGCACAGGCCGATgg                           <  2:175059/150‑1 (MQ=255)
                                                                                                                                     tGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGTTGGTGCAGGCACCCCATCAGTACTGCGCACAGGCCGATGGCCGCGCCTTcagccagc          >  2:355010/1‑150 (MQ=255)
                                                                                                                                             actggccaAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCTGTGGCGCCAGCAACAGTTGGTGCAGGCACCCCATCAGTACTGCGCACAGGCCGATGGCCGCGCCTTCAGCCAGCGTTCGGGt  <  1:221004‑M2/142‑1 (MQ=255)
                                                                                                                                               tggccaAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGt                                                                   <  2:230869‑M2/77‑1 (MQ=255)
                                                                                                                                               tggccaAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGt                                                                   >  1:230869‑M2/7‑83 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GTGCCTTGCGAGCAGGCGCCAGCGGCGTTGGTCAGGTAGCAGGTGGTCACTGGCCAGGCATTGCAGGCACAACAATGCGGGATAGTGGCCCAGGGCAATGGCCAGTGGTAGCAGCCAGGCCCGCAGGGCCAACTGCAGGCAACTGGCCAAGGCTGCCAGCAGCAACAGGCTAGCCAGTACATAGGCCCAGCTGGCCATGTGTTTACGCAGTGGCGCCAGCAACAGTTGGTGCAGGCACCCCATCAGTACTGCGCACAGGCCGATGGCCGCGCCTTCAGCCAGCGTTCGGGT  >  NC_002947/1250911‑1251201

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