New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2121332 (1.150)3 (0.110) 3/268 3.7 9.5% coding (33/219 nt) yaaY uncharacterized protein
?NC_000913 = 21294 27 (1.030)coding (114/219 nt) yaaY uncharacterized protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTAGCCAAGGTCCAACTCCCAAATGTGTTCTATATGGACAATTCAAAGGCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCT  >  NC_000913/21071‑21355
                                                                                                                                              |                                                                                                                                              
gcacccATGGTCCAACTCCCAAATGTGTGCTATATGGACAATTCAAAGGCCGAGGAATATGCCCTTTTAGCCGTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCa                                                                                                                                                <  1:325567/139‑1 (MQ=255)
 tAGCCAAGGTCCAACTCCCAAATGTGTTCTATATGGACAATTCAAAGGCCGAGGAATATGCCCTGTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAg                                                                                                                                               >  2:386108/1‑143 (MQ=255)
   gCCAGGGTCCAACTCCCAAATGTGGTCTATATGGACAAGGGAAAGGTCGAGGAATATGCTCTTTTAGCCTTCTTTTGGCAATGGATTGGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGa                                                                                                                                             <  1:275918/143‑1 (MQ=255)
   gCCAAGGTCCAACTCCCAAATGTGTTCTATATGGACAATTCAAAGGCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAGTGGCTTTGTGCAAATAAGCTCCGTTAATGTGCCGGCACTCGTTACGTAGTAATGGCACAAGa                                                                                                                                             >  1:242937/1‑143 (MQ=255)
   gCCAAGGTCCAACTCCCAAATGTGTTCTATATGGACAATTCAAAGGCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGAATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGa                                                                                                                                             >  2:238890/1‑143 (MQ=255)
    ccAAGGTCCAACTCCCAAATGTGTTCTATATGGACAATTCAAAGGCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAggct                                                                                                                                            >  1:5577/1‑141 (MQ=255)
             aaCTCCCAAATGTGTTCTATATGGACAATTCAAAGGCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTGGTGCAAATAAGCGCCGTTAATGTGCCGGCCCTCGTTACGTAGGGGTGGCGCAGGATTCTACGAGc                                                                                                                                   >  2:127421/1‑143 (MQ=255)
                  ccAAATGTGTTCTATATGGACAATTCAAAGGCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCGTGCg                                                                                                                              >  2:266093/1‑143 (MQ=255)
                    aaaTGTGTTCTATATGGACAATTCAAAGGCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGGATTGGGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTTGTGATGGCGCAGGGTTCTACCCGCTTGgggg                                                                                                                            >  1:24876/1‑143 (MQ=255)
                    aaaTGTGTGCTATATGGACAATTCAAAGGCCGAGGGATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCggg                                                                                                                            <  1:257750/143‑1 (MQ=255)
                      aTGTGTTCTATATGGACAATTCAAAGGCCGATGAATATGCCCTTTTACCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGgt                                                                                                                          <  2:437042/143‑1 (MQ=255)
                        ggcgcttttatGGGCGATTCAAGGGCCGAGGAATATGCCCTTTTAGCCGTCTTTTGTCAAGGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGGGATGGCGCAGGATTCTAGCAGCTTGCGGGgtgt                                                                                                                        <  1:96013/135‑1 (MQ=255)
                              tataTGGACAATTCAAAGGCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATAc                                                                                                                  <  1:273500/143‑1 (MQ=255)
                                                  cGAGGAATATGCCCTTTTAGCCTTCTTTTGTAAATGGATTTGTGCAAATAAGAGCCGTTAATGTGCCGGCACGCGTTACGTAGTGATGGCGCAGGATGGTACCAGCTTCCTGGGTGTGGATACAGCTTTTCcgccacataaac                                                                                              >  2:40111/1‑134 (MQ=255)
                                                    gggAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTg                                                                                            <  2:328260/142‑1 (MQ=255)
                                                    aGGAATATGCCCTTTTAGCCTTCTTTTGTCGATGGGTTTGTACAACTCGGTCGTGTGAATAGGCCGGCACTCGTTGCGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTg                                                                                            <  1:222857/143‑1 (MQ=255)
                                                    aGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTg                                                                                            <  2:331624/143‑1 (MQ=255)
                                                      gAATATGCCCTTTTAGCCTTCTTTTGTCAATGGAGTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTgca                                                                                          <  1:462425/143‑1 (MQ=255)
                                                           tgtccctttAGCCTTCCTTTGGCTAGGGATGTGTGCAAAGAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCATGATTCTACCAACTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAgg                                                                                     <  2:88975/137‑1 (MQ=255)
                                                                 tttAGCCTTCTTTTGTCAATGGATTTGTGCACCTACAGATCGATAGTGTGCCGTCTCTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAACATTGCAGCAGGCGGTCa                                                                               >  1:333054/1‑143 (MQ=255)
                                                                   ttggcttcttTTGTCAATGGGTTGGTGCATAGAAGCGCCTTTAATGTGCCGGCACTCGTGACGTAGTGATGGCGCAGTATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGt                                                                             <  2:179060/139‑1 (MQ=255)
                                                                      ccTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTttc                                                                          >  2:134811/1‑143 (MQ=255)
                                                                         gcttGTGTCAATGGATTTGTGCCAATACGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTttcttc                                                                       <  1:127421/142‑1 (MQ=255)
                                                                               gTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGa                                                                 <  2:24876/143‑1 (MQ=255)
                                                                                              cAAATAAGCGCCGTTAATGTGCCGGCACGCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGGGGGAATACAGCTTTTCCGCGATAAGAATGGCCGGACGAGGGCAGCTTATTCGCGTGAtgagccacatcgcac                                                  >  2:481368/1‑129 (MQ=255)
                                                                                                aaTAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATg                                                >  1:76669/1‑143 (MQ=255)
                                                                                                    ggcgcCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGCTTCTTCCCGTGATTTGCGCCATGGCAATGaaaa                                            <  1:309321/142‑1 (MQ=255)
                                                                                                        acgtagaggggCCGGCACCCGTTACGTAGGCATGGCGCAGGATCCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCa                                        <  1:40111/133‑1 (MQ=255)
                                                                                                                        aCTCGTTACGGAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCg                        <  1:340210/143‑1 (MQ=255)
                                                                                                                            gTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTAc                    >  2:318565/1‑143 (MQ=255)
                                                                                                                                  tAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTTTTCGCGTGATTTGCGTCATGGCAATGAAAAGCCACTTCTTTGTGCTTTCGGTACTCGATc              >  1:90256/1‑143 (MQ=255)
                                                                                                                                     tGATGGTGCAGGACTCTACCAGCTTGCGGGGTGTCAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGcc           <  2:362861/143‑1 (MQ=255)
                                                                                                                                         ggCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTTATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGtt       <  1:134811/143‑1 (MQ=255)
                                                                                                                                              aGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCt  >  1:30938/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                              
TTAGCCAAGGTCCAACTCCCAAATGTGTTCTATATGGACAATTCAAAGGCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCT  >  NC_000913/21071‑21355

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