New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 106135730 (0.860)4 (0.120) 4/252 3.9 12.3% coding (2102/2547 nt) torA trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
?NC_000913 = 1061374 29 (0.890)coding (2119/2547 nt) torA trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGCGCGCGCGGTATTCGTAACGGTGATGTGGTACGCGTC  >  NC_000913/1061215‑1061484
                                                                                                                                              |                                                                                                                               
tACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTc                                                                                                                           >  1:139594/1‑149 (MQ=255)
    cGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAg                                                                                                                       >  3:1418/1‑149 (MQ=255)
        aaCTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTAt                                                                                                                   <  5:98358/149‑1 (MQ=255)
          cTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTAtgt                                                                                                                 <  4:1418/149‑1 (MQ=255)
                          cTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGc                                                                                                 <  3:52313/149‑1 (MQ=255)
                              gacgaTTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGt                                                                                               >  1:141697/1‑147 (MQ=255)
                                cgaTTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTg                                                                                                  <  7:139065/142‑1 (MQ=255)
                                cgaTTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTg                                                                                                  >  8:139065/1‑142 (MQ=255)
                                      tCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAatat                                                                                      <  2:139594/148‑1 (MQ=255)
                                                         tttGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGc                                                                   >  6:28745/1‑148 (MQ=255)
                                                                           tCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCg                                                  >  4:28567/1‑147 (MQ=255)
                                                                           tCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAAACTGTGAATCCGGATTTCCGACTTCCATCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATCGGGTAGCGGGTAAAGAGCCAGTTTTCATTAACCCg                                                  >  7:44233/1‑147 (MQ=255)
                                                                              cACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGc                                                                              >  8:31574/1‑116 (MQ=255)
                                                                              cACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGc                                                                              <  7:31574/116‑1 (MQ=255)
                                                                                        cTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGt                                                                                 >  3:82630/1‑103 (MQ=255)
                                                                                        cTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGt                                                                                 <  4:82630/103‑1 (MQ=255)
                                                                                        cTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAgcgcg                                   >  3:8085/1‑149 (MQ=255)
                                                                                                  aaaGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGCGCGCGCGGTATTc                         <  5:28745/149‑1 (MQ=255)
                                                                                                      tATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGCGCGCGCGGTATTCGTaa                     <  3:28567/149‑1 (MQ=255)
                                                                                                           gTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGcgcgcgcg                               <  5:19277/134‑1 (MQ=255)
                                                                                                           gTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGcgcgcgcg                               >  6:19277/1‑134 (MQ=255)
                                                                                                           gTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGCGCGCGCGGTATTCGTAACGGt                 >  1:114117/1‑148 (MQ=255)
                                                                                                              gCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGCGCGCGCGGTATTCGTAACGGTGAt              <  2:141697/148‑1 (MQ=255)
                                                                                                                aTCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGAt                                            <  4:69575/116‑1 (MQ=255)
                                                                                                                aTCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGAt                                            >  3:69575/1‑116 (MQ=255)
                                                                                                                 tCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGCGCGCGCGGTATTCGTAACGGTGATGTgg          >  3:82400/1‑149 (MQ=255)
                                                                                                                       aTCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGCGCGCGCGGTATTCGTAACGGTGAtgt            <  4:95372/141‑1 (MQ=255)
                                                                                                                       aTCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGCGCGCGCGGTATTCGTAACGGTGAtgt            >  3:95372/1‑141 (MQ=255)
                                                                                                                            tgCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGCGCGCGCGGTATTCGTAACGGTGATGTGGTACGCGTc  <  4:82400/146‑1 (MQ=255)
                                                                                                                            tgCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGCGCGCGCGGTATTCGTAACGGTGATGTGGTACGCGTc  <  8:44233/146‑1 (MQ=255)
                                                                                                                                              |                                                                                                                               
TACTCGAAAACTATCGCCGATATGAACTACGACGATTGTCAGGGGCATCCGATGTGGTTTGAGAAAATCGAACGCTCCCACGGTGGGCCTGGCTCGCAAAAGTATCCGTTGCATCTGCAATCTGTGCATCCGGATTTCCGACTTCACTCGCAGTTATGTGAGTCGGAAACGCTGCGTCAGCAATATACGGTAGCGGGTAAAGAGCCAGTATTCATTAACCCGCAGGATGCCAGCGCGCGCGGTATTCGTAACGGTGATGTGGTACGCGTC  >  NC_000913/1061215‑1061484

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