New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3653230	NA (NA)	23 (0.970)	17/276	0.4	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	= 3766089	0 (0.000)					coding (3907/4134 nt)	rhsA	rhs element protein RhsA

TTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAAT  >  NC_000913/3765955‑3766222                                                                                                                                       |                                                                                                                                      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                                                       <  1:337022‑M2/83‑1 (MQ=255)                                                                            atgtcgccattttcatcaaagcgcgtagcaatgttgccttttttttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGa                                                       >  2:120380‑M2/58‑140 (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:193355‑M2/140‑119 (MQ=255)                                                                                                                  aataatCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccg                 <  1:206811‑M2/140‑119 (MQ=255)                                                                                                                  aataatCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccg                 <  1:89786‑M2/140‑119 (MQ=255)                                                                                                                  aataatCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccg                 <  1:365113‑M2/140‑119 (MQ=255)                                                                                                                        ttcgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCa           <  1:104602‑M2/127‑1 (MQ=255)                                                                                                                        ccGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagac           >  2:359491‑M2/1‑16 (MQ=255)                                                                                                                          cgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAgg         <  2:338840‑M2/129‑1 (MQ=255)                                                                                                                          cgccatttactTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAgg         <  1:283902‑M2/129‑1 (MQ=255)                                                                                                                             aCTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaa      >  2:193355‑M2/1‑11 (MQ=255)                                                                                                                                 gCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaaga  <  2:304386‑M2/140‑134 (MQ=255)                                                                                                                                       |                                                                                                                                      TTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTGTCAATGACATAATAAATGACAGATTCTGGAGTGATGGGAAAAAACCCGACAGATGTGACGTACTTCAGGAACTAAT  >  NC_000913/3765955‑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.