New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA 3190438 =21 (0.880)4 (0.180) 4/274 2.7 17.3% coding (558/750 nt) yqiH putative periplasmic pilin chaperone
?NC_000913_3_hsa_tpiA 3190453 = 18 (0.790)coding (573/750 nt) yqiH putative periplasmic pilin chaperone
Rejected: Frequency below/above cutoff threshold.

TTCTGGCGGCCAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTA  >  NC_000913_3_hsa_tpiA/3190292‑3190566
                                                                                                                                                  |                                                                                                                                
aaatgttatttactacagggcagtttatttcataaaaacttaatccaccaaagtcattgatataacccacctgaaaagttgacggaacatttttaaactgcgtatttgcatagctaaaaggggttaccattaagctttgcgtttgtGGt                                                                                                                                <  1:329222‑M1/3‑1 (MQ=255)
  atgttatttactacagggcagtttatttcataaaaacttaatccaccaaagtcattgatataacccacctgaaaagttgacggaacatttttaaactgcgaatttgcatagctaaaaggggttaccattaagctttgcgtttgtGGTaa                                                                                                                              <  1:189066‑M1/5‑1 (MQ=255)
   tGGCGGCCAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAAc                                                                                                                             >  2:282309/1‑149 (MQ=255)
     gcggcCAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGa                                                                                                                           <  2:386046/149‑1 (MQ=255)
      cggcCAAAAGCCATTGAATTAGAAGACGATGAGGGCATGACTTACGAAAAAGTTGAAATCATTTGTAGAAATGACGGGTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGTGGTAAcatc                                                                                                                          <  1:97077/149‑4 (MQ=255)
                       aTTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCaaa                                                                                                         <  2:116470/149‑1 (MQ=255)
                        ttAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAg                                                                                                        >  1:391836/1‑149 (MQ=255)
                                 gATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACACTAAGCTGCCTTATCACGTTACGCTGGGTTGTATCGGAACAACTGTGGTACCGATGTTGCCACGAACGCAGAGCTTGACGGt                                                                                               >  1:185567/1‑149 (MQ=255)
                                 gATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATACCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGt                                                                                               >  1:184799/1‑149 (MQ=255)
                                 aaaaacttaatccaccaaagtcattgatataacccacctgaaaagttgacggaacatttttaaactgcgtatttgcatagctaaaaggggttaccattaagctttgcgtttggGGTAACGATGTTAc                                                                                                                     <  2:395928‑M1/14‑1 (MQ=255)
                                 aaaaacttaatccaccaaagtcattgatataacccacctgaaaagttgacggaacatttttaaactgcgtatttgcatagctaaaaggggttaccattaagctttgcgtttggGGTAACGATGTTAc                                                                                                                     >  1:395928‑M1/114‑127 (MQ=255)
                                   tGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGtt                                                                                                                       <  2:247376/123‑1 (MQ=255)
                                   tGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGtt                                                                                                                       >  1:247376/1‑123 (MQ=255)
                                            taCTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTa                                                                                    <  1:69748/148‑1 (MQ=255)
                                                     aaaaGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTTTGCaa                                                                           >  1:146930/1‑149 (MQ=255)
                                                     aaaaGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTGATAACGTAAAGCTGGGTTATATCGGAACAAATGGGGTAACAATGTAACAACAAACACAAAGCTTAATGTCAACCCCTTTTAGCTATGCaa                                                                           >  2:362087/1‑149 (MQ=255)
                                                       attgatataacccacctgaaaagttgacggaacatttttaaactgcgtatttgcatagctaaaaggggttaccattaagctttgcgtttgtGGTAACGATGTTACCACAAACGCAAAGCTTa                                                                                                    >  2:219834‑M1/92‑122 (MQ=255)
                                                       attgatataacccacctgaaaagttgacggaacatttttaaactgcgtatttgcatagctaaaaggggttaccattaagctttgcgtttgtGGTAACGATGTTACCACAAACGCAAAGCTTa                                                                                                    <  1:219834‑M1/31‑1 (MQ=255)
                                                         gTTGAAATCATTCGTAGAAATGACGGTTCGAACCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGGTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATAc                                                                       <  1:400916/149‑1 (MQ=255)
                                                                      gTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGACCAACTCGGTTTACGATGGTGCTATAAATTCAAAGCTTAATGGTAACCCCTTTTAACTATCCAAATACGCAGTTTCAAAAt                                                          >  2:361656/1‑149 (MQ=255)
                                                                                       aTCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATAGGCAGTTTAAAAATGTTCCGTCAACTTTGCa                                         >  1:419261/1‑149 (MQ=255)
                                                                                        tCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAGATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAg                                        >  1:277546/1‑149 (MQ=255)
                                                                                           gCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCTTCCACTTTTCAGGTg                                     >  1:35048/1‑149 (MQ=255)
                                                                                                    aTAAGATGCCTTATCACGTTCCGCCGGGTTATATCGGAACAAATGGGGCCACGACGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATc                            <  1:267869/149‑1 (MQ=255)
                                                                                                                ataaCGTTACGCTGGGTTATATCAGAACACATGGGGTAACGATTTTACTACAAACGTAAAGCTTAATGGTAACCCCGTTTATCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTAGCATTTGGGTTATATCAATGACTTtcgt                <  2:260026/146‑1 (MQ=255)
                                                                                                                           cTGGGTTATATCGGAACAAATTGGGTAACGATGTTACCAAAAACGCAAAGCTTAATGGTACCCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGttt     <  2:289600/149‑1 (MQ=255)
                                                                                                                              ggTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGTCTTTCGTGGATTAAGTTTTTa  >  2:234520/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                
TTCTGGCGGCCAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAATACGCAGTTTAAAAATGTTCCGTCAACTTTTCAGGTGGGTTATATCAATGACTTTGGTGGATTAAGTTTTTA  >  NC_000913_3_hsa_tpiA/3190292‑3190566

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.