New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi = 171280325 (0.790)4 (0.140) 4/270 4.0 14.4% coding (35/1503 nt) dtpA dipeptide and tripeptide permease A
?NC_000913_3_vch_pgi = 1712811 24 (0.810)coding (43/1503 nt) dtpA dipeptide and tripeptide permease A
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCT  >  NC_000913_3_vch_pgi/1712675‑1712958
                                                                                                                                        |                                                                                                                                                   
tttCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAaac                                                                                                                                         <  1:162078/149‑1 (MQ=255)
tttCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGCGTCCACTGCAAACCAAAAACCAACTGAACGCGTCAGTTTGACCGCTTTCAaac                                                                                                                                         >  1:127750/1‑149 (MQ=255)
    cccTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAAcc                                                                                                                                     >  1:263122/1‑149 (MQ=255)
                       gcgcGTAGGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCAtc                                                                                                                  <  1:402809/149‑1 (MQ=255)
                            tATGTTTTTTATACATAACGTGAGAGAGAAACCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCTCTGCAAACTAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTGCAAACAACCGAAGGCGTTCTATCTCATCTtctc                                                                                                             <  2:183571/149‑1 (MQ=255)
                                ttttttATACATAACGCGAGAAAGCACCCCCGTTAATATGCGATGTAAAAAAAGAGGTAAAAGTATCCACTGCTAACCAAAAACCAACTGAAAGCGTCAGTTTGCACGCTTTCAAACAACAGAAGGCATTCTATCTCAtcttcttcatt                                                                                                         >  2:160677/1‑144 (MQ=255)
                                                   gaAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGAAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCCCGATTGAATTATGGGAACGTTTTg                                                                                      >  2:16932/1‑149 (MQ=255)
                                                         aCCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACCACCAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCATTCAAACAACAGAAGGCGTTCTATCACATCTCCACGCTTGAGTTATGGGAACGTCTTGCttatt                                                                                >  2:241142/1‑149 (MQ=255)
                                                         aCCCCCGTTAATATGGGATGTAAAAAAAGAGGGAAAAGTGTCCACTGCAAACAAAAAACCAACTGAAAGCGTCAGTTTGAGCGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGGGTTATGGGAACGTTTTGGttatt                                                                                >  1:244695/1‑149 (MQ=255)
                                                          cccccGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGttatta                                                                               >  2:22019/1‑149 (MQ=255)
                                                                       gggATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTAAGTTATGGGAACGTTTTTGTTATTACGGCCTACAAGGa                                                                  >  2:267840/1‑149 (MQ=255)
                                                                              aaaaaaaGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGGTTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATgg                                                           >  2:297591/1‑149 (MQ=255)
                                                                               aaaaaaGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGc                                                          >  1:418296/1‑149 (MQ=255)
                                                                                  aaaGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGt                                                       >  2:497587/1‑149 (MQ=255)
                                                                                        gTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCt                                                 <  1:149368/149‑1 (MQ=255)
                                                                                             aGTGTCCACTGCAAACCAAAAACCAACCGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATTGGAACGTTTTGGTTATCCCGGCCATTTTGTAATTATGGCTGTTTACCTGGTTa                                            <  2:120556/149‑1 (MQ=255)
                                                                                                     cTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCTTTTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTg                                    <  1:140255/149‑1 (MQ=255)
                                                                                                      tGCAAGCCTAAAACCAACTGAAAGCGTCAGTTTGAACGTTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCACGATTGAGTTATGGGAACGTTTTGGTTATCACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTgg                                   <  2:58146/149‑1 (MQ=255)
                                                                                                           accaaaaaccaaCTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATg                              >  1:62886/1‑149 (MQ=255)
                                                                                                                   ccaaCTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGCCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGc                      <  1:103817/149‑1 (MQ=255)
                                                                                                                      aCTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGa                   <  2:440140/149‑1 (MQ=255)
                                                                                                                          aaaGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTAGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGCTTCa               >  1:477109/1‑149 (MQ=255)
                                                                                                                                   gTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGCATTATGGCTGTTTACCTGGTTAAACAACTTGTGATGTCTGAAGCGGATTCAATAACCCtt      >  1:101461/1‑149 (MQ=255)
                                                                                                                                       gAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGCATGTCTGAAGCGGATTCACTCACCCGTTtcc  >  1:36564/1‑148 (MQ=255)
                                                                                                                                        |                                                                                                                                                   
TTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCT  >  NC_000913_3_vch_pgi/1712675‑1712958

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