New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 2775873 =37 (0.830)5 (0.120) 5/270 5.4 12.4% intergenic (+57/‑46) ypjI/yfjX pseudogene, CP4‑57 putative prophage remnant;Phage or Prophage Related/CP4‑57 prophage; putative antirestriction protein
?NC_000913_3_bme_pgi 2775896 = 36 (0.860)intergenic (+80/‑23) ypjI/yfjX pseudogene, CP4‑57 putative prophage remnant;Phage or Prophage Related/CP4‑57 prophage; putative antirestriction protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCC  >  NC_000913_3_bme_pgi/2775725‑2776019
                                                                                                                                                    |                                                                                                                                                  
ggacgctaTGGAACTCTGAACCGGTCGCCGGGGAAGGTGTTAACCGCTATAGCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATa                                                                                                                                                    <  2:151339/142‑1 (MQ=255)
          gAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTCCCTCATCCGCCCGCCAt                                                                                                                                          >  2:90572/1‑149 (MQ=255)
          gAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCAt                                                                                                                                          >  1:257502/1‑149 (MQ=255)
            aCTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTACTACCTCATACGCCAGCCATcg                                                                                                                                        >  2:435870/1‑149 (MQ=255)
              tttGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACCCCAGCCCTcgcc                                                                                                                                      >  2:59752/1‑149 (MQ=255)
                   aCCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTgg                                                                                                                                 >  2:392338/1‑149 (MQ=255)
                          gCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGttttt                                                                                                                          >  1:602615/1‑149 (MQ=255)
                            cGGTGAATGTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTa                                                                                                                        <  2:46627/149‑1 (MQ=255)
                            cGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTa                                                                                                                        <  1:90572/149‑1 (MQ=255)
                                    gTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGga                                                                                                                <  2:602615/149‑1 (MQ=255)
                                    gTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGga                                                                                                                <  2:275101/149‑1 (MQ=255)
                                           ccGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATAc                                                                                                         <  1:392338/149‑1 (MQ=255)
                                             gCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAACCGCGCTTTATTACCCCATACGCCAGCCATCGCCGCTGGCGATTTATATTCACGGAGACATAccc                                                                                                       >  1:707707/1‑149 (MQ=255)
                                               tataTCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACCTTTCACCTTCCCAAACGAGCTTTATTACCTCATAAGACAGCTAACGCCGCTGGCGCTTTTTCTGCACGGAGACATACCCac                                                                                                     >  2:276813/1‑148 (MQ=255)
                                                        aGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAacaca                                                                                            <  1:452440/149‑1 (MQ=255)
                                                         ggTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCCCGGAGAAATACCCATCACAccacag                                                                                           >  2:121434/1‑149 (MQ=255)
                                                          gTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTGCGAATCATTTACCAGCAACATTTCACCTTCCGAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAacacaga                                                                                          >  1:361750/1‑149 (MQ=255)
                                                          gTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAacacaga                                                                                          <  2:587446/149‑1 (MQ=255)
                                                            cagggcgccGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCCCCTTCCCAAACGCGCTTTATTACCTCATACGTCCGGCATCGGCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAacacagaca                                                                                        <  2:198674/142‑1 (MQ=38)
                                                             aGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGCACTACCCATGACAacacagacac                                                                                       >  1:694571/1‑149 (MQ=255)
                                                             aGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTAACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTGACGGAGACATGCCCATGCCAgaacagacac                                                                                       >  1:97366/1‑149 (MQ=38)
                                                                 aTCCGCACCGCCACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTa                                                                                   <  2:182920/149‑1 (MQ=255)
                                                                  tCCGCACCGATAACTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTAc                                                                                  <  1:423461/149‑1 (MQ=255)
                                                                      cACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGAcc                                                                              <  2:531147/149‑1 (MQ=255)
                                                                        ccGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATTACAACACAGACACAGTACGACCTc                                                                            <  1:669599/149‑1 (MQ=255)
                                                                         cGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCGCGGAGACATACCCATGACAacacagacaca                                                                                      <  2:230335/138‑1 (MQ=255)
                                                                         cGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCAACTTCCCAAACGCGCTTTATCACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCg                                                                           >  1:265262/1‑149 (MQ=255)
                                                                         cGACACCCTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCGCGGAGACATACCCATGACAacacagacaca                                                                                      >  1:230335/1‑138 (MQ=255)
                                                                            cacCCTCCCCCGCTGACTCCGAATCAGTTACCAGCTACATTTCACCTTCCCACACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTAc                                                                        <  2:752294/149‑1 (MQ=255)
                                                                                  cccccGCTGACTCCGAATCATTTACCACCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACCCTGACACACTACGCGCTCGTACCCCCTa                                                                  >  2:289697/1‑149 (MQ=37)
                                                                                                     aTTTACCAGCAACATTTCACCTTCCCAAACGCGCGCTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGCATttgagtt                                               >  1:746783/1‑149 (MQ=255)
                                                                                                                           tCCCAAACGCGCTTTAATACCTCATACGCCAGCCATCGCCGGTTCCGGTGTTTATTAACGGAACCACGCCAAATACAACACAGACACAGTACGACCTCGTACCCGCTCACGGATCCGAATTTGAGTTCAGCGTACCGCCCGTACCCCAt                         >  2:164873/1‑149 (MQ=255)
                                                                                                                                               cTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTg     >  2:77351/1‑149 (MQ=255)
                                                                                                                                               cTCATACGCCAGCCATCGCCGCTGGCGTTTTTTAAACAGGGGGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTg     <  1:144117/149‑1 (MQ=255)
                                                                                                                                                  aTACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATCTCCGCCTCTGGcc  >  1:64020/1‑149 (MQ=255)
                                                                                                                                                   tACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACg                                                               >  2:495375/1‑87 (MQ=255)
                                                                                                                                                   tACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACg                                                               <  1:495375/87‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCATCCGCACCGACACCTTCCCCCGCTGACTCCGAATCATTTACCAGCAACATTTCACCTTCCCAAACGCGCTTTATTACCTCATACGCCAGCCATCGCCGCTGGCGTTTTTTATTCACGGAGACATACCCATGACAACACAGACACAGTACGACCTCGTACCCGCTAACGGATCCGAATTTGAGTTGAGCGTAACGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCC  >  NC_000913_3_bme_pgi/2775725‑2776019

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