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.

AACATTACCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGATTCA  >  NC_000913_3_vch_pgi/1712661‑1712945
                                                                                                                                              |                                                                                                                                              
aaCATTACCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGCAAGCACCCCCGTTAATATGGGATGTAAAAAACGTGGTAACAGTGTCCGCTGCACACCAAAAACCAACTTAAAGCGTCAGttt                                                                                                                                          >  1:1428/1‑149 (MQ=255)
              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)
                                                     taCATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGt                                                                                                                                            >  2:236903/1‑94 (MQ=255)
                                                     taCATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGt                                                                                                                                            <  1:236903/94‑1 (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)
                                                                                                                                              |                                                                                                                                              
AACATTACCTGCCGTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGATTCA  >  NC_000913_3_vch_pgi/1712661‑1712945

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.