New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 782739	9 (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

GCAAAAACAGACGGGTGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAG  >  NC_000913/782639‑782910                                                                                                                               |                                                                                                                                                  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)                                                                                          aagaagAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCgg                                    >  2:24990/1‑149 (MQ=255)                                                                                                     ccGGATGCTGCCATACTGGTGCATCCGggtattcttcttgcaagcgggttaacgcctgagtcttaaattcatcatgcacaatacaggcaccctgccagcatagaatgtctccacccgtctgtttttgcacgtaacgccccagatgtttg                         <  1:50399‑M2/149‑123 (MQ=255)                                                                                                                     tGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCg         >  2:72896/1‑149 (MQ=255)                                                                                                                            tCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAg  >  2:102108/1‑149 (MQ=255)                                                                                                                               |                                                                                                                                                  GCAAAAACAGACGGGTGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAG  >  NC_000913/782639‑782910

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