New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4040481NA (NA)37 (0.400) 9/166 NT 100% intergenic (+58/‑36) rrlA/rrfA 23S ribosomal RNA of rrnA operon/5S ribosomal RNA of rrnA operon
?NC_000913 = 4040505 0 (0.000)intergenic (+82/‑12) rrlA/rrfA 23S ribosomal RNA of rrnA operon/5S ribosomal RNA of rrnA operon

AATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCG  >  NC_000913/4040401‑4040580
                                                                                |                                                                                                   
aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGT                                                                                  <  2:1030833‑M2/100‑1 (MQ=255)
aaTGAACCGTGAGGCTCAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGGGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGT                                                                                  <  2:1006456‑M2/100‑1 (MQ=255)
 aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTC                                                                                 >  1:48000‑M2/1‑100 (MQ=255)
 aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTC                                                                                 <  1:160993‑M2/100‑1 (MQ=255)
 aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTC                                                                                 <  1:1738324‑M2/100‑1 (MQ=255)
 aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTC                                                                                 >  2:563386‑M2/1‑100 (MQ=255)
  tGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCT                                                                                <  1:38638‑M2/100‑1 (MQ=255)
   gAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTG                                                                               >  1:1618395‑M2/1‑100 (MQ=255)
   gAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTG                                                                               <  1:299949‑M2/100‑1 (MQ=255)
   gAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTG                                                                               <  2:827574‑M2/100‑1 (MQ=255)
   gAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCT                                                                                <  1:2451414‑M2/99‑1 (MQ=255)
    aaCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA                                                                              >  1:1477818‑M2/1‑100 (MQ=255)
     aCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATGTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGAT                                                                             <  2:2292859‑M2/100‑1 (MQ=255)
                                                   agagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatctgctctcatccgccaaaaccgc                               <  2:4919‑M2/100‑47 (MQ=255)
                                                       aagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacagcttcg                           <  1:680531‑M2/100‑51 (MQ=255)
                                                           ttttCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacagcttcggcg                       >  1:2692498‑M2/1‑46 (MQ=255)
                                                           ttttCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacagcttcggcg                       >  1:2695544‑M2/1‑46 (MQ=255)
                                                                aGCCTGATACAGATTAAATCAGACCGCTTCTGCGTTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacaccttcggcgttgta                  >  1:10489‑M2/1‑41 (MQ=255)
                                                                          aGATTAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacatcttcggcgttgtaaggttaagcct        <  1:410369‑M2/100‑70 (MQ=255)
                                                                          aGATTAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacagcttcggcgttgtaaggttaagcct        <  1:1754649‑M2/100‑70 (MQ=255)
                                                                          aGATTAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacaccttcggcgttgtacggttaagcc        <  2:1079054‑M2/100‑70 (MQ=255)
                                                                          aGATTAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacaccttcggcgttgtaaggttaagcc        <  2:886195‑M2/100‑70 (MQ=255)
                                                                          aGATTAAATCAGAACGCAGAAGCGGCCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacagcttcggcgttgtaaggttaagcct        <  2:351412‑M2/100‑70 (MQ=255)
                                                                           gATTAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacatcttcggcgttgtaaggttaagcctg       <  1:1259907‑M2/100‑71 (MQ=255)
                                                                           gATTAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacagcttcggcgttgtaaggttaagcctc       <  1:914334‑M2/100‑71 (MQ=255)
                                                                           gATTAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacagcttcggcgttgtaaggttaagcctc       <  1:1846703‑M2/100‑71 (MQ=255)
                                                                             ttAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacatcttcggcgttgtaaggttaagcctcac     <  1:2111491‑M2/100‑73 (MQ=255)
                                                                             ttAAATCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacaccttcggcgttgtaaggttaagcctca     <  1:562367‑M2/100‑73 (MQ=255)
                                                                               aaaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacatcttcggcgttgtaaggttaagcctcacgg   <  1:2657870‑M2/100‑75 (MQ=255)
                                                                               aaaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacagcttcggcgttgtaaggttaagcctcacgg   <  1:830474‑M2/100‑75 (MQ=255)
                                                                               aaaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacagcttcggcgttgtaaggttaagcctcacg   <  2:1270113‑M2/100‑75 (MQ=255)
                                                                               aaaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacagcttcggcgttgtaaggttaagcctcacg   <  2:2594269‑M2/100‑75 (MQ=255)
                                                                               aaaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacaccttcggcgttgtaaggttaagcctcacg   <  2:1624810‑M2/100‑75 (MQ=255)
                                                                               aaaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctatcatccgccaaaacatcttcggcgttgtaaggttaagcctcacgg   <  1:426162‑M2/100‑75 (MQ=255)
                                                                                aaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacagcttcggcgttgtaaggttaagcctcacggt  <  1:480317‑M2/100‑76 (MQ=255)
                                                                                aaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacagcttcggcgttgtaaggttaagcctcacggt  <  2:2155778‑M2/100‑76 (MQ=255)
                                                                                aaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcatccgccaaaacagcttcggcgttgtaaggttaagcctcacggt  <  1:2590704‑M2/100‑76 (MQ=255)
                                                                                aaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacagcttcggcgttgtaaggttaagcctcacgg  <  1:1133309‑M2/100‑76 (MQ=255)
                                                                                aaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacaccttcggcgttgtaaggttaagcctcacgg  <  2:1106607‑M2/100‑76 (MQ=255)
                                                                                aaTCAGAACGCAGAAGCGGTCTGATttaatctgtatcaggctgaaaatcttctctcaatccgccaaaacaccttcgacgttgtaaggttaagcctca     <  2:2665173‑M2/97‑73 (MQ=255)
                                                                                |                                                                                                   
AATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCG  >  NC_000913/4040401‑4040580

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