New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2121328 (1.040)3 (0.120) 3/268 3.6 10.3% coding (33/219 nt) yaaY uncharacterized protein
?NC_000913 = 21294 26 (1.030)coding (114/219 nt) yaaY uncharacterized protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATC  >  NC_000913/21167‑21434
                                                                                                                               |                                                                                                                                            
aaTAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATg                                                                                                                               >  1:76669/1‑143 (MQ=255)
    ggcgcCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGCTTCTTCCCGTGATTTGCGCCATGGCAATGaaaa                                                                                                                           <  1:309321/142‑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)
                                               ggATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGCACTCAATCCCCGGTCCACCtt                                                                                >  2:382339/1‑143 (MQ=255)
                                                gATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCATTCGCCGGTTAACCTTg                                                                               >  2:405676/1‑143 (MQ=255)
                                                  ttCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGAc                                                                             >  1:373881/1‑143 (MQ=255)
                                                    ctcccAGCTTGCGGGGTTTGAATACAGTTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCg                                                                           <  2:440787/140‑1 (MQ=255)
                                                      gccAGCTTGCGGGGTGTGCATACAGCTGTGCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCt                                                                         <  1:24617/142‑1 (MQ=255)
                                                             tggggggTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTGCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAg                                                                  <  2:220286/140‑1 (MQ=255)
                                                             tGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTGCCCGTGATTTGCGCCAGGGCAATGAAAAGCTATTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAg                                                                  <  1:336266/143‑1 (MQ=255)
                                                                ggggTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGta                                                               <  1:238890/143‑1 (MQ=255)
                                                                   gtgtGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATAc                                                            <  2:242937/143‑1 (MQ=255)
                                                                                     cgcgATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATAGCCGCATAACCTTGACTGTTGTACAAACTATACTCGGACGGTTTACACTGt                                          >  2:119530/1‑143 (MQ=255)
                                                                                              aaTTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACCCAATCGCCGGTTAGCCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCa                                 <  1:435601/143‑1 (MQ=255)
                                                                                                         gCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCGACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCt                      <  2:213844/143‑1 (MQ=255)
                                                                                                            gTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGCACGATTTTCCCTCTTCTGAGCCAGACATGAGCCTGag                   >  2:270315/1‑142 (MQ=255)
                                                                                                             tCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACGCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTGTGAGCCAGAGAGGAAGcaggta                  >  1:57331/1‑138 (MQ=255)
                                                                                                                  ttctgcCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCgg             <  2:5577/143‑1 (MQ=255)
                                                                                                                   tcttcCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGc            >  2:121690/1‑143 (MQ=255)
                                                                                                                    cttcCCTTTTTTTGCTCTTTGTCAATGCAATTCTACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGca           <  2:362305/143‑1 (MQ=255)
                                                                                                                             aTTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAAtc  <  1:386108/143‑1 (MQ=255)
                                                                                                                               |                                                                                                                                            
AATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATC  >  NC_000913/21167‑21434

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