New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 385820526 (0.680)4 (0.120) 4/240 4.3 13.0% coding (200/1494 nt) yidJ sulfatase/phosphatase superfamily protein
?NC_000913 = 3858228 30 (0.870)coding (177/1494 nt) yidJ sulfatase/phosphatase superfamily protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTC  >  NC_000913/3858101‑3858375
                                                                                                                               |                                                                                                                                                   
gCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTa                                                                                                                                >  8:4750/1‑149 (MQ=255)
      tGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAAc                                                                                                                           <  7:4750/148‑1 (MQ=255)
        ccGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCg                                                                                                                         <  5:109857/148‑1 (MQ=255)
             gTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCg                                                                                                                         >  6:58907/1‑143 (MQ=255)
             gTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCg                                                                                                                         <  5:58907/143‑1 (MQ=255)
                                     cTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGaa                                                                                           >  1:48355/1‑149 (MQ=255)
                                            tAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGgtgt                                                                                                             >  5:37257/1‑124 (MQ=255)
                                            tAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGgtgt                                                                                                             <  6:37257/124‑1 (MQ=255)
                                                           gagaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTc                                                                                    <  5:127314/134‑1 (MQ=255)
                                                           gagaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTc                                                                                    >  6:127314/1‑134 (MQ=255)
                                                           gagaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTAt                                                                      >  5:139244/1‑148 (MQ=255)
                                                             gaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCaa                                                                   >  6:178205/1‑149 (MQ=255)
                                                                 tttttGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGaa                                                                                           <  6:98075/121‑1 (MQ=255)
                                                                        cTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTAt                                                                      <  7:9727/135‑1 (MQ=255)
                                                                        cTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTAt                                                                      >  8:9727/1‑135 (MQ=255)
                                                                          ggCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTAtt                                                      <  2:48355/149‑1 (MQ=255)
                                                                             gcgACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAg                                                   >  3:68392/1‑149 (MQ=255)
                                                                               gACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCg                                                 <  6:139244/149‑1 (MQ=255)
                                                                                      ttgGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCg                                                 >  3:114202/1‑142 (MQ=255)
                                                                                      ttgGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCg                                                 <  4:114202/142‑1 (MQ=255)
                                                                                                     gACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGAc                           >  2:10563/1‑149 (MQ=255)
                                                                                                     gACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGAc                           >  8:65994/1‑149 (MQ=255)
                                                                                                          gTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCata                       <  5:178205/148‑1 (MQ=255)
                                                                                                                   gATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCt             <  1:10563/149‑1 (MQ=255)
                                                                                                                     tACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTgg           <  7:65994/149‑1 (MQ=255)
                                                                                                                          gTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATAtt                                                               >  8:9966/1‑92 (MQ=255)
                                                                                                                          gTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATAtt                                                               <  7:9966/92‑1 (MQ=255)
                                                                                                                          gTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATc      >  2:177527/1‑149 (MQ=255)
                                                                                                                           tAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCg     <  1:177527/149‑1 (MQ=255)
                                                                                                                              aTAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTc  <  4:68392/149‑1 (MQ=255)
                                                                                                                               |                                                                                                                                                   
GCCATTTGCCGATGTAACAGGTGTGATAGCCGGCATCCTTAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTC  >  NC_000913/3858101‑3858375

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