New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 63968338 (1.190)5 (0.160) 5/252 3.0 12.0% intergenic (+175/‑70) ahpC/ahpF alkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
?NC_000913 = 639725 36 (1.150)intergenic (+217/‑28) ahpC/ahpF alkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Rejected: Frequency below/above cutoff threshold.

ACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGG  >  NC_000913/639592‑639864
                                                                                                                                     |                                                                                                                                           
acaGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGc                                                                                                                               >  4:87428/1‑148 (MQ=255)
acaGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGc                                                                                                                               <  3:87428/148‑1 (MQ=255)
         ccgccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGAt                                                                                                                      >  8:119959/1‑148 (MQ=255)
           gccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATaaa                                                                                                                   <  2:100341/149‑1 (MQ=255)
                  tctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGtt                                                                                                                                     >  4:41218/1‑124 (MQ=255)
                  tctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGtt                                                                                                                                     <  3:41218/124‑1 (MQ=255)
                   ctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTc                                                                                                            >  5:22128/1‑148 (MQ=255)
                      aCGCGAAATACTTCGGAAATTCACCTAATTCTTCGGTTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCGAGGAGATAAACATGCTCGaca                                                                                                        <  4:130095/149‑1 (MQ=255)
                                 ttCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGaa                                                                                              >  6:56579/1‑148 (MQ=255)
                                   cGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGaaaa                                                                                            <  6:98994/148‑1 (MQ=255)
                                     gAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactc                                                                                         <  5:56579/149‑1 (MQ=255)
                                        aTTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactcaac                                                                                      <  6:22128/149‑1 (MQ=255)
                                        aTTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactcaac                                                                                      >  4:22481/1‑149 (MQ=255)
                                              cTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGt                                                                                                                                      >  7:41973/1‑95 (MQ=255)
                                              cTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGt                                                                                                                                      <  8:41973/95‑1 (MQ=255)
                                              cTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAgg                                                                                <  3:22481/149‑1 (MQ=255)
                                                aaTTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCt                                                                              >  6:84807/1‑149 (MQ=255)
                                                   aaatcGGGTGCCGCAGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTAc                                                                           >  4:67096/4‑149 (MQ=255)
                                                     atcGGGTGCCGCAGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCt                                                                         >  4:131701/2‑149 (MQ=255)
                                                        gggTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCGTACCTTga                                                                      <  8:66385/149‑1 (MQ=255)
                                                           tgctgcGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGaa                                                                   <  7:119959/149‑1 (MQ=255)
                                                                cagCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGa                                                              <  3:67096/147‑1 (MQ=255)
                                                                         aTTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTg                                                     >  7:123717/1‑149 (MQ=255)
                                                                               ttcCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTg                                                     >  8:43618/1‑143 (MQ=255)
                                                                               ttcCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTg                                                     <  7:43618/143‑1 (MQ=255)
                                                                               ttcCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGt                                               >  1:53945/1‑149 (MQ=255)
                                                                               ttcCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGt                                               >  2:53549/1‑149 (MQ=255)
                                                                                   ccGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGatggtgcggggaagaaaatcgggtgccgcagcacccgaagaattaggtgaatttccgaagtatttcgcgtcagagcaagg                                                             >  5:146160‑M2/1‑51 (MQ=255)
                                                                                   ccGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGatggtgcggggaagaaaatcgggtgccgcagcacccgaagaattaggtgaatttccgaagtatttcgcgtcagagcaagg                                                             <  6:146160‑M2/131‑81 (MQ=255)
                                                                                     gCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGatggtgcggggaagaaaatcgggtgccgcagcacccgaagaattaggtgaatttccgaagtatttcgcgtcagagcaaggcggctagcctgtgaatcccc                                         <  7:84752‑M2/149‑101 (MQ=255)
                                                                                         cATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATCatggtgcggggaagaaaatcgggtgccgcagcacccgaagaattaggtgaatttccgaagtatttcgcgtcagagcaaggcggctagcctgtgaatccccgg                                       <  7:11089‑M2/147‑103 (MQ=255)
                                                                                              tgCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTg                                 >  5:99704/1‑148 (MQ=255)
                                                                                                cAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGAt                              <  5:84807/149‑1 (MQ=255)
                                                                                                       cATCCAGGTAGCCGCAGAGGCCGCTTGCATCatggtgcggggaagaaaatcgggtgccgcagcacccgaagaattaggtgaatttccgaagtatttcgcgtcagagcaaggcggctagcctgtgaatccccgggagcttacacaagtaa                       >  8:11089‑M2/1‑31 (MQ=255)
                                                                                                                  ccGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCa            <  2:53945/149‑1 (MQ=255)
                                                                                                                   cGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCgg              >  3:52858/1‑146 (MQ=255)
                                                                                                                        agGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATc      <  6:99704/149‑1 (MQ=255)
                                                                                                                          gCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCaa    <  1:53549/149‑1 (MQ=255)
                                                                                                                            cGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAgg  <  8:123717/149‑1 (MQ=255)
                                                                                                                             gCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGtt                                                   >  2:137917/1‑99 (MQ=255)
                                                                                                                             gCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGtt                                                   <  1:137917/99‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                           
ACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGG  >  NC_000913/639592‑639864

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