New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA = 403879831 (0.790)4 (0.110) 4/268 4.7 12.4% noncoding (1280/2905 nt) rrlA 23S ribosomal RNA of rrnA operon
?NC_000913_3_pae_tpiA = 4038805 27 (0.720)noncoding (1287/2905 nt) rrlA 23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAG  >  NC_000913_3_pae_tpiA/4038659‑4038872
                                                                                                                                                  |                                                                   
tAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGcc                                                                   <  2:95606/149‑1 (MQ=255)
tAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGcc                                                                   >  1:459761/1‑149 (MQ=255)
       tGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGc                                                            >  1:441453/1‑149 (MQ=255)
           ccGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGa                                                        >  1:432749/1‑149 (MQ=255)
                        aGCGACACTAGGAGGTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTc                                           <  2:38929/149‑1 (MQ=255)
                            acacTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGGCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGt                                       >  1:691360/1‑149 (MQ=17)
                             aacttagtgtTGTTGGGTAGGGGAGGGTTCGTTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATACGTAACGACAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAATTGTTCCTGTc                                      <  2:588553/143‑1 (MQ=11)
                               cTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGGGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATCAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGATCCttcgca                                    >  2:29888/1‑144 (MQ=255)
                                      ttgttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGACATACTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAAACCAAGGGTTCCGGTCCAACGTTaa                             >  2:466412/1‑149 (MQ=255)
                                        gttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTAAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                       <  1:364537/121‑1 (MQ=17)
                                        gttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTAAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                       >  2:364537/1‑121 (MQ=18)
                                          tgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCgg                         <  2:441453/149‑1 (MQ=32)
                                          tgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAAGGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGCCCAAGGGTTCCTGTCCAACGTTAATCgg                         <  1:699326/149‑1 (MQ=17)
                                           gggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCggg                        <  1:271484/149‑1 (MQ=32)
                                           gggCAGTGGAGCGTTGTGTATGCATGTGAAGGTGCGCTGCGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAACTAACGATAAAGCGGGTGAAAAGCGCGCTCGCCGGAAGACCAGGGTTTCCGGGCCAACGTTGATCggg                        <  2:574097/149‑1 (MQ=11)
                                            ggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCgggg                       >  2:47818/1‑149 (MQ=32)
                                            ggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCgggg                       >  2:11411/1‑149 (MQ=32)
                                              tAGGGGAGCGTTCCGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAGGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGAGCAAGGGTTCCTGTCCAACGTTAATCGGGGCa                     >  2:280405/1‑149 (MQ=12)
                                               aGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAg                    <  1:116432/149‑1 (MQ=32)
                                                ggggAGCGTTCTGTAAGCCTGTGAAGGTGTGCAGTGAGGCATACTGGAGGCATCAGAAGTGCGAACGCTGACATAGGTAACGATAAAGCGGGTGAAAATCCCGCTAGCCGGAAGACCGAGGGTTCCTGTCCAACGTTAATCGGGTCAgg                   >  2:662609/1‑149 (MQ=11)
                                                  ggAGTGTTCTGTAAGGCTGTGAAGGGTTGCGGTGAGGCATGTGGGAGGTAGCAGAAGTGCTAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGt                 <  2:487298/149‑1 (MQ=11)
                                                  ggAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGt                 <  2:579757/149‑1 (MQ=32)
                                                   gAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTg                >  1:228268/1‑149 (MQ=32)
                                                   gAGCGTTCTGTAAGCCTGTGAAGGTGTGCCGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTg                >  1:592555/1‑149 (MQ=17)
                                                         tCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGa          <  1:47699/149‑1 (MQ=32)
                                                           tGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGACGACGCAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAcc        >  1:306674/1‑149 (MQ=12)
                                                           tGTAAGCCTGTGAAGGTGTGCTGTGAGCCCTGCGTCGGTTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAcc        >  2:502790/1‑149 (MQ=11)
                                                            gTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAccc       <  2:459761/149‑1 (MQ=32)
                                                                 ccTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAg  <  1:11411/149‑1 (MQ=32)
                                                                                                                                                  |                                                                   
TAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAG  >  NC_000913_3_pae_tpiA/4038659‑4038872

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