New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5224548 =5 (0.060)26 (0.310) 20/252 4.3 66.1% intergenic (‑204/+175) PP_4603/PP_4604 transposase/DMT superfamily permease
?NC_002947 5224584 = 22 (0.270)intergenic (‑240/+139) PP_4603/PP_4604 transposase/DMT superfamily permease
Rejected: Coverage evenness skew score above cutoff.

CCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGCCCAGGCAACCATCATCTATTT  >  NC_002947/5224537‑5224728
                                               |                                                                                                                                                
ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCgg                                            >  1:802505/1‑150 (MQ=31)
ccTAACAAAGCAAAAGACGTGTGGCTCTACCTGGCCAACCGCCGCGGGTGTTTCAAAATGGCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCgg                                            >  2:306183/1‑150 (MQ=11)
 cTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGt                                           <  2:143342/150‑1 (MQ=32)
        aGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTAATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCcg                                    >  2:1861685/1‑150 (MQ=38)
         gCAACAGGCGGCTGGCCCTACCGGGCCAAGCGCCTTTTCTGTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCcgc                                   <  2:470953/150‑1 (MQ=255)
             aaGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGaaa                               <  1:441595/150‑1 (MQ=255)
              aGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAg                              >  1:1522579/1‑150 (MQ=255)
              aGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAg                              <  1:980036/150‑1 (MQ=255)
              aGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTACTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAg                              <  2:87035/150‑1 (MQ=255)
               ggcggcTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAgg                             <  1:2130489/150‑1 (MQ=255)
                 cggcTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGggc                           <  2:759086/150‑1 (MQ=255)
                    cTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGggccgg                        >  2:435585/1‑150 (MQ=255)
                       gCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGccc                     <  1:435585/150‑1 (MQ=255)
                         ccTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCcgc                                   <  2:1187760/134‑1 (MQ=39)
                         ccTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCcgc                                   >  1:1187760/1‑134 (MQ=39)
                          cTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGCCCAgg                  >  2:52898/1‑150 (MQ=255)
                               gggCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGCCCAGGCAAcc             >  1:742804/1‑150 (MQ=255)
                                  ccAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGCCCAGGCAACcatc          <  2:1836217/150‑1 (MQ=255)
                                  ccAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGCCCAGGCAACcatc          <  2:802505/150‑1 (MQ=255)
                                    aaCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGCCCAGGCAACcatcat        <  2:1678926/150‑1 (MQ=255)
                                     aCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGCCCa                    <  2:1979746/137‑1 (MQ=255)
                                     aCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGCCCa                    >  1:1979746/1‑137 (MQ=255)
                                        gCCTTTTCTTTTTCAAAATGTCAGTGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGccc                     <  2:1482914/133‑1 (MQ=255)
                                        gCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGccc                     >  1:1482914/1‑133 (MQ=255)
                                          cttttcttttTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGCCCAGGCAACCATCATCTAttt  <  2:335732/150‑1 (MQ=255)
                                               |                                                                                                                                                
CCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGCGGTTTACCCGCGAAAGGGCCGGCCCAGGCAACCATCATCTATTT  >  NC_002947/5224537‑5224728

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