New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? 1 31460 =0 (0.000)11 (0.140)
+GCTTATTCTGTGCGGCACTGTTGGCA
8/378 4.8 20.9% intergenic (‑161/+48) unknown/unknown FIGfam050825/Cytochrome c‑type heme lyase subunit nrfF, nitrite reductase complex assembly / Cytochrome c‑type heme lyase subunit nrfG, nitrite reductase complex assembly
?1 31616 = 95 (1.080)coding (942/1050 nt) unknown Cytochrome c‑type heme lyase subunit nrfF, nitrite reductase complex assembly / Cytochrome c‑type heme lyase subunit nrfG, nitrite reductase complex assembly
Rejected: Coverage evenness skew score above cutoff.

TTATTTGTACGTTTTACACTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATCCTGCGTCAAAACC  >  1/31442‑31703
                  |                                                                                                                                                                                                                                                   
gtgccgcacagaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATc               >  1:56838/19‑249 (MQ=255)
 tgccgcacagaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGAt                <  1:653252/230‑1 (MQ=255)
 tgccgcacagaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATcc              <  1:691336/232‑1 (MQ=255)
 tgccgcacagaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATcc              <  1:233080/232‑1 (MQ=255)
    cgcacagaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATCCTGCGt         <  1:398962/237‑1 (MQ=255)
     gcacagaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATCCTGCGt         >  1:205554/14‑250 (MQ=255)
     gcacagaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATCCTGCGt         >  1:324817/14‑250 (MQ=255)
          gaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACcattttca                                                                                                                                                             <  1:559786/89‑1 (MQ=255)
          gaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCAt                                                                                                                         >  1:508765/9‑133 (MQ=255)
          gaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGaa                                                                                                          <  1:54309/140‑1 (MQ=255)
          gaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAg                                                                                       <  1:40642/159‑1 (MQ=255)
          gaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGaat                                                                            <  1:16261/170‑1 (MQ=255)
          gaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATCCTGCGt         <  1:200001/237‑1 (MQ=255)
          gaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATCCTGCGTCaaaa    >  1:418523/9‑250 (MQ=255)
          gaataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATCCTGCGTCaaaa    >  1:127057/9‑250 (MQ=255)
           aataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATCCTGCGTCAAAAcc  <  1:332372/244‑1 (MQ=255)
           aataagccTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATCCTGCGTCAAAAc   <  1:406124/243‑1 (MQ=255)
                  |                                                                                                                                                                                                                                                   
TTATTTGTACGTTTTACACTGCCGCCCCGGCGCCGATAATGAGGGTCTCATGGCTGACGGTTCGACTCATTGCATATTCCCCTTCTGCATAAAATCCACCATTTTCATTTTCTGAATGATCGCACGACGATCCACCGTGGCATTGCCCGAATCTAGAATTTGCTGCCAATAATCAAGTGCGGTGGAATAATCCCGCGTTTTTAAGGCTTCCGCCGCCAAAAGGGACAGAACCGATACCTCTTTGCGATCCTGCGTCAAAACC  >  1/31442‑31703

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