New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	= 1822923	13 (0.800)	3 (0.190)	3/242	1.9	21.1%	coding (466/828 nt)	nadE	NAD synthetase, NH3/glutamine‑dependent
	?	NC_000913_3_pae_pgi	= 1822933	10 (0.640)					coding (476/828 nt)	nadE	NAD synthetase, NH3/glutamine‑dependent

TTAATATCAAGGGCGCGGTATTGGCCAGCGAACAGGCATTGCGGGAAGCAGGCATTGAACTGAGCGATTTTGTCCGTGGCAATGAAAAAGCGCGTGAGCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTGGGCACCGATCATGCAGCAGAAGCCATTACCGGATTCTTCACTAAATATGGTGACGGCGGTACGGACATTAATCCGCTGTATCGTCTCAACAAACGTCAGGGTAAACAGTTACTGGCGGCAT  >  NC_000913_3_pae_pgi/1822783‑1823052                                                                                                                                             |                                                                                                                                  ttAATACCAAGGGCGCGGTATTGGCCAGCGAACAGGCATTGCGGGAAGCAGGCATTGAACTGAGCGATTTTGTCTGTGGCAATGAAAAAGCGCGTGAGCGGATGAAAGCACAATATAGCATTGCGGGTATGCCCAGCGGAGTCGtggtgg                                                                                                                          >  2:184916/1‑150 (MQ=255)                   tATTGGCCAGCGAACAGGCATTGCGGGAAGCAGGCATTGAACTGAGCGATTTTGTCCGTGGCAATGAAAAAGCGCGTGAGCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTggg                                                                                                                         >  1:140000/1‑133 (MQ=255)                   tATTGGCCAGCGAACAGGCATTGCGGGAAGCAGGCATTGAACTGAGCGATTTTGTCCGTGGCAATGAAAAAGCGCGTGAGCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTggg                                                                                                                         <  2:140000/133‑1 (MQ=255)                                aCAGGCATTGCGGGAAGCAGGCATTGAACTGAGCGATTTTGTCCGTGGCAATGAAAAAGCGCGTGAGCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGcccaccacgacaccgctggtcatacccgcaatgctatatt                                                                                           <  4:81554‑M1/150‑41 (MQ=255)                                  aGGCATTGCGGGAAGCAGGCATTGAACTTAGCGATTTTGTCCGTGGCAATGAAAAAGCGCGTGAGCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTGGGCACCGATCATGCAGCAGAAGCCATTACCGGAt                                                                                         <  3:89187/150‑1 (MQ=255)                                                         gAACTGAGCGATTTTGTCCGTGGCAATGAAAAAGCGCGTGAGCGGATGAAAGCTCAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTGGGCACCGATGATGCAGCAGAAGCCATTACCGGATTCTTGACTAAATATGGTGAcggc                                                                  <  2:125444/150‑1 (MQ=255)                                                             tGAGCGATTTTGTCCGTGGCAATGAAAAAGCGCGTGAGCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTGGGCACCGATCATGCAGCAGAAGCCATTACCGGATTCTTCACTAAATATGGTGACGGCGGTa                                                              >  1:77707/1‑150 (MQ=255)                                                                gCGATTTTGTCCGTGGCAATGAAAAAGCGCGTGAGCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTGGGCACCGATCATGCAGCAGAAGCCATTACCGGATTCTTCACTAAATATGGTGACGGCGGTACgg                                                           <  2:77707/150‑1 (MQ=255)                                                                           cGTGGCAATGAAAAAGCGCGTGAGCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTGGGCACCGATCATGCAGCAGAAGCCATTACCGGATTCTTCACTAAATATGGTGACGGCGGTACGGACATTAATCCg                                                <  3:96983/150‑1 (MQ=255)                                                                                       aaaGCGCGTGAGCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTGGGCACCGATCATGCAGCAGAAGCCATTACCGGATTCTTCACTAAATATGGTGACGGCGGTACGGACATTAATCCGCTGTATCGtctc                                    >  3:31830/1‑150 (MQ=255)                                                                                                  gCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGcccaccacgacaccg                                                                                                                    <  1:125454‑M1/59‑16 (MQ=255)                                                                                                  gCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTGGGCACCGATCATgcag                                                                                                           <  3:26570/68‑1 (MQ=255)                                                                                                  gCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTGGGCACCGATCATgcag                                                                                                           >  4:26570/1‑68 (MQ=255)                                                                                                  gCGGATGAAAGCACAATACAGCATTGCGGGTATGACCAGCGGTGcccaccacgacaccg                                                                                                                    >  2:125454‑M1/1‑44 (MQ=255)                                                                                                                aaTATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTGGGcaccgctggtcatacccgcaatgctatattgtgctttcatccgctcacgcgctttttcattgccacggacaaaatcgctcagttcaatgcctgcttcccgcaatgcctgt           <  4:81554‑M2/150‑111 (MQ=255)                                                                                                                                   gACCAGCGGTGTCGTGGTGGGCACCGATCATGCAGCAGAAGCCATTACCGGATTCTTCACTAAATATGGTGACGGCGGTACGGACATTAATCCGCTGTATCGTCTCAACAAACGTCAGGGTAAACAGTTACTGGCGGCAt  >  3:49124/1‑140 (MQ=255)                                                                                                                                   gACCAGCGGTGTCGTGGTGGGCACCGATCATGCAGCAGAAGCCATTACCGGATTCTTCACTAAATATGGTGACGGCGGTACGGACATTAATCCGCTGTATCGTCTCAACAAACGTCAGGGTAAACAGTTACTGGCGGCAt  <  4:49124/140‑1 (MQ=255)                                                                                                                                         cGGTGTCGTGGTGGGcaccgctggtcatacccgcaatgctgtattgtgctttcatccgc                                                                             >  2:125454‑M2/1‑15 (MQ=255)                                                                                                                                         cGGTGTCGTGGTGGGcaccgctggtcatacccgcaatgctatattgtgctttcatccgc                                                                             <  1:125454‑M2/59‑45 (MQ=255)                                                                                                                                             |                                                                                                                                  TTAATATCAAGGGCGCGGTATTGGCCAGCGAACAGGCATTGCGGGAAGCAGGCATTGAACTGAGCGATTTTGTCCGTGGCAATGAAAAAGCGCGTGAGCGGATGAAAGCACAATATAGCATTGCGGGTATGACCAGCGGTGTCGTGGTGGGCACCGATCATGCAGCAGAAGCCATTACCGGATTCTTCACTAAATATGGTGACGGCGGTACGGACATTAATCCGCTGTATCGTCTCAACAAACGTCAGGGTAAACAGTTACTGGCGGCAT  >  NC_000913_3_pae_pgi/1822783‑1823052

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