New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 335560636 (1.040)5 (0.150) 5/280 2.4 12.5% coding (882/4461 nt) gltB glutamate synthase, large subunit
?NC_000913 = 3355642 36 (1.100)coding (918/4461 nt) gltB glutamate synthase, large subunit
Rejected: Frequency below cutoff threshold.

CGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGG  >  NC_000913/3355482‑3355755
                                                                                                                            |                                                                                                                                                     
cgcgtACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTc                                                                                                                              <  2:465656/150‑1 (MQ=255)
      cTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGt                                                                                                                       >  2:485828/1‑151 (MQ=255)
      cTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGt                                                                                                                       >  1:318767/1‑151 (MQ=255)
        tataAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTac                                                                                                                     <  2:369256/151‑1 (MQ=255)
            aaTTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcc                                                                                                                 >  2:188719‑M2/1‑149 (MQ=255)
              ttCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcccg                                                                                                               <  1:548610‑M2/151‑5 (MQ=255)
              ttCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcccg                                                                                                               <  1:113647‑M2/151‑5 (MQ=255)
                     cacCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATatcatc                                                                                                                                           <  1:270655/116‑1 (MQ=255)
                     cacCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATatcatc                                                                                                                                           >  2:270660/1‑116 (MQ=255)
                          cTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCGCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGaac                                                                                                    <  1:154927/150‑1 (MQ=255)
                              tCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAAcc                                                                                                >  2:243155/1‑150 (MQ=255)
                              tCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAAcc                                                                                                >  2:60029/1‑150 (MQ=255)
                                       tGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcccgcctgccagcagcagttccagcatgt                                                                                      >  2:321414‑M2/1‑122 (MQ=255)
                                        gCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGtggtactaatagacgcatggcgcggatgatatccatcccgcctgccagcagcagttccagcatgtt                                                                                     >  2:818‑M1/1‑85 (MQ=255)
                                        gCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTac                                                                                                                     <  1:189142/119‑1 (MQ=255)
                                        gCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTac                                                                                                                     >  2:189146/1‑119 (MQ=255)
                                         cACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGAc                                                                                     >  2:382654/1‑150 (MQ=255)
                                             acgCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcccgcctggcagcagcagttccagcatgttatcc                                                                                 >  2:216953‑M2/1‑116 (MQ=255)
                                             acgCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGa                                                                                >  1:1662/1‑151 (MQ=255)
                                             acgCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGa                                                                                >  2:505438/1‑151 (MQ=255)
                                             acgCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGa                                                                                >  1:464073/1‑151 (MQ=255)
                                             acgCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATACCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGa                                                                                >  1:507370/1‑151 (MQ=255)
                                                ccgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGaa                                                                               >  2:241462/1‑149 (MQ=255)
                                                ccgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGaa                                                                               >  2:496908/1‑149 (MQ=255)
                                                ccgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAAc                                                                              <  1:138644/150‑1 (MQ=255)
                                                ccgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAAc                                                                              <  2:139003/150‑1 (MQ=255)
                                                ccgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAAc                                                                              <  2:207677/150‑1 (MQ=255)
                                                     ccGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTg                                                                        >  2:190257/1‑151 (MQ=255)
                                                        ttCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcccgcctgccagcagcagttccagcatgttatccatcgaactggag                                                                     >  2:315787‑M2/1‑105 (MQ=255)
                                                             cAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCttt                                                                <  2:113649/151‑1 (MQ=255)
                                                               aCGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcccgcctggcagcagcagttccagcatgttatccatcgaactggagtcagagc                                                              >  2:186020‑M2/1‑98 (MQ=255)
                                                                cGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGAc                                                             <  1:8914/151‑1 (MQ=255)
                                                                  aaaCCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACt                                                            <  1:343318/150‑1 (MQ=255)
                                                                     ccGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTaa                                                        >  1:448599/1‑151 (MQ=255)
                                                                     ccGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTaa                                                        >  1:76564/1‑151 (MQ=255)
                                                                       ggCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTaa                                                        >  1:416806/1‑149 (MQ=255)
                                                                        gCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGCCAGGCGGGtggtactaatagacgcatggcgcggatgatatccatcccgcctgccagcagcagttccagcatgttatccatcgaactggagtcagagccggtttcgt                                                     >  2:186020‑M1/1‑53 (MQ=255)
                                                                         ctctGACTCCAGTTCGATGGATAACATGCTGTAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTc                                                     <  1:525147/150‑1 (MQ=255)
                                                                               cTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGtggtactaatagacgcatggcgcggatgatatccatcccgcctgccagcagcagttccagcatgttatccatcgaactggagtcagagccggtttcgttgacgaa                                              >  2:315787‑M1/1‑46 (MQ=255)
                                                                               cTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCat                                              >  1:496906/1‑151 (MQ=255)
                                                                                    gTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGa                                         >  1:171845/1‑151 (MQ=255)
                                                                                           ggATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTgg                                  >  1:80296/1‑151 (MQ=255)
                                                                                           ggATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTgg                                  >  1:54848/1‑151 (MQ=255)
                                                                                           ggATAACATGCTGGAACTGCTGCTGCCAGGCGGGtggtactaatagacgcatggcgcggatgatatccatcccgcctgccagcagcagttccagcatgttatccatcgaactggagtcagagccggtttcgttgacgaacggtgcggcgt                                   >  2:216953‑M1/1‑34 (MQ=255)
                                                                                               aaCATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcccgcctgccagcagcagttccagcatgttatccatcgaactggagtcagagccggtttcgttgacgaacggtgcggcgtcgtgc                              >  2:818‑M2/1‑66 (MQ=255)
                                                                                                aCATGCTGGAACTGCTGCTGGCAGGCGGGtggtactaatagacgcatggcgcggatgatatccatcccgcctgccagcagcagttccagcatgttatccatcgaactggagtcagagccggtttcgttgacgaacggtgcggcgtcgtgca                             >  2:321414‑M1/1‑29 (MQ=255)
                                                                                                             gctgctgGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCAt                 >  1:398414/1‑150 (MQ=255)
                                                                                                                   ggCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATg          >  2:359273/1‑151 (MQ=255)
                                                                                                                         cGGGtggtactaatagacgcatggcgcggatgatatccatcccgcctgccagcagcagttccagcatgttatccatcgaactggagtcagagccggtttcgttgacgaacggtgcggcgtcgtgcaggtcagggataagcggtgtctggaa    <  1:548610‑M1/151‑148 (MQ=255)
                                                                                                                         cGGGtggtactaatagacgcatggcgcggatgatatccatcccgcctgccagcagcagttccagcatgttatccatcgaactggagtcagagccggtttcgttgacgaacggtgcggcgtcgtgcaggtcagggataagcggtgtctggaa    <  1:113647‑M1/151‑148 (MQ=255)
                                                                                                                         cGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGAATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGa     <  1:496974/150‑1 (MQ=255)
                                                                                                                           ggtggtactaatagacgcatggcgcggatgatatccatcccgcctgccagcagcagttccagcatgttatccatcgaactggagtcagagccggtttcgttgacgaacggtgcggcgtcgtgcaggtcagggataagcggtgtctggaatt  >  2:188719‑M1/1‑2 (MQ=255)
                                                                                                                            |                                                                                                                                                     
CGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGG  >  NC_000913/3355482‑3355755

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