New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi 3375666 =19 (1.040)4 (0.230) 4/272 1.9 18.6% coding (798/1368 nt) dcuD putative transporter
?NC_000913_3_vch_pgi 3375686 = 17 (0.980)coding (818/1368 nt) dcuD putative transporter
Rejected: Frequency below/above cutoff threshold.

GCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCAA  >  NC_000913_3_vch_pgi/3375520‑3375812
                                                                                                                                                  |                                                                                                                                                  
gacacacagcgttaacacgtgggcaaactgctgacccatgccgtcaaaaaacgcctgcacatcgtccattgtttcgcgcaagttatgcttgcggaagaactcaacaaacatcgtcacagtcaaactcagtaacatcaccaccaccaGAt                                                                                                                                                  >  2:48005‑M1/147‑149 (MQ=255)
   tttGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATACTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCa                                                                                                                                               <  1:208487/149‑1 (MQ=255)
     ccagcgttaccacgttggcaaactgcgtacccatgccgtcaaaaaacgcctgcacatcgtccattgtttcgcgcaagtaatgcttgcggaagaactgaacaaacatcgtcacagtcaaactcagtaacatcaccaccaccaGATGCAGt                                                                                                                                             <  1:299042‑M1/8‑1 (MQ=255)
            aaaGATAAGAATATCGATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACt                                                                                                                                      <  2:193683/149‑1 (MQ=255)
                 tAAAACTATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATc                                                                                                                                 <  1:21403/149‑1 (MQ=255)
                           aatcaCGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTggtggtggt                                                                                                                       <  2:111458/149‑1 (MQ=255)
                                      ggCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGa                                                                                                            <  1:197447/149‑1 (MQ=255)
                                        cAGAGCAAAAAGCTCTCGATAATGTCCCGGCGCTCTATTACGCCATTTTACCTGCGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGt                                                                                                          <  1:278390/149‑1 (MQ=255)
                                             aaaaaacgcctgcacatcgtccattgtttcgcgcaagttatgcttgcggaagaactcaacaaacatcgtcacagtcaaactcagtaacatcaccaccaccaGATGCAGTCAGAACTGCATCtgg                                                                                                                              >  1:243026‑M1/102‑124 (MQ=255)
                                             aaaaaacgcctgcacatcgtccattgtttcgcgcaagttatgcttgcggaagaactcaacaaacatcgtcacagtcaaactcagtaacatcaccaccaccaGATGCAGTCAGAACTGCATCtgg                                                                                                                              <  2:243026‑M1/23‑1 (MQ=255)
                                                                  ccgccgCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGtt                                                                                >  2:186159/1‑149 (MQ=255)
                                                                                cGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCCGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCGTa                                                                  >  1:303368/1‑149 (MQ=255)
                                                                                   cATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGAGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTGTGTTGAGTTCTTCCGCAAGCATAACt                                                               >  1:210261/1‑149 (MQ=255)
                                                                                       ttATCTGTGATGCCGTTAAGCCGGATGCTCGGCTCGCTGTTCCTCGTCCACGTCGGGCTTATGCAGATAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGCCGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTgcg                                                           <  2:159430/149‑1 (MQ=255)
                                                                                            tgtgATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTGTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGaaa                                                      >  1:200949/1‑149 (MQ=255)
                                                                                                           ccTGATGCTAGTTTCGATGGACCTCGCCAACGTCGGGCCGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGc                                       <  1:140703/149‑1 (MQ=255)
                                                                                                            cTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGcg                                      >  1:99042/1‑148 (MQ=255)
                                                                                                            cTGATGCCCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCa                                      <  1:255360/149‑1 (MQ=255)
                                                                                                                           actcagtaacatcaccaccaccaGATGCAGTCAGAACTGCATCtggtg                                                                                                                            <  1:25605‑M1/25‑1 (MQ=255)
                                                                                                                           actcagtaacatcaccaccaccaGATGCAGTCAGAACTGCATCtggtg                                                                                                                            >  2:25605‑M1/24‑48 (MQ=255)
                                                                                                                              ttCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCAt                    <  1:186159/149‑1 (MQ=255)
                                                                                                                                 cTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATggg                 >  1:117323/1‑149 (MQ=255)
                                                                                                                                        aCGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAg          <  1:304682/149‑1 (MQ=255)
                                                                                                                                               gCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCATCCTACTATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCa   <  1:144385/149‑1 (MQ=255)
                                                                                                                                                cTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGCGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCCGTTTGccca  >  2:185917/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATGTGCAGGCGTTTTTTGACGGCATGGGTACGCAGTTTGCCAA  >  NC_000913_3_vch_pgi/3375520‑3375812

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.