New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2941507	28 (0.710)	3 (0.080)	3/258	3.5	9.5%	coding (125/396 nt)	ygdD	conserved inner membrane protein YgdD
	?	NC_000913	= 2941815	31 (0.840)					coding (753/918 nt)	gcvA	DNA‑binding transcriptional dual regulator GcvA
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AAAGGCTGCCGCTGAACAACACCGTGCCTAACGCGAGGAAAACGCTACTCCAGTAAAACCAGATACTGATGCGACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATAAAACGGCTGGTCATGGTCGTACCCTAC  >  NC_000913/2941370‑2941644                                                                                                                                          |                                                                                                                                          aaaGGCTGCCGCTGAACAACACCGTGCCTAACGCGAGGAAAACGCTACTCCAGTAAAACCAGATACTGATGCGACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTc                                                                                                                                         <  1:416637/140‑1 (MQ=255)       tGCCGCTGAACAACACCGTGCCTAACGCGAGGAAAACGCTACTCCAGTAAAACCAGATACTGATGCGACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGTATc                                                                                                                                   >  1:163148/1‑140 (MQ=255)                 caacaCCGTGCCTAACGCGAGGAAAACGCTACTCCAGTAAAACCAGATACTGATGCGACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCAtct                                                                                                                         <  1:662292/140‑1 (MQ=255)                         tGCCTAACGCGAGGAAAACGCTACTCCAGTAAAACCAGATACTGATGCGACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGcc                                                                                                                 >  2:421486/1‑140 (MQ=255)                            cTAACGCGAGGAAAACGCTACTCCAGTAAAACCAGATACTGATGCGACGCTGCGTTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACACCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGccccc                                                                                                              <  2:333703/140‑1 (MQ=255)                            cTAACGCGAGGAAAACGCTACTCCAGTAAAACCAGATACTGATGCGACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGccccc                                                                                                              >  1:625901/1‑140 (MQ=255)                             tAACGCGAGGAAAACGCTACTCCAGTAAAACCAGATACTGATGCGACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCa                                                                                                             <  1:421486/140‑1 (MQ=255)                                  cgAGGAAAACGCTCCTCTAGTAAAACCAGATACTGATGCGACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGtg                                                                                                        <  2:478007/140‑2 (MQ=255)                                  cgAGGAAAACGCTACTCCAGTAAAACCAGATACTGATGCGACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGtt                                                                                                        <  2:94380/140‑1 (MQ=255)                                  cgAGGAAAACGCTACTCCAGTAAAACCAGATACTGATGCGACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGtt                                                                                                        <  1:440246/140‑1 (MQ=255)                                  cgAGGAAAACGCTACTCAAGTACAACCAGATTCTGATGCTACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGtt                                                                                                        <  2:199648/140‑1 (MQ=255)                                                                                 aTTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGaa                                                         <  2:625901/140‑1 (MQ=255)                                                                                          gCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAAt                                                >  2:335035/1‑140 (MQ=255)                                                                                                      aTCGCCAGCGTATGAAACGCCTTGTATTCGATGCCCGTCGGGATCCAGTCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAAt                                    <  2:239547/140‑1 (MQ=255)                                                                                                          ccAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGc                                >  2:646632/1‑140 (MQ=255)                                                                                                              cGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATaa                            >  1:638512/1‑140 (MQ=255)                                                                                                              cGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATaa                            >  2:530541/1‑140 (MQ=255)                                                                                                                tATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCAAGAGCCACAAAAGTGAAGCCGCTAATGGCGCCGAAAATCAGCATaaaa                          >  2:233077/1‑140 (MQ=255)                                                                                                                   gAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATAAAACgg                       <  1:344623/140‑1 (MQ=255)                                                                                                                   gAAACGCCTGGCATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATAAAACgg                       <  2:163148/140‑1 (MQ=255)                                                                                                                        gCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATAAAACGGCTGGt                  <  1:530541/140‑1 (MQ=255)                                                                                                                          cTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATAAAACGGCTGGTCa                <  1:335035/140‑1 (MQ=255)                                                                                                                          cTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATAAAACGGCTGGTCa                <  1:493397/140‑1 (MQ=255)                                                                                                                             gTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATAAAACGGCTGGTCATgg             >  1:383077/1‑140 (MQ=255)                                                                                                                                  cGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATAAAACGGCTGGTCATGGTCGTa        >  2:448319/1‑140 (MQ=255)                                                                                                                                   gAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGGCGCTAATGGCGGAGGAGATCAGCATAAAACTGCTGGTCAGGGTCggtc       >  2:312546/1‑137 (MQ=255)                                                                                                                                     ggCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGAGAAAATCTGCAAAACACGGCTGGCCATGGTCGTAccc     >  1:268448/1‑140 (MQ=255)                                                                                                                                        cGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATAAAACGGCTGGTCATGGTCGTACCCTAc  <  2:383077/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                          AAAGGCTGCCGCTGAACAACACCGTGCCTAACGCGAGGAAAACGCTACTCCAGTAAAACCAGATACTGATGCGACGCTGCATTGCCACCGCCAGACCTAAGATCGCCAGCGTATGAAACGCCTGGTATTCGAGGCCGGTCTGGATCCAGCCCATCTCAACGGCCCCCATGGTTTTACTTAACACATGCGCGCCAAAAGCGCCCAGAGCCACAAAAATGAAGCCGCTAATGGCGGCGAAAATCAGCATAAAACGGCTGGTCATGGTCGTACCCTAC  >  NC_000913/2941370‑2941644

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