New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 346713833 (1.100)4 (0.140) 4/206 2.9 12.2% intergenic (‑147/+22) bfd/chiA bacterioferritin‑associated ferredoxin/periplasmic endochitinase
?NC_000913 = 3467144 27 (0.970)intergenic (‑153/+16) bfd/chiA bacterioferritin‑associated ferredoxin/periplasmic endochitinase
Rejected: Frequency below/above cutoff threshold.

GACTCCCGTATACTTTCTTCAAAAATAAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTCGCGCGCGATTGAC  >  NC_000913/3466995‑3467280
                                                                                                                                               |                                                                                                                                              
gaCTCCCGTATACTTTCTTCAAAAATAAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGaa                                                                                                                                           <  4:133828/149‑1 (MQ=255)
         ataCTTTCTTCAAAAATAAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGgtatcgc                                                                                                                                   <  2:115868‑M2/148‑8 (MQ=255)
                    aaaaaTAAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTa                                                                                                                        <  1:115868/148‑1 (MQ=255)
                     aaaaTAAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTAtt                                                                                                                      >  8:58722/1‑149 (MQ=255)
                                         ttATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTa                                                                                                                 >  6:91801/1‑134 (MQ=255)
                                         ttATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTa                                                                                                                 <  5:91801/134‑1 (MQ=255)
                                                     cGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATAAGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGaa                                                                                                   <  1:165608/136‑1 (MQ=255)
                                                     cGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATAAGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGaa                                                                                                   >  2:165608/1‑136 (MQ=255)
                                                      ggAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAg                                                                                                                                          <  8:108134/96‑1 (MQ=255)
                                                        aaaGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGaa                                                                                                                                           >  5:58987/1‑93 (MQ=255)
                                                        aaaGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGaa                                                                                                                                           <  6:58987/93‑1 (MQ=255)
                                                        aaaGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCg                                                                                   >  4:42655/1‑149 (MQ=255)
                                                              tctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGaa                                                                                                                                           >  3:132361/1‑87 (MQ=255)
                                                              tctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGaa                                                                                                                                           <  4:132361/87‑1 (MQ=255)
                                                                ttttttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGc                                                                                    <  3:49371/140‑1 (MQ=255)
                                                                ttttttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGc                                                                                    >  4:49371/1‑140 (MQ=255)
                                                                 tttttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTa                                                                                 >  5:111274/1‑142 (MQ=255)
                                                                 tttttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTa                                                                                 <  6:111274/142‑1 (MQ=255)
                                                                        gCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAaca                                                                                                                              >  6:43166/1‑90 (MQ=255)
                                                                        gCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAaca                                                                                                                              <  5:43166/90‑1 (MQ=255)
                                                                              aaTAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACcttccagtatcgccgcatttttt                                                                                                                         >  6:149342‑M1/1‑66 (MQ=255)
                                                                              aaTAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACcttccagtatcgccgcatttttt                                                                                                                         <  5:149342‑M1/89‑24 (MQ=255)
                                                                                        tGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGt                                                    <  3:42655/148‑1 (MQ=255)
                                                                                          aaTAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTaaa                                                                                                            <  2:70421/90‑1 (MQ=255)
                                                                                               aaaTAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACcttccagtatcgcc                                                                                                                                  >  6:66679‑M1/1‑49 (MQ=255)
                                                                                               aaaTAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACcttccagtatcgcc                                                                                                                                  <  5:66679‑M1/63‑15 (MQ=255)
                                                                                                          tACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAg                                    >  3:84202/1‑146 (MQ=255)
                                                                                                          tACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAg                                    <  4:84202/146‑1 (MQ=255)
                                                                                                                         aTTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGAttt                                                                                              >  1:123830/1‑73 (MQ=255)
                                                                                                                         aTTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGAttt                                                                                              <  2:123830/73‑1 (MQ=255)
                                                                                                                         aTTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCg                                                                            >  7:59636/1‑91 (MQ=255)
                                                                                                                         aTTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCg                                                                            <  8:59636/91‑1 (MQ=255)
                                                                                                                             aCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCg                                                                            <  4:147907/87‑1 (MQ=255)
                                                                                                                             aCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCg                                                                            >  3:147907/1‑87 (MQ=255)
                                                                                                                             aCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTc               <  5:165471/148‑1 (MQ=255)
                                                                                                                               aaaaaaTGCGGCGATACTGGAAGgtatcgccgcattttttgtcaatgggctacttttggtagcccattatttttattcatatttctatt                                                                        <  5:149342‑M2/89‑67 (MQ=255)
                                                                                                                               aaaaaaTGCGGCGATACTGGAAGgtatcgccgcattttttgtcaatgggctacttttggtagcccattatttttattcatatttctatt                                                                        >  6:149342‑M2/1‑23 (MQ=255)
                                                                                                                                        ggcgATACTGGAAGgtatcgccgcattttttgtcaatgggctacttttggtagcccattattt                                                                                         <  5:66679‑M2/63‑50 (MQ=255)
                                                                                                                                        ggcgATACTGGAAGgtatcgccgcattttttgtcaatgggctacttttggtagcccattattt                                                                                         >  6:66679‑M2/1‑14 (MQ=255)
                                                                                                                                        ggcgATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTCGCGCGCGa       >  2:169057/1‑145 (MQ=255)
                                                                                                                                        ggcgATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTCGCGCGCGa       <  1:169057/145‑1 (MQ=255)
                                                                                                                                         gcgATACcttccagtatcgccgcattttttgtcaatgggctacttttggtagcccattatttttattcatatttctattctttgcttaaaaaagagtctttccgtaattgaaataagaactattttcatttatttttgaagaaagtat   <  2:115868‑M1/148‑142 (MQ=255)
                                                                                                                                           gATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTCGCGCGCGATTGAc  >  1:139381/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                              
GACTCCCGTATACTTTCTTCAAAAATAAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTTAGTCAGGCCGTGGAATTCCGGAGACAGGTTAGTGCCGCCCGGCAGGTCGCGCGCGATTGAC  >  NC_000913/3466995‑3467280

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