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.

AATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTC  >  NC_000913/3355494‑3355790
                                                                                                                                                    |                                                                                                                                                    
aaTTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcc                                                                                                                                                    >  2:188719‑M2/1‑149 (MQ=255)
  ttCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcccg                                                                                                                                                  <  1:548610‑M2/151‑5 (MQ=255)
  ttCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcccg                                                                                                                                                  <  1:113647‑M2/151‑5 (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)
                             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:496908/1‑149 (MQ=255)
                                    ccgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGaa                                                                                                                  >  2:241462/1‑149 (MQ=255)
                                    ccgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAAc                                                                                                                 <  2:207677/150‑1 (MQ=255)
                                    ccgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAAc                                                                                                                 <  2:139003/150‑1 (MQ=255)
                                    ccgcACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAAc                                                                                                                 <  1:138644/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:76564/1‑151 (MQ=255)
                                                         ccGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTaa                                                                                           >  1:448599/1‑151 (MQ=255)
                                                           ggCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTaa                                                                                           >  1:416806/1‑149 (MQ=255)
                                                             ctctGACTCCAGTTCGATGGATAACATGCTGTAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTc                                                                                        <  1:525147/150‑1 (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)
                                                                                   aaCATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCAcccgcctgccagcagcagttccagcatgttatccatcgaactggagtcagagccggtttcgttgacgaacggtgcggcgtcgtgc                                                                 >  2:818‑M2/1‑66 (MQ=255)
                                                                                                 gctgctgGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCAt                                                    >  1:398414/1‑150 (MQ=255)
                                                                                                       ggCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATg                                             >  2:359273/1‑151 (MQ=255)
                                                                                                             cGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGAATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGa                                        <  1:496974/150‑1 (MQ=255)
                                                                                                                 atggatATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGt                                    <  2:548620/150‑1 (MQ=255)
                                                                                                                                 ccATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTaa                   >  2:399352/1‑151 (MQ=255)
                                                                                                                                    tGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCt                <  1:41932/151‑1 (MQ=255)
                                                                                                                                     gCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCt                >  2:83359/1‑150 (MQ=255)
                                                                                                                                         cTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACgg                                     <  1:266214/125‑1 (MQ=255)
                                                                                                                                         cTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACgg                                     >  2:266219/1‑125 (MQ=255)
                                                                                                                                          tATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCg          >  2:90599/1‑151 (MQ=255)
                                                                                                                                                  ccaccCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTc  <  2:507380/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCATGCATATGGAGCCGTGGGATGGCCCGGCGGGCATCGTGATGTCCGACGGTCGTTTTGCCGCCTGTAACCTCGACCGTAACGGTC  >  NC_000913/3355494‑3355790

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.