New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA 2729630 =NA (NA)11 (0.290) 7/274 4.0 NA noncoding (1528/1542 nt) rrsG 16S ribosomal RNA of rrnG operon
?NC_000913_3_pae_tpiA 2729652 = NA (NA)noncoding (1506/1542 nt) rrsG 16S ribosomal RNA of rrnG operon
Rejected: Coverage evenness skew score above cutoff.

GTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAA  >  NC_000913_3_pae_tpiA/2729483‑2729714
                                                                                                                                                    |                                                                                    
gtAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAGGGAGGTGATCCaa                                                                                     >  1:78717/1‑149 (MQ=255)
  aGGGTGAGCTTTCATTAATAGATAGCGAACTGTCTTATGCTCTTCAGCGTCCCTCCCAACGTGTAAACGCCTTG‑CTATTCACTTTTCATCATACAAGCTGTGTGAGCACTGCCAAGTACGTTTCTTTAATGTAAGGAGGTGATCCAAcc                                                                                   <  2:284338/149‑1 (MQ=14)
  aGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGCGGTGATCCAAcc                                                                                   >  1:151391/1‑149 (MQ=255)
     gTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCa                                                                                <  1:579071/149‑1 (MQ=40)
        aGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGCTGATCCAACCGCAGGt                                                                             >  2:182624/1‑149 (MQ=255)
           tttCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCAAACGCGTAAAAGCCTTG‑CTATTAACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGCTGACCAAACCGCAGGTTcc                                                                          >  2:99492/1‑149 (MQ=14)
            ttCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTccc                                                                         <  2:183311/149‑1 (MQ=255)
                ttAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTAc                                                                     >  1:656476/1‑149 (MQ=255)
                             aaCGGCCTTATTCTCTTCAGCCTCACTCCCAAAGCGTCACAGCCATGTCCACT‑ATGATTCGTCAGTGAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTa                                                        <  2:665482/149‑1 (MQ=11)
                              aCGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTAc                                                       <  1:352429/149‑1 (MQ=40)
                                    ttATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGGTCCCCTACGGTTACCTTGTTACGACTTc                                                 >  2:371796/1‑149 (MQ=255)
                                            ttCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGt                                         >  1:295984/1‑149 (MQ=255)
                                               aGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCAt                                      <  1:544938/149‑1 (MQ=40)
                                               aGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTAAt                                      >  1:122487/1‑149 (MQ=255)
                                                   tCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAAt                                  >  2:519700/1‑149 (MQ=255)
                                                    cACTCCCAACGCGTAAACGACTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCCGGTTCCCCTACGGTTACCTTGTTACGTCTTCACCCCAGGCAAGacgt                                 >  2:83467/1‑146 (MQ=14)
                                                         ccAACGCGTAAACGCCTTG‑CTATTCGCTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACaaa                            <  1:535962/149‑1 (MQ=255)
                                                         ccAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACaaa                            <  1:519700/149‑1 (MQ=40)
                                                              gcgTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGt                       >  1:377952/1‑149 (MQ=39)
                                                                 tAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAg                    >  2:271161/1‑149 (MQ=255)
                                                                        cTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGCTAGGCGCCCTc             <  1:426990/149‑1 (MQ=255)
                                                                              tATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGaa        <  1:632761/149‑1 (MQ=35)
                                                                                 tCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGCAGGt     >  2:312496/1‑149 (MQ=17)
                                                                                  cACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCGTCTTTAATGTAAGGAGGGGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGtt    <  2:247853/149‑1 (MQ=12)
                                                                                    cTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTaa  <  2:553267/149‑1 (MQ=11)
                                                                                                                                                    |                                                                                    
GTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTG‑CTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAA  >  NC_000913_3_pae_tpiA/2729483‑2729714

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