New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3209524 =19 (0.650)3 (0.120) 3/150 2.0 12.8% coding (327/765 nt) nfeF NADPH‑dependent ferric chelate reductase
?NZ_CP009273 3209606 = 24 (0.930)coding (245/765 nt) nfeF NADPH‑dependent ferric chelate reductase
Rejected: Frequency below/above cutoff threshold.

CCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGA  >  NZ_CP009273/3209523‑3209695
                                                                                   |                                                                                         
ccccACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGt                                                                                      <  1:566783/89‑1 (MQ=255)
ccccACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGt                                                                                      >  2:566783/1‑89 (MQ=255)
     CCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTcc                                                                                <  1:441954‑M1/90‑1 (MQ=255)
               tAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGa                                                                      <  1:619105/90‑1 (MQ=255)
                      atcaatcGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCctc                                                               >  1:737421/1‑90 (MQ=255)
                       ccgctgtatgacgaactacgccatgacctggcgattgtttttttttttcacgacggtgggGGTGGTCCTTCCGGCCAGACGATGCCctct                                                              <  1:554671‑M2/30‑1 (MQ=255)
                            cGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGt                                                         <  2:788917/90‑1 (MQ=255)
                                 gTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCg                                                    >  2:328869/1‑90 (MQ=255)
                                      tGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTggc                                               <  1:156547/90‑1 (MQ=255)
                                      tGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTggc                                               >  2:285613/1‑90 (MQ=255)
                                      ctacgccatgaactggcgattgatttctttattcccgacggtgggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTggc                                               <  2:433599‑M2/45‑1 (MQ=255)
                                       tgcgccatgaactggcgattgatttctttattcacggcggtgggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTggcg                                              <  2:219899‑M2/46‑1 (MQ=255)
                                        gCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTggcgg                                             >  1:137049/1‑90 (MQ=255)
                                           tAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGcac                                          >  2:524656/1‑90 (MQ=255)
                                                    aTACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGc                                 >  1:622950/1‑90 (MQ=255)
                                                       cAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTc                              >  2:771971/1‑90 (MQ=255)
                                                         gcggcgTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAg                            >  1:307154/1‑90 (MQ=255)
                                                           ggcgTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGt                          >  1:70808/1‑90 (MQ=255)
                                                           ggcgTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGt                          <  1:621982/90‑1 (MQ=255)
                                                            gcgTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGtt                         >  1:595092/1‑90 (MQ=255)
                                                              gTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGa                       >  1:110518/1‑90 (MQ=255)
                                                              gTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGa                       >  1:462124/1‑90 (MQ=255)
                                                                 taGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGa                    >  1:571108/1‑90 (MQ=255)
                                                                   gTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGaaa                  >  1:501812/1‑90 (MQ=255)
                                                                      aCGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAaagaa               >  2:27248/1‑90 (MQ=255)
                                                                        gcgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCGGGTTGAGGAAagaaga             >  2:331886/1‑90 (MQ=255)
                                                                        gcgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGGGGAAagaaga             >  2:377694/1‑90 (MQ=255)
                                                                        gcgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAagaaga             >  1:114620/1‑90 (MQ=255)
                                                                         cgacggtgggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAagag            >  2:264315‑M2/11‑90 (MQ=255)
                                                                         cgacggtgggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAagag            <  2:307154‑M2/80‑1 (MQ=255)
                                                                                   cGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGa  <  2:431157/90‑1 (MQ=255)
                                                                                   cGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGa  <  1:141933/90‑1 (MQ=255)
                                                                                   |                                                                                         
CCCCACCGTCGTGAATAAAGAAATCAATCGCCAGTTCATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGA  >  NZ_CP009273/3209523‑3209695

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