New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 4001760 =34 (0.850)3 (0.080) 3/300 NT 9.6% intergenic (+77/‑57) dmoA‑II/allE dimethyl‑sulfide monooxygenase/putative S‑ureidoglycine aminohydrolase
?AE015451 4001778 = 26 (0.730)intergenic (+95/‑39) dmoA‑II/allE dimethyl‑sulfide monooxygenase/putative S‑ureidoglycine aminohydrolase

ACCTGCTGGTGCCCGAGCTGCAACGGCGCGGCCGCTACCCGCTGAGCTACACCGAGGGCCCCCTGCGCCACAAGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGG  >  AE015451/4001547‑4001981
                                                                                                                                                                                                                                       |                                                                                                                                                                                                           
aCCTGCTGGTGCCCGAGCTGCAACGGCGCGGCCGCTACCCGCTGAGCTACACCGAGGGCCCCCTGCGCCACAAGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGtt                                                                                                                                                                                                            >  2:321251/1‑233 (MQ=255)
                 cTGCAACGGCGCGGCCGCTACCCGCTGAGCTACACCGAGGGCCCCCTGCGCCACAAGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATgggg                                                                                                                                                                                                <  1:229380/228‑1 (MQ=255)
                                 gCTACCCGCTGAGCTACACCGAGGGCCCCCTGCGCCACAAGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTc                                                                                                                                                                        >  1:742031/1‑236 (MQ=255)
                                 gCTACCCGCTGAGATACACCGAGGGCCCCCTGCGCCACAAGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTc                                                                                                                                                                                 <  1:93169/227‑1 (MQ=255)
                                                             ccTGCGCCACAAGCTGGTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAacaaca                                                                                                                                                                                     <  1:125875/195‑1 (MQ=255)
                                                              cTGCGCCACAAGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGtt                                                                                                                                                                                                            <  2:63189/171‑1 (MQ=255)
                                                              cTGCGCCACAAGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGtt                                                                                                                                                                                                            >  1:63189/1‑171 (MQ=255)
                                                               tGCGCCACAAGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCa                                                                                                                                          >  1:64544/1‑236 (MQ=255)
                                                                      caAGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTggg                                                                                                                                   <  1:197093/236‑1 (MQ=255)
                                                                       aaGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTggg                                                                                                                                   <  1:321252/235‑1 (MQ=255)
                                                                       aaGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGGTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGc                                                                                                                                  <  1:429849/236‑1 (MQ=255)
                                                                                ggCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCaa                                                                                                                                                                                                         <  2:534762/156‑1 (MQ=255)
                                                                                                                   aCGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACCTTTAATTTGTGGCTGCAACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGga                                                                                                                                                                                               <  2:571230/131‑1 (MQ=255)
                                                                                                                        aggtcagCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGa                                                                                    >  2:497006/1‑233 (MQ=255)
                                                                                                                         tcggttggagccgggtgcccaccgtggggggcgaagtatgaatggttcgacatggacgatctcgacttgttgttctctccccatgccttgaaCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGt                                                                                                                                                                   >  2:202463‑M1/93‑153 (MQ=255)
                                                                                                                             agCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTc                                                                                                                                                         >  2:612504/1‑159 (MQ=255)
                                                                                                                               cATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGcc                                                                          <  1:182866/236‑1 (MQ=255)
                                                                                                                               cATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGcc                                                                          <  1:497007/236‑1 (MQ=255)
                                                                                                                                            aaTACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGCTCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCAATGTTCACCGAAGCCTATGCCGTCATCCCCa                                                                >  1:5913/1‑233 (MQ=255)
                                                                                                                                                   ggTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGCTCGTCCATGTCGAACCATTCATACTTCGCCCCCCACgg                                                                                                                                       >  2:539196/1‑155 (MQ=255)
                                                                                                                                                     tATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCa                                                                                                                                          >  1:267012/1‑150 (MQ=255)
                                                                                                                                                     tATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCa                                                                                                                                          <  2:267011/150‑1 (MQ=255)
                                                                                                                                                     tATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACgg                                                                                                                                       <  1:689867/153‑1 (MQ=255)
                                                                                                                                                     tATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACgg                                                                                                                                       >  2:689866/1‑153 (MQ=255)
                                                                                                                                                      aTTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGt                                                                                                                                                                   >  1:606399/1‑124 (MQ=255)
                                                                                                                                                      aTTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGt                                                                                                                                                                   <  2:606398/124‑1 (MQ=255)
                                                                                                                                                      aTTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCAcc                                                                                                                               >  1:322196/1‑160 (MQ=255)
                                                                                                                                                                ttAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATcc                                                                   <  2:742030/210‑1 (MQ=255)
                                                                                                                                                                gaatggttcgacatggacgttctcgacttgttgttctctccccatgccttgaaCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTc                                                                                                                                                                  >  2:379410‑M1/54‑115 (MQ=255)
                                                                                                                                                                         gacatggacgatctcgacttgttgttctctccccatgccttgaaccgcaccggccggtgcggTTCAAGGCATGGGGAGAGAACAACAAGTCGAGAACGTCCATGTCGAACCATTc                                                                                                                                                         >  2:379410‑M2/63‑115 (MQ=255)
                                                                                                                                                                          acatggacgatctcgacttgttgttctctccccatgccttgaaccgcaccggccggtgcggTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCCAACCGa                                                                                                                  >  2:202463‑M2/62‑153 (MQ=255)
                                                                                                                                                                           aTGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCAcc                               >  1:421014/1‑235 (MQ=255)
                                                                                                                                                                                gcgatctcgacttgttgttctctccccatgccttgaaccgccccggccggtgcggTTCAAGGCATGGGGAGAGAACAACAAGTCg                                                                                                                                                                                <  1:186707‑M2/30‑1 (MQ=255)
                                                                                                                                                                                       cgacttgttgttctctccccatgccttgaaCCGCACCGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATcgc                                                                                                                                                                         <  1:186707‑M1/55‑2 (MQ=255)
                                                                                                                                                                                                       ccACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTgg  >  1:671274/1‑236 (MQ=255)
                                                                                                                                                                                                         acaTCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCCCCGAAGCCTATGCCGTCATCCCCAAAGGCGTGACGCGTGACATCGTCACCAccccc                                >  2:78506/1‑202 (MQ=255)
                                                                                                                                                                                                                   cgccgcACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTCCTGACTGCCCGTGCCag                                                                                             >  2:107046/1‑132 (MQ=255)
                                                                                                                                                                                                                          ccggccggtgcggTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGt                                                                                                  <  1:495428‑M2/108‑1 (MQ=255)
                                                                                                                                                                                                                          CCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGt                                                                                                  <  2:495427‑M1/121‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                           
ACCTGCTGGTGCCCGAGCTGCAACGGCGCGGCCGCTACCCGCTGAGCTACACCGAGGGCCCCCTGCGCCACAAGCTGTTCGGCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCCGTCATCCCCAAAGGCGTGATGCGTGACATCGTCACCAGCCACCTGCCCTTCTGGGACAAGATGCGCATGTGG  >  AE015451/4001547‑4001981

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