New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 7827399 (0.800)3 (0.270) 3/260 NT 25.1% coding (655/1044 nt) nadA quinolinate synthase, subunit A
?NC_000913 = 782765 9 (0.810)coding (681/1044 nt) nadA quinolinate synthase, subunit A

GCACCCGACAAACATCTGGGGCGTTACGTGCAAAAACAGACGGGTGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTA  >  NC_000913/782610‑782887
                                                                                                                                 |                                                                                                                                                    
gCACCCGACAAACATCTGGGGCGTTACGTGCAAAAACAGACGGGTGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACTCGGATGCTGCCATACTGGt                                                                                                                                   <  3:61628/149‑1 (MQ=255)
        cAAACATCTGGGGCGTTACGTGCAAAAACAGACGGGTGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCcggatgcaccagtatggcagcatccgg                                                                                                                           <  1:50399‑M1/149‑28 (MQ=255)
                             gCAAAAACAGACGGGTGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCg                                                                                                      >  4:39073/1‑149 (MQ=255)
                                                      cTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCTGGATGCACCAGTATGGCAGCATCCGggtattcttcttgcaagcgg                                                                                                        >  3:9081‑M2/1‑102 (MQ=255)
                                                      cTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCTGGATGCACCAGTATGGCAGCATCCGggtattcttcttgcaagcgg                                                                                                        <  4:9081‑M2/122‑21 (MQ=255)
                                                                        tGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCgctgc                                                           <  2:44008/149‑1 (MQ=255)
                                                                             tgtgCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGa                                                      <  3:39073/149‑1 (MQ=255)
                                                                                  atgatgAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACa                                                 >  3:31007/1‑149 (MQ=255)
                                                                                          tttAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATCTGGCGGATGCGGTCGGTTCCACCAGTCAATTGATCGCTGCTGCGAAAACATTGCcaca                                         >  4:87040/1‑149 (MQ=255)
                                                                                                  tCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGc                                 >  3:77604/1‑149 (MQ=255)
                                                                                                        gTTAACCCGCTTACAAGAACAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATtg                           <  2:77056/149‑1 (MQ=255)
                                                                                                              ccGCTTGCAAGAAGAATACCcggatgctgccatactggtgcatccaggtattcttcttgcaagcgggttaacgcctgagtcttaaattcatcatgcacaatacaggcaccctgccagcatag                                                >  3:9081‑M1/1‑20 (MQ=255)
                                                                                                              ccGCTTGCAAGAAGAATACCcggatgctgccatactggtgcatccaggtattcttcttgcaagcgggttaacgcctgagtcttaaattcatcatgcacaatacaggcaccctgccagcatag                                                <  4:9081‑M1/122‑103 (MQ=255)
                                                                                                                      aagaagAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCgg             >  2:24990/1‑149 (MQ=255)
                                                                                                                                 ccGGATGCTGCCATACTGGTGCATCCGggtattcttcttgcaagcgggttaacgcctgagtcttaaattcatcatgcacaatacaggcaccctgccagcatagaatgtctccacccgtctgtttttgcacgtaacgccccagatgtttg  <  1:50399‑M2/149‑123 (MQ=255)
                                                                                                                                 |                                                                                                                                                    
GCACCCGACAAACATCTGGGGCGTTACGTGCAAAAACAGACGGGTGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTA  >  NC_000913/782610‑782887

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