New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 2489400 =10 (0.580)4 (0.240) 4/268 1.6 28.9% coding (354/1578 nt) PP_2184 formate dehydrogenase subunit beta
?NC_002947 2489439 = 10 (0.600)coding (393/1578 nt) PP_2184 formate dehydrogenase subunit beta

ACCCGGTGGCCACTCGCTGGCACTGGGGCCGGTCGAAGAAATACCCTATCTCAAGACCCAGCAACGCTTGCTGTTCGCTCGCGCTGGCATTACCCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGAT  >  NC_002947/2489307‑2489568
                                                                                                                                    |                                                                                                                                 
gccgcaccaccacggccacgcaggccggagtcgagtactgcggcaaccacttcggcaccgtccagcgccacagcggcctccaggcctgcgaaaCCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAg                                                                                                                    >  2:189265‑M1/94‑148 (MQ=255)
gccgcaccaccacggccacgcaggccggagtcgagtactgcggcaaccacttcggcaccgtccagcgccacagcggcctccaggcctgcgaaaCCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAg                                                                                                                    <  1:189265‑M1/55‑1 (MQ=255)
   cGGTGGCCACTCGCTGGCACTGGGGCCGGTCGAAGAAATACCCTATCTCAAGACCCAGCAACGCTTGCTTTTCGCTCGCGCTGGCATTACCCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGc                                                                                                               >  1:115251/1‑150 (MQ=255)
    ggTGGCCGCTCGCTGGCACTGGGGTGTGGCGAAGAAATACCCTATCTCAAGACCCGGCAACGCTTGCTGTTCGCTCGCGCTGGCATTACCCGCCCGTTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCt                                                                                                              <  2:879/150‑1 (MQ=255)
                  ggCACTGGGGCCGGTCGAAGAAATACCCTATCTCAAGACCCAGCAACGCTTGCTGTTCGCTCGCGCTGGCATTACCCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACgg                                                                                                <  2:31731/150‑1 (MQ=255)
                                     gaaATACCCTATCTCAAGACCCAGCAACGCTTGCTTTTCGCTCGCGCTGGCATTACCCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCa                                                                             <  2:115251/150‑1 (MQ=255)
                                               aTCTCAAGAACCAGCAACGCTTGCTGTTCGCTCGCGCTGGCATTACCCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACt                                                                   <  1:123680/150‑1 (MQ=255)
                                                tctcAAGACCCAGCAACGCTTGCTGTTCGCTCGCGCTGGCATTACCCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGTCGCTGTGGCGCTGGACGGTGGCGAAGTGGTTGCCGCAGTACTTGACTc                                                                  >  2:398526/1‑150 (MQ=255)
                                                          cAGCAACGCTTGCTGTTCGCTCGCGCTGGCATTACCCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGt                                                        >  2:177499/1‑150 (MQ=255)
                                                                             ctccaggcctgcgaaaccgccgttggcgcgatagtcgtccagcgacaatgggcggTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTgcggc                                       <  1:189265‑M2/93‑1 (MQ=255)
                                                                             ctccaggcctgcgaaaccgccgttggcgcgatagtcgtccagcgacaatgggcggTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTgcggc                                       >  2:189265‑M2/56‑148 (MQ=255)
                                                                                             CCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACt                                                                   >  2:149726‑M1/1‑104 (MQ=255)
                                                                                             CCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACt                                                                   <  1:149726‑M1/104‑1 (MQ=255)
                                                                                              cgccgttggcgcgatagtcgtccagcgacaatgggcggTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAg                    >  1:232261‑M2/39‑150 (MQ=255)
                                                                                                  gttggcgcgatagtcgtccagcgacaatgggcggTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGc                <  2:232261‑M2/116‑1 (MQ=255)
                                                                                                       tCGCTGGACGACTATCGCTCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTAcc           <  1:184761/150‑1 (MQ=255)
                                                                                                       tCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTAcc           <  1:177499/150‑1 (MQ=255)
                                                                                                                gacTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGAt  >  1:168570/1‑150 (MQ=255)
                                                                                                                                    |                                                                                                                                 
ACCCGGTGGCCACTCGCTGGCACTGGGGCCGGTCGAAGAAATACCCTATCTCAAGACCCAGCAACGCTTGCTGTTCGCTCGCGCTGGCATTACCCGCCCATTGTCGCTGGACGACTATCGCGCCAACGGCGGTTTCGCAGGCCTGGAGGCCGCTGTGGCGCTGGACGGTGCCGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGAT  >  NC_002947/2489307‑2489568

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