New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 253783716 (0.340)16 (0.340) 14/484 1.7 50.0% coding (89/735 nt) yfeS PF05406 family protein YfeS
?NC_000913 2996361 = NA (NA)noncoding (1331/1331 nt) IS2 repeat region

TATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACT  >  NC_000913/2996138‑2996602
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
tgtaGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATc                                                                                                                                                                                                                         >  2:471977/3‑251 (MQ=255)
      gTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTAt                                                                                                                                                                                                                   <  1:358303/251‑1 (MQ=255)
                      aCCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTAc                                                                                                                                                                                                    >  1:28723/1‑250 (MQ=255)
                                       tgtgATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTAt                                                                                                                                                                                                                   >  1:107533/1‑218 (MQ=255)
                                       tgtgATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTAt                                                                                                                                                                                                                   <  2:107533/218‑1 (MQ=255)
                                                                                     ttCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAtt                                                                                                                                     >  1:248637/1‑251 (MQ=255)
                                                                                      tCCAGCGGTATAGCTGGTAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACACGCCTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGAAAACAGCTCAATGGAAGTCAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCg                                                                                                                                   >  2:76788/1‑251 (MQ=255)
                                                                                      tCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTc                                                                                                                                    <  2:374244/250‑1 (MQ=255)
                                                                                          gCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTtat                                                                                                                               <  1:371930/251‑1 (MQ=255)
                                                                                                                                             aCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGtt                                                                           <  2:372720/251‑1 (MQ=255)
                                                                                                                                                                  ctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGt                                                      >  2:48956/1‑251 (MQ=255)
                                                                                                                                                                    cGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTc                                                     >  1:345510/1‑250 (MQ=255)
                                                                                                                                                                                 aaCCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGttt                                       >  2:107988/1‑251 (MQ=255)
                                                                                                                                                                                  aCCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGatat                                                                                                                                                                                <  1:422882/114‑1 (MQ=255)
                                                                                                                                                                                  aCCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGatat                                                                                                                                                                                >  2:422882/1‑114 (MQ=255)
                                                                                                                                                                                    cACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTc                                     >  1:311356/1‑250 (MQ=255)
                                                                                                                                                                                        gAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTGATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGACCTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCAGGATACTGATGTAGTCACGCTTTATCGTTTTCACGaa                                <  1:152939/251‑1 (MQ=255)
                                                                                                                                                                                         aGTAAACAGATCAATGGAAGTACATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCCGACATCTGTTATCACTTAACCCATTACCAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGCACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑AGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGa                                 >  2:44328/1‑249 (MQ=255)
                                                                                                                                                                                          gTAAACAGATGAATGGAAGTAAATTCCAATTGTTTATGGGATTT‑GCCCCTATAGTGCCCGCCAGCGGTTATCACTTAACACAATTCAAGCGCGCTGCCGCAGATATTCCCGTGTCGATCGATAACCCAGCGGACTATGCGGATGCCATTCGTTATGAGGGTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCTCG‑CTGTTTTGGTCATTATACTGATATAGTCACGCTGTATCGTTGTCACGaata                              <  2:277838/251‑3 (MQ=9)
                                                                                                                                                                                                aGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTg                        >  2:338075/1‑251 (MQ=255)
                                                                                                                                                                                                  aTCAATGGAATTAAATTCCAATTGTTTATTTGATTT‑GCCCCTATATTTCCAGACATCTTTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGGGCGATAACCCAGCTCACTATGCCGATGCGATTCGTTATAATGCTCGATCGCCTCTGCAATCTTCTTTGCTGCCGTTAACC‑CGTCTCGTTTGTGCATGATACGTATGTAGTCACGCGTTATCCTTTTCCCGAATGTCTCTGCt                      <  2:284764/251‑1 (MQ=14)
                                                                                                                                                                                                     aaTGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCTCTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTGTCACGAAGCTCTCTGCTAtt                   <  2:35172/251‑1 (MQ=255)
                                                                                                                                                                                                     aaTGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTAtt                   <  1:323239/251‑1 (MQ=40)
                                                                                                                                                                                                        ggAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCg                >  2:174933/1‑251 (MQ=39)
                                                                                                                                                                                                               aaTTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTActct         >  1:383426/1‑251 (MQ=35)
                                                                                                                                                                                                               aaTTCCAATTGGGTATTGGAGTT‑ACCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACAAAGCTCTCTGCTATTGCTTTActct         <  2:311030/251‑1 (MQ=11)
                                                                                                                                                                                                                aTTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGCCATCTGTTATCACTTAACCCATTACAAGACCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGCCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACtctc        >  2:333035/1‑251 (MQ=25)
                                                                                                                                                                                                                        tGTTTATTGGATTTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACt  <  2:248637/250‑1 (MQ=17)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
TATAGGGTAACAACGTTGAAAAACCACATGTGATGATATTGTGATAGTTAGCGATGCCGGGTTTGTATCGTGCCTTTTCATTTATTTCCAGCGGTATAGCTGATAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTT‑GCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACC‑CGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACT  >  NC_000913/2996138‑2996602

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