New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA = 3944991NA (NA)4 (0.110) 4/268 4.7 11.9% noncoding (1288/2904 nt) rrlC 23S ribosomal RNA of rrnC operon
?NC_000913_3_pae_tpiA = 4038798 31 (0.790)noncoding (1280/2905 nt) rrlA 23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAG  >  NC_000913_3_pae_tpiA/4038651‑4038872
                                                                                                                                                   |                                                                          
aaCGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGt                                                                           <  2:469085/149‑1 (MQ=255)
     ggCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGaaaa                                                                      <  2:53280/149‑1 (MQ=40)
        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                                                       >  2:364537/1‑121 (MQ=18)
                                                gttgGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTAAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                       <  1:364537/121‑1 (MQ=17)
                                                  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:11411/1‑149 (MQ=32)
                                                    ggTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCgggg                       >  2:47818/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)
                                                                                                                                                   |                                                                          
AACGGGGCTAAACCATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAG  >  NC_000913_3_pae_tpiA/4038651‑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.