New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 330399622 (1.080)4 (0.210) 4/276 1.7 17.6% coding (549/1008 nt) yhbW putative luciferase‑like monooxygenase
?NC_000913 = 3304188 17 (0.900)coding (741/1008 nt) yhbW putative luciferase‑like monooxygenase
Rejected: Frequency below cutoff threshold.

AATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGCCCGCGATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCGCAGCAACTTCAAACCGTCAGCACGGCTGGAAAAACCATACGCGATGGTGTGCATCAATATTATCGCCGCCGACAGCAAC  >  NC_000913/3303847‑3304143
                                                                                                                                                     |                                                                                                                                                   
aaTTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGCCCGCGATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTg                                                                                                                                                    >  2:309091/1‑151 (MQ=255)
    tCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGCCCGCGATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCt                                                                                                                                                <  1:295919/151‑1 (MQ=255)
               gtgGCGGAGCTGGTGGACTGGTTTGACGCCCGCGATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTcacaaaggcttggtg                                                                                                                                      >  1:210922‑M1/1‑135 (MQ=255)
                        cTGGTGGACTGGTTTGACGCCCGCGATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTcacaaaggcttgctgcattgagg                                                                                                                              <  1:92364‑M1/149‑24 (MQ=255)
                           gtggACTGGTTTGACGCCCGCGATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCacac                                                                                                                          <  1:274569/150‑1 (MQ=255)
                                           cccGCGATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTcacaaaggcttgctgcattgaggtaaacagaaattcagcgtcg                                                                                                          >  1:269635‑M1/1‑107 (MQ=255)
                                           cccGCGATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCt                                                                                                         <  1:90073/151‑1 (MQ=255)
                                            ccGCGATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTg                                                                                                        <  1:191673/151‑1 (MQ=255)
                                            ccGCGATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTg                                                                                                        <  2:138349/151‑1 (MQ=255)
                                               cgATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGtt                                                                                                      >  2:8514/1‑150 (MQ=255)
                                                 aTCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCa                                                                                                   >  1:136954/1‑151 (MQ=255)
                                                 aTCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCa                                                                                                   >  1:230226/1‑151 (MQ=255)
                                                        tCCGCATGTGCGCCCGGTACCAGGCTATGGCGGGAGAATCCACGTGTGGTTGGTAGGCGCCCGCCTTTACAGCGCGCAACTCGCGGCGCAGCTTGGTCTGCCGCTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGgcg                                                                                               >  2:203280/1‑148 (MQ=255)
                                                           gCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCAt                                                                                         >  2:281876/1‑151 (MQ=255)
                                                                        gTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAATTGGCGGCGCAGCTTcacaaaggcttgctgcattgaggtaaacagaaattcagcgtcgcggttgctgtcggcggcgataatattgatg                                                                            >  1:105415‑M1/1‑78 (MQ=255)
                                                                        gTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCgcagca                                                                            >  2:120357/1‑151 (MQ=255)
                                                                        gTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCgcagca                                                                            >  2:148523/1‑151 (MQ=255)
                                                                        gTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCgcagca                                                                            >  2:165004/1‑151 (MQ=255)
                                                                        gTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCgcagca                                                                            >  2:201799/1‑151 (MQ=255)
                                                                        gTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCgcagca                                                                            >  1:228676/1‑151 (MQ=255)
                                                                                                               ggCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCGCAGCAACTTCAAACCGTCAGCACGGCTGGAAAAACCATACGCGa                                     >  2:242614/1‑151 (MQ=255)
                                                                                                                  tCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCGCAGCAACTTCAAACCGTCAGCACGGCTGGAAAAACCATACGCGATgg                                  <  2:48353/151‑1 (MQ=255)
                                                                                                                      gCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCGCAGCAACTTCAAACCGTCAGCACGGCTGGAAAAACCATACGCGATGGtgtg                              >  1:144063/1‑151 (MQ=255)
                                                                                                                       ccTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATTGCAGCAACTTCAAACCGTCAGCACGGCTGGAAAAACCATACGCGATGGtgtg                              >  2:307360/1‑150 (MQ=255)
                                                                                                                          ttACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCGCAGCAACTTCAAACCGTCAGCACGGCTGGAAAAACCATACGCGATGGTGTGCATc                          <  1:120402/151‑1 (MQ=255)
                                                                                                                                                   cTTGGTCTGCCGGTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCGCAGCAACTTCAAACCGTCAGCACGGCTGGAAAAACCATACGCGATGGTGTGCATCAATATTATCGCCGCCGACAGCAAc  <  2:200773/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AATTTCCCCCGCGATGTGGCGGAGCTGGTGGACTGGTTTGACGCCCGCGATCCCAATCCGCATGTGCGCCCGGTACCAGGCTATGGCGAGAAAATCCCCGTGTGGTTGTTAGGCTCCAGCCTTTACAGCGCGCAACTGGCGGCGCAGCTTGGTCTGCCGTTTGCGTTTGCCTCACACTTCGCGCCGGATATGCTGTTCCAGGCGCTGCATCTTTATCGCAGCAACTTCAAACCGTCAGCACGGCTGGAAAAACCATACGCGATGGTGTGCATCAATATTATCGCCGCCGACAGCAAC  >  NC_000913/3303847‑3304143

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