New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4449524132 (1.030)6 (0.050) 4/168 NT 4.6% coding (129/531 nt) ppa inorganic pyrophosphatase
?NC_000913 = 4449539 127 (1.050)coding (114/531 nt) ppa inorganic pyrophosphatase
Rejected: Frequency below cutoff threshold.

CGGAGTTGGGACCAGTACGTCAACCGGGTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4449422‑4449524
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcggCGCACTGTTCGTTGACCGCTTCATGTCCACCGCGATGTTCTATCCATGCAACTACGGTTACATCAACCACACCCTGTCTCTGGACGGTGACCCGGTTG  <  NC_000913/4449539‑4449442
                                                                                                                                                                                                         
CGGAGTTGGGACCAGTACGTCAACCGGGTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAG                                                                                                       >  1:1091162/1‑100
CGGAGTTGGGACCAGTACGTCAACCGGGTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAG                                                                                                       >  2:1858463/1‑100
CGGAGTTGGGACCAGTACGTCAACCGGGTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAG                                                                                                       <  2:3332927/100‑1
 GGAGTTGGGACCAGTACGTCAACCGGGTCACCGTCCAGAGACCGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGC                                                                                                      <  2:490683/100‑1
 GGAGTTGGGACCAGTACGTCAACCGGGTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGC                                                                                                      <  1:871386/100‑1
 GGAGTTGGGACCAGTACGTCAACCGGGTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGC                                                                                                      <  2:2812203/100‑1
  GAGTTGGGACCAGTACGTCAACCGGGTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGCG                                                                                                     >  1:2184921/1‑100
   AGTTGGGACCAGTACGTCAACCGGGTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGCGG                                                                                                    >  1:2466022/1‑100
   AGTTGGGACCAGTACGTCAACCGGGTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGCGG                                                                                                    <  1:2580413/100‑1
                          GGTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGCGGCGCACTGTTCGTTGACCGCTTCA                                                                             <  1:2982769/100‑1
                           GTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGCGGCGCACTGTTCGTTGACCGCTTCAT                                                                            <  1:279333/100‑1
                            TCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGCGGCGCACTGTTCGTTGACCGCTTCATG                                                                           <  2:965600/100‑1
                            TCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGCGGCGCACTGTTCGTTGACCGCTTCATG                                                                           <  2:602824/100‑1
                            TCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGCGGCGCACTGTTCGTTGACCGCTTCATG                                                                           <  2:3004891/100‑1
                             CACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGCGGCGCACTGTTCGTTGACCGCTTCATGT                                                                          <  1:1175479/100‑1
                                                                                                   GCGGCGCACTGTTCGTTGACCGCTTCATGTCCACCGCGATGTTCTATCCATGCAACTACGGTTACATCAACCACACCCTGTCTCTGGACGGTGACCCGGT    >  2:1014782/1‑100
                                                                                                    CGGCGCACTGTTCGTTGACCGCTTCATGTCCACCGCGATGTTCTATCCATGCAACTACGGTTACATCAACCACACCCTGTCTCTGGACGGTGACCCGGTT   <  1:1104687/100‑1
                                                                                                     GGCGCACTGTTCGTTGACCGCTTCATGTCCACCGCGATGTTCTATCCATGCAACTACGGTTACATCAACCACACCCTGTCTCTGGACGGTGACCCGGTTG  >  2:147041/1‑100
                                                                                                      GCGCATTGTTCGTTGACCGCTTCATGTCCACCGCGATGTTCTATCCATGCAACTACGGTTACATCAACCACACCCTGTCTCTGGACGGTGACCCGGTTG  <  1:1283250/99‑1
                                                                                                      GCGCACTGTTCGTTGACCGCTTCATGTCCACCGCGATGTTCTATCCATGCAACTACGGTTACATCAACCACACCCTGTCTCTGGACGGTGACCCGGTTG  <  2:2490147/99‑1
                                                                                                                                                                                                         
CGGAGTTGGGACCAGTACGTCAACCGGGTCACCGTCCAGAGACAGGGTGTGGTTGATGTAACCGTAGTTGCATGGATAGAACATCGCGGTGGACATGAAGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4449422‑4449524
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agcggCGCACTGTTCGTTGACCGCTTCATGTCCACCGCGATGTTCTATCCATGCAACTACGGTTACATCAACCACACCCTGTCTCTGGACGGTGACCCGGTTG  <  NC_000913/4449539‑4449442

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