New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 22497328 (0.450)10 (0.170) 8/268 6.5 27.2% noncoding (1203/1542 nt) rrsH 16S ribosomal RNA of rrnH operon
?NC_000913_3_hsa_tpiA 2729972 = NA (NA)noncoding (1186/1542 nt) rrsG 16S ribosomal RNA of rrnG operon
Rejected: Coverage evenness skew score above cutoff.

GCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCG  >  NC_000913_3_hsa_tpiA/224828‑225038
                                                                                                                                                  |                                                                 
gCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCGAAGGAGATTGCCAGTGAGAAACTGACGGAAGGTCGCGACGACGTGCAGTCatgatg                                                                >  2:674683/1‑149 (MQ=11)
       tcAGCTCATGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTAATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAACGGGACTGCCAGTGACAAACCGGGGGAAGGCGGGGAGGACGGCAAGTCATCATGGCCCtaa                                                         >  2:890351/1‑147 (MQ=11)
            tCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAACCTGGAGGAAGGTGGGGCTGACGTCACGTCATCATGGCACTTCGGAcc                                                    >  2:665772/1‑149 (MQ=11)
             cGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCa                                                   <  1:545157/149‑1 (MQ=32)
              gggttgtGAAAGGGTAGGTTAAGGCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAg                                                  <  2:890686/147‑1 (MQ=255)
                    tgAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTac                                            >  2:402421/1‑149 (MQ=32)
                      aaaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCCTCATGGCCCTTACGACCAGGGCTccac                                          >  2:133006/1‑149 (MQ=17)
                      aaaTGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTacac                                          >  2:366550/1‑149 (MQ=32)
                         tGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACg                                       <  2:1087272/149‑1 (MQ=32)
                           ttGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCTAGTGATATTCTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTg                                     <  1:58310/149‑1 (MQ=12)
                           ttGGGGTGAGTCCCGCAACGAGCGCACACTCTTTCTTTTGTTGCCAGGGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTg                                     <  2:883027/149‑1 (MQ=11)
                                    gTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGACGGGAACTCAAAGGAGACTGACAGTGATAAACTGGAGAGAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGATCAGGGCTACACCCGGGCTACCATgg                            >  1:989913/1‑149 (MQ=11)
                                     tCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGgc                           <  1:402421/149‑1 (MQ=32)
                                       ccGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAACCTGGAGGAAGGTGGGGATGACGTCAAGTCATCCTGGCCCTTACGACCAGGGCTACACACGTGCCACAATGgcgc                         >  2:669202/1‑149 (MQ=12)
                                       ccGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGCCTACACACGTGCTACacggccgc                         >  2:536031/1‑142 (MQ=12)
                                        cGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACCCGTGCTACAAGGGCGCa                        >  2:565414/1‑149 (MQ=17)
                                        cGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCa                        <  2:146086/149‑1 (MQ=32)
                                              gAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACaaa                  <  2:572060/149‑1 (MQ=32)
                                              gAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACaaa                  <  1:366550/149‑1 (MQ=32)
                                               aGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAag                 <  1:440758/149‑1 (MQ=32)
                                                 cgcAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACAAACGTGCTACAATGGCACATACAAagag               >  2:349849/1‑149 (MQ=17)
                                                  gcAACCCTTATCTTTTATTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGGCTGCCAGTGATAAACTGGAGGAGGGTGGGGATGACGTCAAGTCATCATGGCCCTTGCGACCAGGCCTACACACCTTCTACCCTTGCGCATACAAAgaga              >  2:289082/1‑149 (MQ=11)
                                                      atccctATTTTTTATTGTCAGCGCGCCTTCCGGGAATTTAAAGGAGACTGCCAGTGTTAAACTGGAGTAAGGTGTCTATGACGTCAAGTCATCATGGCCCTTACTACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGc           <  2:306693/144‑1 (MQ=11)
                                                      aCCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGc           <  2:697325/149‑1 (MQ=32)
                                                        ccTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGa         >  1:515800/1‑149 (MQ=32)
                                                         cTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGAc        <  1:906537/149‑1 (MQ=32)
                                                           atacTTTTGTGCCGGCGGTCCGGCCGGTAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGCCATCCCGATCTTCACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTc     <  1:213834/146‑1 (MQ=11)
                                                               ttttGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTcgcg  >  1:1000238/1‑149 (MQ=11)
                                                                                                                                                  |                                                                 
GCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCTTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCG  >  NC_000913_3_hsa_tpiA/224828‑225038

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