New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3988386 (0.890)6 (0.060) 6/252 NT 6.8% coding (504/1143 nt) caiA crotonobetaine reductase subunit II, FAD‑binding
?NC_000913 = 39906 83 (0.900)coding (481/1143 nt) caiA crotonobetaine reductase subunit II, FAD‑binding

TCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/39731‑39883
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgccTACACCCCGTACATCGTGGTGATGGCGCGCGACGGGGCTTCTCCGGACAAACCTGTCTACACCGAATGGTTTGTTGATATGAGCAAACCGGGCATCAAAGTGACCAAACTTGAAAAGCTCGGTCTGCGTATGGATAGCTGCTGTGAAA  <  NC_000913/39906‑39759
                                                                                                                                                                                                                                                                                                             
TCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCG                                                                                                                                                          >  2:643697/1‑149
 CGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGC                                                                                                                                                         <  1:516315/149‑1
   TCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGC                                                                                                                                                       >  2:180892/1‑149
    CCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCC                                                                                                                                                      <  1:180892/149‑1
                     AAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCTACACCCCGTACATCGT                                                                                                                                     >  1:200285/1‑149
                                ACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCTACACCCCGTACATCATGGTGATGGCGC                                                                                                                          >  1:203328/1‑149
                                ACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCTACACCCCGTACATCATGGTGATGGCGC                                                                                                                          <  2:203328/149‑1
                                        TATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCTACACCCCGTACATCGTGGTGATGGCGCGCGACGGG                                                                                                                  <  2:200285/149‑1
                                                                                                             ATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCTACACCCCGTACATCGTGGTGATGGCGCGCGACG                                                                                                                    >  3:140286/1‑78
                                                                                                             ATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCTACACCCCGTACATCGTGGTGATGGCGCGCGACG                                                                                                                    <  4:140286/78‑1
                                                                                                                                                        CTACACCCCGTACATCGTGGTGATGGCGCGCGACGGGGCTTCTCCGGACAAACCTGTCTACACCGAATGGTTTGTTGATATGAGCAAACCGGGCATCAAAGTGACCAAACTTGAAAAGCTCGGTCTGCGTATGGATAGCTGCTGTGAAA  <  2:860150/149‑1
                                                                                                                                                                                                                                                                                                             
TCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/39731‑39883
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgccTACACCCCGTACATCGTGGTGATGGCGCGCGACGGGGCTTCTCCGGACAAACCTGTCTACACCGAATGGTTTGTTGATATGAGCAAACCGGGCATCAAAGTGACCAAACTTGAAAAGCTCGGTCTGCGTATGGATAGCTGCTGTGAAA  <  NC_000913/39906‑39759

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