New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3363001 =	27 (1.140)	34 (1.450)	23/274	0.2	71.7%	coding (195/2382 nt)	yhcD	putative fimbrial usher protein YhcD
	?	NC_000913	3766087 =	0 (0.000)					coding (3905/4134 nt)	rhsA	rhs element protein RhsA

TTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAAT  >  NC_000913/3765948‑3766222                                                                                                                                            |                                                                                                                                        caatatatagtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactt                                                                                                                                         <  1:250397‑M2/1‑1 (MQ=255)        ttCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaag                                                                                                                                  <  2:306407‑M2/140‑6 (MQ=255)        ttCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaag                                                                                                                                  <  1:18476‑M2/140‑6 (MQ=255)         agtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATaa                                                                                                                                 >  2:337022‑M2/132‑140 (MQ=255)         agtcgtaacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATaa                                                                                                                                 >  1:181305‑M2/132‑140 (MQ=255)              taacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATg                                                                                                                            <  2:181305‑M2/14‑1 (MQ=255)              taacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATg                                                                                                                            >  2:208869‑M2/127‑140 (MQ=255)              taacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATg                                                                                                                            >  1:180772‑M2/127‑140 (MQ=255)              taacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATg                                                                                                                            >  1:301219‑M2/127‑140 (MQ=255)              taacagccatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATg                                                                                                                            >  1:24748‑M2/127‑140 (MQ=255)                 ccAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaaca                                                                                                                         >  1:306407‑M2/1‑126 (MQ=255)                   aaaaGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaag                                                                                                                       <  2:254780‑M2/140‑17 (MQ=255)                    acatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAg                                                                                                                                       >  2:39703‑M2/121‑123 (MQ=255)                    acatttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAg                                                                                                                                       <  1:39703‑M2/3‑1 (MQ=255)                     catttttttcttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCtt                                                                                                                     >  2:283902‑M2/120‑140 (MQ=255)                      aGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtcc                                                                                                                    >  2:264494‑M2/1‑121 (MQ=255)                      aGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtcc                                                                                                                    >  2:337122‑M2/1‑121 (MQ=255)                      aGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtcc                                                                                                                    >  1:394678‑M2/1‑121 (MQ=255)                       gCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccc                                                                                                                   >  2:257630‑M2/1‑120 (MQ=255)                              cttcagaaggaatttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGaa                                                                                                            <  2:175001‑M2/30‑1 (MQ=255)                                           ttttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGacac                                                                                               <  2:24748‑M2/43‑1 (MQ=255)                                             ttaaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAg                                                                                             <  1:99338‑M2/45‑1 (MQ=255)                                               aaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAggg                                                                                           <  2:301219‑M2/47‑1 (MQ=255)                                               aaaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAggg                                                                                           <  2:180772‑M2/47‑1 (MQ=255)                                                aaccggcttgttgcagaaatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGa                                                                                          >  1:20568‑M2/93‑140 (MQ=255)                                                   tgctgcAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                                                                       <  2:156825‑M2/140‑49 (MQ=255)                                                    gctgcAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggag                                                                                      >  2:18476‑M2/1‑91 (MQ=255)                                                              cATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaaca                                                                            >  1:139550‑M2/1‑81 (MQ=255)                                                                 tGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacaggg                                                                         >  2:129375‑M2/1‑78 (MQ=255)                                                                 aatgcctgatcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAAt                                                                         >  1:124810‑M2/76‑140 (MQ=255)                                                                          tcaagacaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACataata                                                                >  1:344600‑M2/67‑140 (MQ=255)                                                                               acaaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGa                                                           >  2:373142‑M2/62‑140 (MQ=255)                                                                                caaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGAc                                                          <  1:208869‑M2/80‑1 (MQ=255)                                                                                caaatgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGAc                                                          <  2:20568‑M2/80‑1 (MQ=255)                                                                                   atgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGAc                                                          <  1:222568‑M2/80‑1 (MQ=255)                                                                                   atgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGAc                                                          >  2:222568‑M2/58‑137 (MQ=255)                                                                                   atgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGa                                                       >  2:120380‑M2/58‑140 (MQ=255)                                                                                   atgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGa                                                       <  1:337022‑M2/83‑1 (MQ=255)                                                                                         ccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTgg                                                 >  2:135531‑M2/52‑140 (MQ=255)                                                                                            cacaAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc                                              <  2:344840‑M2/140‑90 (MQ=255)                                                                                              caAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgc                                            <  2:145004‑M2/140‑92 (MQ=255)                                                                                                   caaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGa                                       <  2:344600‑M2/99‑1 (MQ=255)                                                                                                             tctctAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcag                             <  1:359491‑M2/140‑107 (MQ=255)                                                                                                                   gttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGAtg                       <  1:135531‑M2/115‑1 (MQ=255)                                                                                                                     atTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagc                     >  1:184619‑M2/1‑26 (MQ=255)                                                                                                                       ccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGAc                   >  1:77166‑M2/22‑140 (MQ=255)                                                                                                                       ccttctttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGAc                   >  1:148260‑M2/22‑140 (MQ=255)                                                                                                                         aataatCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccg                 <  1:365113‑M2/140‑119 (MQ=255)                                                                                                                         aataatCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccg                 <  1:206811‑M2/140‑119 (MQ=255)                                                                                                                         aataatCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccg                 <  1:193355‑M2/140‑119 (MQ=255)                                                                                                                         aataatCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccg                 <  1:89786‑M2/140‑119 (MQ=255)                                                                                                                               ttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCa           <  1:104602‑M2/127‑1 (MQ=255)                                                                                                                               ccGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagac           >  2:359491‑M2/1‑16 (MQ=255)                                                                                                                                 cgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAgg         <  1:283902‑M2/129‑1 (MQ=255)                                                                                                                                 cgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAgg         <  2:338840‑M2/129‑1 (MQ=255)                                                                                                                                    aCTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaa      >  2:193355‑M2/1‑11 (MQ=255)                                                                                                                                        gCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaaga  <  2:304386‑M2/140‑134 (MQ=255)                                                                                                                                            |                                                                                                                                        TTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAAT  >  NC_000913/3765948‑3766222

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