New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1570938	22 (0.210)	8 (0.080)	6/170	NT	16.7%	coding (1108/1401 nt)	gadB	glutamate decarboxylase B, PLP‑dependent
	?	NC_000913	= 3666449	59 (0.600)					coding (1132/1401 nt)	gadA	glutamate decarboxylase A, PLP‑dependent

CGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGCCCCAGTTTGGCGATTTCATCCGCCAGATAAGCGGCAACCTGGTAAGAGGCGTTCTGTACTTTGGTATAGCCTTC  >  NC_000913/1570841‑1571037                                                                                                  |                                                                                                    cgcagacgcagaCGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGAc                                                                                                   >  1:1383088/1‑100 (MQ=33)  gcagacgcagacGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGAc                                                                                                   <  1:1338659/99‑1 (MQ=32)     gacgcagacgTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGt                                                                                               >  1:1873563/1‑100 (MQ=33)       cgcagacgTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTAc                                                                                             >  1:2805947/1‑100 (MQ=33)        gcagacgTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACa                                                                                            <  1:182530/100‑1 (MQ=33)         cagacgTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAg                                                                                           <  1:1252164/100‑1 (MQ=33)         cagacgTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAg                                                                                           <  2:1461159/100‑1 (MQ=33)             cgTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGa                                                                                       >  2:1949965/1‑100 (MQ=33)             cgTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGa                                                                                       >  2:823107/1‑100 (MQ=33)               ttCAGAGAGGTCATACAGGGTNTATCCCGGATCTTCACNATCTTTCAGTTTGAAGCAAACCGCCGGNATGCCTTCGTCCGGGCGACCCGTACAGATGAAc                                                                                     >  2:2558638/1‑100 (MQ=18)               ttCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAAc                                                                                     <  1:187032/100‑1 (MQ=32)               ttCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAAc                                                                                     <  1:1382660/100‑1 (MQ=32)                tCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGTTGCCTTCGTCCGGGCGACCCGTACCGATGAAc                                                                                     >  1:2693872/1‑99 (MQ=16)                 cAGAGAGGTCATACAGGGTGTATCCCGGCTCTTCACCATCTTTCAGTTTGAAGCACACCGCCGGGATGCCTTCGTCCTGGCGACCCGTACAGATGAACTc                                                                                   >  1:913653/1‑100 (MQ=11)                   gagagGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCAt                                                                                 >  1:2693367/1‑100 (MQ=33)                   gagagGTCATACAGGGTGTAGCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCAt                                                                                 >  1:1772370/1‑100 (MQ=17)                    agagGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATa                                                                                >  2:474348/1‑100 (MQ=33)                     gagGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATAc                                                                               <  1:1677791/100‑1 (MQ=33)                     gagGTCATACAGGGTGTATCCCAGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATAc                                                                               >  1:1684537/1‑100 (MQ=18)                       ggTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACgg                                                                             >  1:2609061/1‑100 (MQ=32)                       ggTCATACAGGGTGTATCCAGGATCTTCACGATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGCACTCCTACg                                                                              >  1:1106312/1‑99 (MQ=11)                          cATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGccc                                                                          >  1:454330/1‑100 (MQ=31)                          cATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGccc                                                                          >  1:1383188/1‑100 (MQ=31)                          cATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGGTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGgtttgcttcaaactgaaagatggtga                                                                           >  1:1560935‑M1/1‑73 (MQ=255)                           aTACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGcccc                                                                         <  1:55992/100‑1 (MQ=11)                                                             aGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGgtttgcttcaaactgaaagatggtgaagatccgggatacaccctgtacgacctctctgaacg                                       >  2:2014519‑M1/1‑38 (MQ=255)                                                              gTTTGAAGCAAACCGCCGGGATGCATTCGTCCGGGCGgtttgcttcaaactgaaagatggtgaagatccgggatacaccctctatcagcgctctgaacg                                       >  1:708932‑M1/1‑37 (MQ=255)                                                                   aaGCAAACCGCCGGGATGCCTTCGTCCGGGCGgtttgcttcaaactgaaagatggtgaagatccgggagacaccctgtacaacctctctgaacgtctgcg                                 >  2:1238764‑M1/1‑32 (MQ=255)                                                                     gCAAACCGCCGGGATGCCTTCGTCCGGGCGgtttgcttcacactgaaagatggtgaagatccgggatacacactgtatgacctctctgaacgtctgcg                                 >  1:679343‑M1/1‑30 (MQ=255)                                                                      cAAACCGCCGGGATGCGTTCGTCCGGGCGgtttgcttcaaactgaaagctggtgacgatccgggatacaccctgtacgacctctctgaacgtctgcatct                              >  1:763324‑M1/1‑29 (MQ=255)                                                                          ccgccgGGATGCCTTCGTCCGGGCGgtttgcttcaaactgaaagatggtgaagatccgggatacaccctgtacgacctctctgaacgtctgcgtctgcgc                          >  2:1850564‑M1/1‑25 (MQ=255)                                                                             ccgGGATGCCTTCGTCCGGGCGgtttgcttcaaactgaaagatggtgaagatccgggatacaccctgtacgacctctctgaacgtctgcgtctgcgcggc                       >  2:853787‑M1/1‑22 (MQ=255)                                                                                               ggCGgtttgcttcaaactgaaagatggtgaagatccgggatacaccctgtacgacctctctgaacgtctgcgtctgcgcggctggcaggttccggccttc     <  2:2167230‑M1/100‑97 (MQ=255)                                                                                               ggCGgtttgcttcaaactgaaagatggtgaagatccgggatacaccctgtacgacctctctgaacgtctgcgtctgcgcggctggcaggttccggccttc     <  2:2889784‑M1/100‑97 (MQ=255)                                                                                               ggCGgtttgcttcaaactgaaagatggtgaagatccgggatacaccctgtacgacctctctgaacgtctgcgtctgcgcggctggcaggttccggccttc     <  1:2855651‑M1/100‑97 (MQ=255)                                                                                               ggCGgtttgcttcaaactgaaagatggtgaagatccgggatacaccctgtacgacctctctgaacgtctgcgtctgcgcggctggcaggttccggccttc     <  2:586112‑M1/100‑97 (MQ=255)                                                                                               ggCGgtttgcttcaaactgaaagatggtgaagatccgggatacaccctgtacgacctctctgaacgtctgcgtctgcgcggctggcaggttccggccttc     >  1:1574757‑M1/1‑4 (MQ=255)                                                                                                 cGgtttgcttcaaactgaaagatggtgaagatccgggatacaccctgtacgacctctctgaacgtctgcgtctgcgcggctggcaggttccggccttcac   <  2:1732347‑M1/100‑99 (MQ=255)                                                                                                  ggtttgcttcaaactgaaagatggtgaagatccgggatacaccctgtacgacctctctgaacgtctgcgtctgcgcggctggcaggttccggccttcact  <  2:1605440‑M1/100‑100 (MQ=255)                                                                                                  ggtttgcttcaaactgaaagatggtgaagatccgggatacaccctgtacgacctctctgaacgtctgcgtctgcgcggctggcaggttccggccttcact  <  2:2527840‑M1/100‑100 (MQ=255)                                                                                                  |                                                                                                    CGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAACCGCCGGGATGCCTTCGTCCGGGCGACCCGTACAGATGAACTCATACGGCCCCAGTTTGGCGATTTCATCCGCCAGATAAGCGGCAACCTGGTAAGAGGCGTTCTGTACTTTGGTATAGCCTTC  >  NC_000913/1570841‑1571037

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