New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1964259 =34 (0.960)4 (0.130) 3/114 2.7 11.3% coding (792/2079 nt) flhA putative flagellar export pore protein
?NC_000913 1964341 = 33 (1.060)coding (710/2079 nt) flhA putative flagellar export pore protein
Rejected: Frequency below/above cutoff threshold.

CCAGCGCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATG  >  NC_000913/1964276‑1964406
                                                                 |                                                                 
ccAGCGCCTGAATTTGTTCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCgtgt                                                                 <  1:107330/68‑1 (MQ=255)
ccAGCGCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCgtgt                                                                 <  2:111120/68‑1 (MQ=255)
ccAGCGCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCgtgt                                                                 <  2:1096124/68‑1 (MQ=255)
ccAGCGCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGAgtgt                                                                 >  1:1247873/1‑68 (MQ=255)
  ctattgaccattggcgactgtcttgtggcacaaattccggcgctcggtatttctaccgcctcgTGTcc                                                               <  1:1097989‑M2/5‑1 (MQ=255)
      ccGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATg                                                           <  1:1011558/68‑1 (MQ=255)
            tttGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATg                                                     <  1:775298/68‑1 (MQ=255)
              ggcgacggtctggtggcacaaattccgtcgctgggtatttctaccgccgcgTGTCCCATGCTCATGcc                                                   <  1:907694‑M2/17‑1 (MQ=255)
                tgCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCAt                                                 <  1:923041/68‑1 (MQ=255)
                cgacggtcttgtgtcacaaattccggcgctggtgatttctaccgcctcgTGTCCCATGCTCATGCCAt                                                 <  1:1309874‑M2/19‑1 (MQ=255)
                cgacggtctggcggcacaatttccggcgctggtgatttctaccgccgcgTGTCCCATGCTCATGCCAt                                                 <  2:779753‑M2/19‑1 (MQ=255)
                  ccaccaGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGt                                               <  1:398470/68‑1 (MQ=255)
                  ccaccaGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGt                                               <  1:980980/68‑1 (MQ=255)
                         aCCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTgcagcc                                        >  1:200019/1‑67 (MQ=255)
                         aCCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTgcagca                                        >  1:1317720/1‑68 (MQ=255)
                         aCCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTgcagca                                        >  1:897403/1‑68 (MQ=255)
                         aCCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCACATGCTCAAGCAATATTgcagca                                        >  2:1213333/1‑68 (MQ=255)
                           cGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACg                                      >  1:71750/1‑68 (MQ=255)
                             tCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGcc                                    <  2:707716/68‑1 (MQ=255)
                                     ggTCAATAGCGTATAACTTTCCGCAGCGTGTCCCATGATCATGCCAAGTTGCAGCACGACGACCAggg                            >  1:964351/1‑66 (MQ=255)
                                        cAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAAcc                         <  2:1059194/68‑1 (MQ=255)
                                        cAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAAcc                         <  1:1093134/68‑1 (MQ=255)
                                          aTAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGCCCAGCAACccg                       >  1:352388/1‑68 (MQ=255)
                                          aTAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACccg                       >  1:160888/1‑68 (MQ=255)
                                            aGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACccgcc                     >  2:259507/1‑68 (MQ=255)
                                            aGCGTATAACTTTCCGCCGCGTGACCCATGCTCATCCCATGTTGAAGCACGCCGACCAGCAACccgac                     >  2:641874/1‑66 (MQ=255)
                                            aGCGTATAACAATCCGCCGCGTGTCCCAAGCTCATGCCGTGTTGCAGCACGCCGACCAGCACCccgac                     >  2:392052/1‑66 (MQ=37)
                                             gCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACccgccg                    <  2:131248/68‑1 (MQ=255)
                                                  taACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAAt               <  1:98536/68‑1 (MQ=255)
                                                  taACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAAt               <  1:686748/68‑1 (MQ=255)
                                                    aCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGt             <  2:159070/68‑1 (MQ=255)
                                                    aCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGt             <  2:71750/68‑1 (MQ=255)
                                                    aCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGt             >  2:228613/1‑68 (MQ=255)
                                                    aCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGt             >  1:832801/1‑68 (MQ=255)
                                                    aCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGt             <  2:26160/68‑1 (MQ=255)
                                                    aCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGt             >  2:590577/1‑68 (MQ=255)
                                                    aCTTTCCGCCGAGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGAACAGCACCCCGCCGCCAATGt             >  2:531544/1‑68 (MQ=255)
                                                         ccgccgcgTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATg        >  2:949048‑M2/9‑68 (MQ=255)
                                                           gccgcgTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGAc      <  2:1158686‑M2/62‑1 (MQ=255)
                                                           gccgcgTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGAc      <  1:228614‑M2/62‑1 (MQ=255)
                                                            ccgcgTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGAcc     <  1:903145‑M2/63‑1 (MQ=255)
                                                            ccgcgTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGAcc     <  1:1068385‑M2/63‑1 (MQ=255)
                                                              gcgTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCAt   <  2:160888‑M2/65‑1 (MQ=255)
                                                               cgTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATg  <  2:832799‑M2/66‑1 (MQ=255)
                                                                 |                                                                 
CCAGCGCCGGAATTTGTGCCACCAGACCGTCGCCAATGGTCAATAGCGTATAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATG  >  NC_000913/1964276‑1964406

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