New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 390790526 (0.700)4 (0.120) 4/238 4.3 14.4% coding (462/774 nt) pstB phosphate ABC transporter ATPase
?NC_000913 = 3907917 23 (0.660)coding (450/774 nt) pstB phosphate ABC transporter ATPase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCG  >  NC_000913/3907762‑3908043
                                                                                                                                               |                                                                                                                                          
tCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCagagaga                                                                                                                                         >  1:44633/1‑147 (MQ=255)
               gTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTg                                                                                                                        >  2:132416/1‑149 (MQ=255)
                              aTACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTg                                                                                                                           >  4:93702/1‑131 (MQ=255)
                              aTACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTg                                                                                                                           <  3:93702/131‑1 (MQ=255)
                                                gggTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCgg                                                                                        <  2:44633/148‑1 (MQ=255)
                                                  gTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCaa                                                                                                                     >  3:137516/1‑117 (MQ=255)
                                                  gTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCaa                                                                                                                     <  4:137516/117‑1 (MQ=255)
                                                  gTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCtt                                                                                                             >  8:96173/1‑125 (MQ=255)
                                                  gTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCtt                                                                                                             <  7:96173/125‑1 (MQ=255)
                                                  gTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCtt                                                                                     <  1:132416/149‑1 (MQ=255)
                                                      aGCGCCGATCACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGt                                                                                 <  4:79956/149‑1 (MQ=255)
                                                       gcgcCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCagagaga                                                                                                                                         <  8:57765/92‑1 (MQ=255)
                                                       gcgcCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCagagaga                                                                                                                                         >  7:57765/1‑92 (MQ=255)
                                                                    ttCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCa                                                                   >  3:131101/1‑149 (MQ=255)
                                                                       gtcgAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCAcgc                                                                >  3:4789/1‑149 (MQ=255)
                                                                            gcagcagcaCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCg                                                           >  3:72066/1‑149 (MQ=255)
                                                                              agcagcaCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTgg                                                                                  >  5:160667/1‑124 (MQ=255)
                                                                              agcagcaCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTgg                                                                                  <  6:160667/124‑1 (MQ=255)
                                                                               gcagcaCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTc                                                         >  2:93706/1‑148 (MQ=255)
                                                                                    aCTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTc                                                   <  4:4789/149‑1 (MQ=255)
                                                                                     cTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCg                                                  >  8:106707/1‑149 (MQ=255)
                                                                                     cTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCg                                                  <  4:131101/149‑1 (MQ=255)
                                                                                                aTGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGagag                                       <  4:72066/149‑1 (MQ=255)
                                                                                                                            gCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGga                                          <  4:21615/118‑1 (MQ=255)
                                                                                                                            gCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGga                                          >  3:21615/1‑118 (MQ=255)
                                                                                                                                        caccaGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCg  >  7:23117/1‑146 (MQ=255)
                                                                                                                                               |                                                                                                                                          
TCCTGCTTCAGTTCGGTGATCAGCTCTTCAATACGCCCGGTAGAGATAGGGTCGAGCGCCGAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGAGAGAGTAACCGCTCTGGTGCAATTTATCTTTGGTTTCGTTCCACAATGCGGCTTTGGTCAATGCCCACTGCACGCGCTCGTCCATGTCGGCACGGGAGAGCTTCTCAAACAGACGAACGCCAAAAGCGATGTTGTCG  >  NC_000913/3907762‑3908043

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