New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3919777 =38 (1.500)3 (0.130) 3/240 2.8 8.7% coding (81/1542 nt) atpA F1 sector of membrane‑bound ATP synthase, alpha subunit
?NC_000913 3919798 = 27 (1.140)coding (60/1542 nt) atpA F1 sector of membrane‑bound ATP synthase, alpha subunit
Rejected: Frequency below/above cutoff threshold.

GTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGACGTCTGCAAGGCGCTCAAGACGACCGCGTACGCTGCCATCAATGACCA  >  NC_000913/3919651‑3919931
                                                                                                                                                   |                                                                                                                                     
ggcgcgCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACAt                                                                                                                                      <  6:73101/147‑1 (MQ=255)
gatggcagcgtacgcggtcgtcttgagcgccttgcagacgtcttgcagtcttaaggggactggagcatgcaactgaattccaccgaaatcagcgaactgatcaagcagcgcattgctcagttcaatGTTGTGAGCTTCACTCacaaca                                                                                                                                       >  3:21481‑M1/127‑148 (MQ=255)
    cgcTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGaa                                                                                                                                  <  6:16163/149‑1 (MQ=255)
     gcTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAAc                                                                                                                                 >  7:104769/1‑149 (MQ=255)
         gAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGaa                                                                                                                                  <  6:14711/144‑1 (MQ=255)
         gAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGaa                                                                                                                                  >  5:14711/1‑144 (MQ=255)
           ggTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCa                                                                                                                           <  2:27725/149‑1 (MQ=255)
                  tGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCgct                                                                                                                    >  1:56792/1‑149 (MQ=255)
                    cGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCgctgc                                                                                                                  <  2:56792/149‑1 (MQ=255)
                     gATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCgctgct                                                                                                                 >  3:32138/1‑149 (MQ=255)
                                          cAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCgg                                                                                            <  4:32138/149‑1 (MQ=255)
                                                 aTCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAAt                                                                                      <  8:104769/148‑1 (MQ=255)
                                                                  aCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGAt                                                                                                             >  3:67693/1‑108 (MQ=255)
                                                                  aCAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGAt                                                                                                             <  4:67693/108‑1 (MQ=255)
                                                                          ccAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCtt                                                              >  8:9274/1‑147 (MQ=255)
                                                                            aGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGa                                                                                                              >  2:120483/1‑97 (MQ=255)
                                                                            aGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGa                                                                                                              <  1:120483/97‑1 (MQ=255)
                                                                            aGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAAtt                                                                                     >  8:2517/1‑122 (MQ=255)
                                                                            aGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAAtt                                                                                     <  7:2517/122‑1 (MQ=255)
                                                                                caccgaaatcagcgaactgatcaagcagcgcattgctcagttcaatGTTGTGAGCTTCACTCACAACATTGaa                                                                                                                                  <  8:25693‑M1/27‑1 (MQ=255)
                                                                                caccgaaatcagcgaactgatcaagcagcgcattgctcagttcaatGTTGTGAGCTTCACTCACAACATTGaa                                                                                                                                  >  7:25693‑M1/47‑73 (MQ=255)
                                                                                 gTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGc                                                      >  5:94493/1‑148 (MQ=255)
                                                                                    aaTGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGa                                                  <  2:10875/149‑1 (MQ=255)
                                                                                           aTAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGACGTCTGc                                           <  6:44975/149‑1 (MQ=255)
                                                                                                  cGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGACGTCTGCAAGgcgc                                    >  2:21063/1‑149 (MQ=255)
                                                                                                        ttACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGACGTCTGCAAGGCGCTCAAga                              >  8:82820/1‑149 (MQ=255)
                                                                                                              aaaCAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTc                                                                        <  7:5986/101‑1 (MQ=255)
                                                                                                              aaaCAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTc                                                                        >  8:5986/1‑101 (MQ=255)
                                                                                                                      gTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGACGTCTGCAAGGCGCTCAAGACGACCGCGTACGCt                <  6:94493/149‑1 (MQ=255)
                                                                                                                      gTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGACGTCTGCAAGGCGCTCAAGACGACCGCGTACGCt                >  6:71417/1‑149 (MQ=255)
                                                                                                                        ttcaatgttgtgagtgaagctcacaacATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTg                                                                                          <  8:25693‑M2/46‑1 (MQ=255)
                                                                                                                        ttcaatgttgtgagtgaagctcacaacATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTg                                                                                          >  7:25693‑M2/28‑73 (MQ=255)
                                                                                                                             tgttgtgagtgaagctcacaacATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGACGTCTGCAAGGCGCTCAAGACGACCGCGTACGCTGCCATc          >  3:21481‑M2/23‑148 (MQ=255)
                                                                                                                                    AGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGACGTCTGCAAGGCGCTCAAGACGACCGCGTACGCTGCCATCAATGACCa  >  5:71417‑M1/1‑149 (MQ=255)
                                                                                                                                            tcacaacATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGACGTCTGCAAGGCGCTCAAGACGACCGCGTACGCTGCCATCAAt       >  4:32465‑M2/8‑136 (MQ=255)
                                                                                                                                            tcacaacATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGACGTCTGCAAGGCGCTCAAGACGACCGCGTACGCTGCCATCAAt       <  3:32465‑M2/129‑1 (MQ=255)
                                                                                                                                               caacATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTc                                                                        >  3:29477‑M2/5‑68 (MQ=255)
                                                                                                                                               caacATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTc                                                                        <  4:29477‑M2/64‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                     
GTCGCGCTCGAGGTTCAGTGCGATAGCGTAACGGTTACCCGGCAGGGAGATCATTTCACCCTGCATACAATCGGCCAGGCCGTGAATGCGGATAACACCGTCACTTACAGAAACAATAGTACCTTCGTTGTGAGCTTCACTCACAACATTGAACTGAGCAATGCGCTGCTTGATCAGTTCGCTGATTTCGGTGGAATTCAGTTGCATGCTCCAGTCCCCTTAAGACTGCAAGACGTCTGCAAGGCGCTCAAGACGACCGCGTACGCTGCCATCAATGACCA  >  NC_000913/3919651‑3919931

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