New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 335704333 (0.980)4 (0.130) 4/264 4.1 11.1% coding (2319/4461 nt) gltB glutamate synthase, large subunit
?NC_000913_3_bme_pgi = 3357057 33 (1.040)coding (2333/4461 nt) gltB glutamate synthase, large subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCGCTGCTCGAAACTGTTTGAAGCGGTCGGTCTACACGATGATGTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGACTATCAGGAATACGCGAAGCTGGT  >  NC_000913_3_bme_pgi/3356899‑3357189
                                                                                                                                                |                                                                                                                                                  
ccGCTGCTCGAAACTGTTTGAAGCGGTCGGTCTACACGATGATGTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGgcgc                                                                                                                                                >  2:168018/1‑149 (MQ=255)
  gctgctCGAAACTGTTTGAAGCGGTCGGTCTACACGATGATGTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGGGGAGCAAGCTTTGAATACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGt                                                                                                                                              >  1:48395/1‑149 (MQ=255)
                tttGAAGCGGTCGGTCTACACGATGATGTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCa                                                                                                                                <  1:549625/149‑1 (MQ=255)
                   gAAGCGGTCGGTCTACACGATGATGTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCACCCAggg                                                                                                                             >  1:421388/1‑149 (MQ=255)
                     aGCGGTCGGTCTACACGATGATGTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCg                                                                                                                           >  1:72400/1‑149 (MQ=255)
                             gtcTACACGATGATGTAGTGGGCCTGTGCTTCCAGGTGGCGGGCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTctgctg                                                                                                                   <  1:8695/149‑1 (MQ=255)
                              tcTACACGATGATGTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTCCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGa                                                                                                                  >  1:313637/1‑149 (MQ=255)
                                tACACGATGATGTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTTCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGaaa                                                                                                                <  1:168018/149‑1 (MQ=255)
                                       tgatgTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTcc                                                                                                         <  2:72400/149‑1 (MQ=255)
                                         atgTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCAc                                                                                                       >  2:461726/1‑149 (MQ=255)
                                          tgTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGCCTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGACCATCAACCAGGGCGGTCTGCTGAAATACGTCCAcg                                                                                                      >  1:58287/1‑149 (MQ=255)
                                                gggCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGTGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGgcggcg                                                                                                <  2:590209/149‑1 (MQ=255)
                                                        gCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGGGAATAcccc                                                                                        >  1:128665/1‑147 (MQ=255)
                                                               gggggCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACa                                                                                 >  1:363689/1‑149 (MQ=255)
                                                                       tCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGAGAGCTTGGCAATTACGTCCACGGCTGCGAATACCAAGCCTACAACCCAGAc                                                                         >  1:495207/1‑149 (MQ=255)
                                                                           ccGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATAAGTCCAAGGGGGCGGATACAACGCCTACAACTCGTATgtgg                                                                     >  1:87151/1‑149 (MQ=255)
                                                                            cGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACgtggt                                                                    <  1:461726/149‑1 (MQ=255)
                                                                                     ggAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCt                                                           <  2:363689/149‑1 (MQ=255)
                                                                                       aGCAAGCTTTGAAGCCTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCTGc                                                         <  2:48395/149‑1 (MQ=255)
                                                                                        gCAAGCTTTGAAGACTTGCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTTAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCTGca                                                        <  2:288025/149‑1 (MQ=255)
                                                                                         cAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCTGcaa                                                       >  1:325833/1‑149 (MQ=255)
                                                                                                 gAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCTGCAACAAGCGGt                                               >  2:484064/1‑149 (MQ=255)
                                                                                                   aGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCTGCAACAAGCGGTAc                                             >  1:170877/1‑149 (MQ=255)
                                                                                                       ttCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGc                                                                                      <  1:146945/104‑1 (MQ=255)
                                                                                                       ttCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGc                                                                                      >  2:146945/1‑104 (MQ=255)
                                                                                                             cagGATCTGCTCAATCTGTCTAAACGTGCCTGGCTGGCGCGTAAGTCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGa                                   <  1:94402/149‑1 (MQ=255)
                                                                                                                      ctgAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTTGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGa                          <  1:89091/149‑1 (MQ=255)
                                                                                                                       tgAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGAc                         <  2:170877/149‑1 (MQ=255)
                                                                                                                            cTGTCGAAACGTGCGTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGTCGTGGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGACTATCa                    <  2:58287/149‑1 (MQ=255)
                                                                                                                             tGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGTCGGCGAATACCACGCCTACAACCCGGACGAGGTACGCACGCTGCAAAAAGCGGTACAAAACGGCGAGTACAGCGACTATCAg                   >  1:8437/1‑149 (MQ=255)
                                                                                                                               tCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGACTATCAGGa                 >  2:339016/1‑149 (MQ=255)
                                                                                                                                         ccTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCTGCAAAAAGCGGTACAAAGCGGCGAGTACAGAGACTAACAGGAATACGCGAAg       >  1:197505/1‑149 (MQ=255)
                                                                                                                                              ctggaGCGTAAGCTACGCAGCCAGTGCGGGCTGCGCGAATGCGTCCCCGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACTCTTCGACAAGCGGTACAAAGCGGCGAGTACATCGACTATCAGGAATACGCGAAGCTGGt  <  2:83773/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCGCTGCTCGAAACTGTTTGAAGCGGTCGGTCTACACGATGATGTAGTGGGCCTGTGCTTCCAGGGGGCGGTCAGCCGCATTGGTGGAGCAAGCTTTGAAGACTTCCAGCAGGATCTGCTGAATCTGTCGAAACGTGCCTGGCTGGCGCGTAAGCCCATCAGCCAGGGCGGTCTGCTGAAATACGTCCACGGCGGCGAATACCACGCCTACAACCCGGACGTGGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGACTATCAGGAATACGCGAAGCTGGT  >  NC_000913_3_bme_pgi/3356899‑3357189

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