New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA 3804946 =22 (0.760)4 (0.140) 4/278 3.7 16.2% coding (171/936 nt) waaS lipopolysaccharide rhamnose:KdoIII transferase; lipopolysaccharide core biosynthesis protein
?NC_000913_3_pae_tpiA 3804980 = 20 (0.710)coding (137/936 nt) waaS lipopolysaccharide rhamnose:KdoIII transferase; lipopolysaccharide core biosynthesis protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAAAACTTTATGACAGTTAATATTTTTTAATATTTGTTTATTGATACCTTTATTAACAACACTGTCATTAATGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTT  >  NC_000913_3_pae_tpiA/3804832‑3805125
                                                                                                                                                    |                                                                                                                                                  
cAAAACTTTATGACAGTTAATATTTTTTAATATTTGTTTATTGATACCTTTATTAACAACACTGTCATTAATGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTatccatcc                                                                                                                                                    <  2:151735/149‑1 (MQ=255)
    aCTTTATGACAGTTAATATTTTTTAATATTTGTTTATTGATACCTTTATTAACAACACTGTCATTAATGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTAt                                                                                                                                                >  2:97182/1‑149 (MQ=255)
                  aaTATTTTTTAATATTTGTTTATTGATACCTTTATTAACAACACTGTCATTAATGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATAtttttt                                                                                                                                  >  1:165679/1‑149 (MQ=255)
                   atatTTTTTAATATTTGTTTATTGATACCTTTATTAACAACACTGTCATTAATGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATAttttttt                                                                                                                                 <  2:402033/149‑1 (MQ=255)
                                      tATTGATACCTTTATTAACAACACTGTCATTAATGATAATTATGTCATCTTTTTTTATATTATTGAGGTAATGCTTGATGAATAATCTCCCCAGCCATGCATTATCCATCCCTCTTTTTGATATTTTTTTATTAATTCACTTATAtttt                                                                                                              >  2:282456/1‑149 (MQ=255)
                                               cTTTATTAACAACACTGTCATTAATGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGggtt                                                                                                     >  1:483940/1‑149 (MQ=255)
                                                tttATTAACAACACTGTCATTAATGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGACATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGttttg                                                                                                    >  2:405231/1‑149 (MQ=255)
                                                                   ttAATGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGTTTTTTTTATATTTTTTATCTTGTa                                                                                 >  2:105120/1‑149 (MQ=255)
                                                                        gATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCAATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTATTTGGTGTTATTATATTTTTTATCTTTTACtttt                                                                           >  2:405477/1‑149 (MQ=255)
                                                                                    attttttttAATATCATTGAAGTAAGGCTTGATGAATACTCTCCCCAGCCATTCATTATCTATCCCTATTTTTGTTATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACt                                                               <  2:359287/146‑1 (MQ=255)
                                                                                             aaaatCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAGTCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTTCAGGTACTGAATCAtt                                                       <  2:414698/146‑1 (MQ=255)
                                                                                                atCATTGAGGTAATGCTTGAGGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTGTATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTa                                                   >  1:75258/1‑149 (MQ=255)
                                                                                                   aTTGAGGTAGTGGTTGATGAATAATCGCCCCAGTCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAAtt                                                <  2:94506/149‑1 (MQ=255)
                                                                                                   aTTGAGGTAATTCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAAtt                                                <  2:296144/149‑1 (MQ=255)
                                                                                                               cTTGATGAATAATTGCCCCAGCCTTGAATTATCCATCCCTGTTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCaa                                    <  2:111366/149‑1 (MQ=255)
                                                                                                                ttGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTCTCAGCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAAc                                   >  1:494983/1‑149 (MQ=255)
                                                                                                                   atgaATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAGTCCAGTGATATTTTTTTGGTGTTGTTATATTTTTTTTCTTGTACTTTTTTTTCAGTTACTGAATCATTTTTAATTCATAATCTCCAACCCa                                >  1:476354/1‑149 (MQ=255)
                                                                                                                        taatCGCCCCAGCCATGTATTATCTATTCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTta                           <  2:117314/149‑1 (MQ=255)
                                                                                                                             gCCCCAGTCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTataaaa                      <  2:187743/149‑1 (MQ=255)
                                                                                                                                             attcatccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCAtaataa      <  1:288980/146‑1 (MQ=255)
                                                                                                                                              tccatgcCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATaaa     <  2:486457/149‑1 (MQ=255)
                                                                                                                                                 atccCTATTTTTGTTTTTTTTTTATTAAATCAGTGATATTTTTTGTGTTTTGTTATATTTTTTATCTTGTACTTTTTTTTTAGTTACTGAATCATTTTTAATTCATAATTTGCAACCCAATTTATAAAATAAATAGTCATAATAAattt  >  2:311842/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAAAACTTTATGACAGTTAATATTTTTTAATATTTGTTTATTGATACCTTTATTAACAACACTGTCATTAATGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTT  >  NC_000913_3_pae_tpiA/3804832‑3805125

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.