New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi 1411899 =20 (1.040)4 (0.250) 4/238 1.8 24.0% coding (50/936 nt) ttcA tRNA s(2)C32 thioltransferase, iron sulfur cluster protein
?NC_000913_3_vch_pgi = 1434958 9 (0.570)pseudogene (51/51 nt) ttcC pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related

CTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTT  >  NC_000913_3_vch_pgi/1411760‑1412041
                                                                                                                                           |                                                                                                                                              
ctgcAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACACAGACGGTTTTGTAATTTGTTCAGGt                                                                                                                                       >  1:305477/1‑149 (MQ=255)
         gcgcAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGtt                                                                                                                              <  1:292777/149‑1 (MQ=255)
          cgcAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTGTTGttc                                                                                                                             >  1:126957/1‑149 (MQ=255)
                                     agagtggattttctgagcttattgttgtcgctgatgatatttaacacattgaaaaataagtaaaatacttatgagtcagagagttgtgatttttgcccttacTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTg                                                                                                  >  1:151121‑M1/103‑149 (MQ=255)
                                       tttACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTCTAATTTGTTCCGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGcc                                                                                                >  2:289556/1‑148 (MQ=255)
                                            ccccGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCattatt                                                                                           <  1:108883/149‑1 (MQ=255)
                                                ggAGAGGCAAACCATGATGTGATCGCCTTCTTCAATCTTATTGAAGTTAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAg                                                                                       <  2:134159/149‑1 (MQ=255)
                                                   tagcttattgttgtcgctgatgatatttaacacattgaaaaataagtaaaatacttatgagtcagagagttgtgatttttgcccttacTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTtc                                                                                    <  1:238449‑M1/61‑1 (MQ=255)
                                                    gggaaaaCCATGATGCTATCGCCTTTTTCAATCAGATTGAAGTCAGCAATTGCTTTGCCCACGTTACGACGCAGTCGTTTTTTTAATTTGTTTAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTtct                                                                                   <  2:225990/145‑1 (MQ=255)
                                                      gCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGAAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTctct                                                                                 >  2:272538/1‑149 (MQ=255)
                                                                   tgagcGCCTTCTTCAATCATATTGAAGTCAGCGATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCTGGTTTTATTGTTCTTTCTTTGTGATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATggg                                                                    <  2:251274/145‑1 (MQ=255)
                                                                   cGATCGCTTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGGACTAAATggg                                                                    <  2:211159/149‑1 (MQ=255)
                                                                   cGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATggg                                                                    <  1:9749/149‑1 (MQ=255)
                                                                   cGATCGCCTTCTTCAATCATATTGAAGGCAGCAATGGCTTCGACCACAACAAAACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATggg                                                                    <  1:10059/149‑1 (MQ=255)
                                                                        gCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCaaa                                                               <  1:340660/149‑1 (MQ=255)
                                                                            tttaacacattgaaaaataagtaaaatacttatgagtcagagagttgtgatttttgcccttacTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTgg                                                           <  2:20569‑M1/86‑1 (MQ=255)
                                                                               taaatattatTTAAGTCATTAATGGTTTCGCCCCCGTTACGTCGCTGCTGTTTTTGTTATTTGTTCAGGTTGTATTGTTCTTTCTTTTTAATTTTTTGTTTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTggggg                                                        <  2:247996/142‑1 (MQ=255)
                                                                                 aaTCATATTGAAGTCAGCAATGACTCCGTCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTGCTAAATGGGGCAGATTTGGGCCa                                                      >  2:278573/1‑149 (MQ=255)
                                                                                                      tacttatgagtcagagagttgtgatttttgcccttacTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAAcc                                 >  2:238449‑M1/38‑149 (MQ=255)
                                                                                                             cccACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCCTTGTAATTTGTGCATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTCACAACTCCCATAACCCCGCAtt                          >  2:18538/1‑149 (MQ=255)
                                                                                                                                gTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCCTTCAACATGGCTATCTATtcg       >  1:89796/1‑149 (MQ=255)
                                                                                                                                gTTTTTGTAATTTGTTCAGGGTGGATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTtcg       >  2:183273/1‑149 (MQ=255)
                                                                                                                                    ttGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTcgtcgt   >  2:203643/1‑149 (MQ=255)
                                                                                                                                     tGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCCACTACCCTATCCCCGCATTCAACATGGCTATCTGTTCGTCGtt  <  2:252608/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
CTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTT  >  NC_000913_3_vch_pgi/1411760‑1412041

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