New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 1900904 =20 (1.060)5 (0.270) 5/240 NT 22.5% coding (1413/4026 nt) PP_1703 assimilatory nitrate reductase/sulfite reductase
?NC_002947 1900952 = 15 (0.820)coding (1461/4026 nt) PP_1703 assimilatory nitrate reductase/sulfite reductase

TACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGACCTGCGCGGTATCGGCTATGCCGAATTGCGGACACAGCCGCGCCAGTGGCCT  >  NC_002947/1900806‑1901099
                                                                                                                                                  |                                                                                                                                                   
ggtcataaccggttgaggggttgtggaagcggcgtatttcctcatacaccgcttcggcattggggtagtcgaaggcgtcggcatagcccatggcacacGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCgtgt                                                                                                                                                  >  2:53230‑M1/99‑150 (MQ=255)
   gCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGtgtg                                                                                                                                                 >  1:152088/1‑148 (MQ=255)
   gCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGtgtg                                                                                                                                                 <  2:152088/148‑1 (MQ=255)
       accggttgaggggttgtggaagcggcgtatttcctcatacaccgcttcggcattggggtagtcgaaggcgtcggcatagcccatggcacacGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCTTGGCGTGTGCCATgg                                                                                                                                           >  2:33492‑M1/92‑150 (MQ=255)
                                  ccGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCa                                                                                                                >  1:352629/1‑150 (MQ=255)
                                           aGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAg                                                                                                       <  2:352629/150‑1 (MQ=255)
                                                   atgatCAACAGCGAGCGCAACCTGACCCTCATGCCACGCGCTGTCGAGCCCCCCGGTCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTAt                                                                                               >  1:140367/1‑150 (MQ=255)
                                                         aaCAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGaa                                                                                         <  1:302867/150‑1 (MQ=255)
                                                                               tGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAAccccg                                                                   <  2:140367/150‑2 (MQ=255)
                                                                                     cACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGc                                                                                                                                                <  2:16242/67‑1 (MQ=255)
                                                                                     cACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGc                                                                                                                                                >  1:16242/1‑67 (MQ=255)
                                                                                       cgcgCTGTCGAGCGCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCAGGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGTGGAAATACGCCGCTTCCACATCCCCTCAACCGGt                                                           <  2:139251/150‑1 (MQ=255)
                                                                                       ccatggcacacgccaagcgggcaatgatctgccagtccggcaagctttgcccggggggcGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGt                                                           >  2:33492‑M2/60‑150 (MQ=255)
                                                                                           cTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATg                                                       <  1:508720/150‑1 (MQ=255)
                                                                                              acacgccacgcgggcaatgatctgccagtccggcaagctttgcccggggggcGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGAcc                                                    >  2:53230‑M2/53‑150 (MQ=255)
                                                                                                   ccacgcgggcaatgatctgccagtccggcaagctttgcccggggggcGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGACCTgcgc                                               <  1:53230‑M2/103‑1 (MQ=255)
                                                                                                     ccccGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGAc                                                     <  1:319342/142‑1 (MQ=255)
                                                                                                     ccccGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGAc                                                     >  2:319342/1‑142 (MQ=255)
                                                                                                      cccGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGACCTGCGCGGt                                            >  1:159060/1‑150 (MQ=255)
                                                                                                         gggCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGACCTGCGCGGTATc                                         >  1:68384/1‑150 (MQ=255)
                                                                                                             aatgatctgccagtccggcaagctttgcccggggggcGTGTGCCATggg                                                                                                                                          <  1:361783‑M2/12‑1 (MQ=255)
                                                                                                             aatgatctgccagtccggcaagctttgcccggggggcGTGTGCCATggg                                                                                                                                          >  2:361783‑M2/38‑49 (MQ=255)
                                                                                                             aaaGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGACCTGCGCGGTATCGGCt                                     <  2:68384/150‑1 (MQ=255)
                                                                                                                     ccGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGACCTGCGCGGTATCGGCTATGCCGaa                             >  1:204706/1‑150 (MQ=255)
                                                                                                                                            cgtggcgtgTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGACCTGCGCGGTATCGGCTATGCCGAATTGCGGACACAGCCGCGCCAGta      <  2:159060/150‑2 (MQ=255)
                                                                                                                                                gcGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGACCTGCGCGGTATCGGCTATGCCGAATTGCGGACACAGCCGCGCCAGTGGCCt  >  2:202863‑M2/3‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TACGCCGACATCCTGCTGCCGGCCGCGCTGTGGGCCGAAGGCGAGGGTGTGATGATCAACAGCGAGCGCAACCTGACCCTGATGCCACGCGCTGTCGAGCCCCCCGGGCAAAGCTTGCCGGACTGGCAGATCATTGCCCGCGTGGCGTGTGCCATGGGCTATGCCGACGCCTTCGACTACCCCAATGCCGAAGCGGTGTATGAGGAAATACGCCGCTTCCACAACCCCTCAACCGGTTATGACCTGCGCGGTATCGGCTATGCCGAATTGCGGACACAGCCGCGCCAGTGGCCT  >  NC_002947/1900806‑1901099

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.