New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 41921776 (0.370)4 (0.260) 4/242 1.6 35.1% coding (459/771 nt) thiG thiamine biosynthesis ThiGH complex subunit
?NC_000913_3_pae_pgi = 4192203 9 (0.580)coding (433/771 nt) thiG thiamine biosynthesis ThiGH complex subunit

AGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAG  >  NC_000913_3_pae_pgi/4192055‑4192337
                                                                                                                                                    |                                                                                                                                      
aGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAg                                                                                                                                       <  1:83849‑M2/150‑2 (MQ=255)
                          gCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACaa                                                                                                             >  2:6009/1‑150 (MQ=255)
                                                        tgtgGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCAc                                                                               >  2:147316/1‑150 (MQ=255)
                                                                    gATAATAAACGCCAGCATGGCGCGCGTTTCCAGACTCTGATTCGAGCCAAACGGCGCGCCGAGCGGCATAACCGCTGCCCAGCCGACCTCATCCAGACGATTCCCCAATACCGGAACGACCCCACAGCAAGACAGCACGACAAATccccg                                                                   >  2:138399/1‑148 (MQ=255)
                                                                                    aTGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTcg                                                   >  4:108002/1‑150 (MQ=255)
                                                                                                        cTGATTCGACCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGAt                               >  1:137038/1‑150 (MQ=255)
                                                                                                               gAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAgccgattggctcgaatcagggactggaaacccgcgccatgctggagattattatccagcaggccacagtgccggtggttgtcgatgctggc                                             <  1:169889‑M2/129‑92 (MQ=255)
                                                                                                               gAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAgccgattggctcgaatcagggactggaaacccgcgccatgctggagattattatccagcaggccacagtgccggtggttgtcgatgctggc                                             >  2:169889‑M2/1‑38 (MQ=255)
                                                                                                                       cGGCGCGCCGAGCGGCATCACCGCTGCACAgccgattggctcgaatcagggactggaaacccgcgccatgctggagattattatccagcaggccacagtgccggtggttgtc                                                      <  3:60315‑M2/112‑83 (MQ=255)
                                                                                                                       cGGCGCGCCGAGCGGCATCACCGCTGCACAgccgattggctcgaatcagggactggaaacccgcgccatgctggagattattatccagcaggccacagtgccggtggttgtc                                                      >  4:60315‑M2/1‑30 (MQ=255)
                                                                                                                         gcgcgcCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAAcagcc              <  1:147316/150‑1 (MQ=255)
                                                                                                                                     gCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAg  >  4:2397/1‑150 (MQ=255)
                                                                                                                                        tCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCa                                         >  1:128262/1‑108 (MQ=255)
                                                                                                                                        tCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCa                                         <  2:128262/108‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                      
AGCGCCTGCGCGGCATGGCTGGGAACGCCGATGCCAGCATCGACAACCACCGGCACTGTGGCCTGCTGGATAATAATCTCCAGCATGGCGCGGGTTTCCAGTCCCTGATTCGAGCCAATCGGCGCGCCGAGCGGCATCACCGCTGCACAGCCGACTTCTTCCAGACGTTTACACAATACCGGATCGGCCCCGCAGTAAGGCAGCACGACAAATCCCTGTTGTACCAGCGTTTCGGCGGCTTTCAGGGTTTCGATGGGATCGGGCAACAGCCAGCGGGCGTCAG  >  NC_000913_3_pae_pgi/4192055‑4192337

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