New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 444080139 (0.990)4 (0.110) 4/246 4.6 10.7% coding (415/1344 nt) ytfL UPF0053 family inner membrane protein
?NC_000913 = 4440841 30 (0.810)coding (375/1344 nt) ytfL UPF0053 family inner membrane protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAA  >  NC_000913/4440694‑4440967
                                                                                                                                                   |                                                                                                                              
gTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATAc                                                                                                                               >  5:178778/1‑149 (MQ=255)
 tCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATAc                                                                                                                               <  4:27139/148‑1 (MQ=255)
              ccATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCggg                                                                                                                 <  6:5832/149‑1 (MQ=255)
              ccATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCggg                                                                                                                 <  4:104111/149‑1 (MQ=255)
                          tGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGc                                                                                                      <  7:165083/148‑1 (MQ=255)
                           gAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCa                                                                                                     >  4:153303/1‑148 (MQ=255)
                                tACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCgggggt                                                                                                              <  1:206061/134‑3 (MQ=255)
                                tACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCgggggt                                                                                                              >  2:206061/1‑132 (MQ=255)
                                tACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAgg                                                                                               <  3:153303/149‑1 (MQ=255)
                                  cGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGa                                                                                              >  2:61540/1‑148 (MQ=255)
                                    gAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATa                                                                                            >  4:70914/1‑148 (MQ=255)
                                     aaGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAg                                                                                                <  1:193968/143‑1 (MQ=255)
                                     aaGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAg                                                                                                >  2:193968/1‑143 (MQ=255)
                                        aTTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACAt                                                                                       <  6:178778/149‑1 (MQ=255)
                                                 gCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCa                                                                               <  5:4302/148‑1 (MQ=255)
                                                 gCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCa                                                                               >  6:4302/1‑148 (MQ=255)
                                                          ttGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGtt                                                                                                              <  6:58349/108‑1 (MQ=255)
                                                          ttGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGtt                                                                                                              >  5:58349/1‑108 (MQ=255)
                                                                      acCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATaaa                                                         <  3:70914/149‑1 (MQ=255)
                                                                              ggTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGtt                                                 <  1:61540/149‑1 (MQ=255)
                                                                                         aCAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCagagag                                       >  6:49880/1‑148 (MQ=255)
                                                                                         aCAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCagagag                                       <  5:49880/148‑1 (MQ=255)
                                                                                             gCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTg                                  >  6:26902/1‑149 (MQ=255)
                                                                                                      cATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACata                          >  6:39655/1‑148 (MQ=255)
                                                                                                       aTCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAg                                                        >  2:78917/1‑117 (MQ=255)
                                                                                                       aTCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAg                                                        <  1:78917/117‑1 (MQ=255)
                                                                                                            gTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGc                                                                                     <  2:124210/83‑1 (MQ=255)
                                                                                                            gTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGc                                                                                     >  1:124210/1‑83 (MQ=255)
                                                                                                             ttGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGga                  <  8:43980/149‑1 (MQ=255)
                                                                                                                             ccACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGaaaa  <  5:39655/149‑1 (MQ=255)
                                                                                                                                                cATACCGATGCGTTTCGGGGTTAAATCCGCAAACAg                                                                                                <  7:56233/36‑1 (MQ=255)
                                                                                                                                                cATACCGATGCGTTTCGGGGTTAAATCCGCAAACAg                                                                                                >  8:56233/1‑36 (MQ=255)
                                                                                                                                                   |                                                                                                                              
GTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAA  >  NC_000913/4440694‑4440967

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