Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 168,628 IS2 (+) +5 bp coding (1145‑1149/2244 nt) fhuA → ferrichrome outer membrane transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 168628 =3 (0.120)20 (0.790) 18/256 0.5 86.5% coding (1145/2244 nt) fhuA ferrichrome outer membrane transporter
?NC_000913 1467910 = NA (NA)noncoding (1331/1331 nt) IS2 repeat region
* ? NC_000913 = 1686323 (0.120)26 (1.030) 24/256 0.2 89.3% coding (1149/2244 nt) fhuA ferrichrome outer membrane transporter
?NC_000913 381260 = NA (NA)noncoding (1/1331 nt) IS2 repeat region

TGCTGAAAACAAAACCTCGCAAAACAGCGTTTATGGTTACGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCA  >  NC_000913/168488‑168776
                                                                                                                                                |                                                                                                                                                
tGCTGAAAACAAAACCTCGCAAAACAGCGTTTATGGTTACGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAg                                                                                                                                               <  6:24123/148‑1 (MQ=255)
    gAAAACAAAACCTCGCAAAACAGCGTTTATGGTTACGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactgg                                                                                                                                          >  1:132725‑M1/1‑141 (MQ=255)
    cagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaATGAGAAGCTGc                                                                                                                                           >  5:49841‑M1/137‑148 (MQ=255)
            aacCTCGCAAAACAGCGTTTATGGTTACGGCGTCGGCTCCGATCCGGCGAATGCTTACAGCAAACAGGGTGCGGCATTAGCGCCAGCGGATACAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctga                                                                                                                                  <  7:73049‑M1/149‑17 (MQ=255)
               gagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaATGAGAAGCTGCAAAACTTCTcc                                                                                                                                <  8:4694‑M1/23‑1 (MQ=255)
               gagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaATGAGAAGCTGCAAAACTTCTcc                                                                                                                                >  7:4694‑M1/126‑148 (MQ=255)
                                gcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatggggtaagtgataacagatgtctggaaatataggggcaaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTg                                                                                                                 <  2:128057‑M1/38‑1 (MQ=255)
                                 tGGTTACGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCag                                                                                                              <  5:76641/148‑1 (MQ=255)
                                       cGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctcca                                                                                                                           >  6:114297‑M1/1‑106 (MQ=255)
                                       cGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctcca                                                                                                                           <  5:114297‑M1/129‑24 (MQ=255)
                                           gTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataag                                                                                                   <  2:132725‑M1/149‑48 (MQ=255)
                                                cTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaat                                                                                                      >  4:1234‑M1/1‑97 (MQ=255)
                                                cTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaat                                                                                                      <  3:1234‑M1/141‑45 (MQ=255)
                                                 tCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatag                                                                                             >  4:66466‑M1/1‑96 (MQ=255)
                                                        cGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaat                                                                                      <  3:66466‑M1/149‑61 (MQ=255)
                                                                  ttACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttt                                                                             >  7:62459‑M1/1‑79 (MQ=255)
                                                                     cAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtctta                                                                                        >  5:75005‑M1/1‑76 (MQ=255)
                                                                     cAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtctta                                                                                        <  6:75005‑M1/134‑59 (MQ=255)
                                                                        cAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagacta                                                                      >  4:35452‑M1/1‑73 (MQ=255)
                                                                           ctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCaa                                                                                                         <  4:33412‑M1/46‑1 (MQ=255)
                                                                           ctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCaa                                                                                                         >  3:33412‑M1/66‑111 (MQ=255)
                                                                           aCAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagt                                                                    >  1:52994‑M1/1‑70 (MQ=255)
                                                                              gtgtGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcatt                                                                <  1:132943‑M1/149‑83 (MQ=255)
                                                                                 tgCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattgga                                                             <  3:35452‑M1/149‑86 (MQ=255)
                                                                                  gCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttatt                                                                                                       >  3:77521‑M1/1‑63 (MQ=255)
                                                                                  gCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttatt                                                                                                       <  4:77521‑M1/106‑44 (MQ=255)
                                                                                   cGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaat                                                                                      >  1:107050‑M1/1‑62 (MQ=255)
                                                                                   cGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaat                                                                                      <  2:107050‑M1/122‑61 (MQ=255)
                                                                                       aTTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacag                                                       <  2:52994‑M1/149‑92 (MQ=255)
                                                                                                   ggATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGACGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatga                                                   >  4:39290‑M1/1‑46 (MQ=255)
                                                                                                   ggATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGACGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatga                                                   <  3:39290‑M1/141‑96 (MQ=255)
                                                                                                      agtgataacagatgtctggaaatataggggcaaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAAt                                         >  1:128057‑M1/39‑148 (MQ=255)
                                                                                                       aaaGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttag                                       >  1:36613‑M1/1‑42 (MQ=255)
                                                                                                            aacagatgtctggaaatataggggcaaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGttt                                                                                                     >  2:129417‑M1/33‑82 (MQ=255)
                                                                                                            aacagatgtctggaaatataggggcaaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGttt                                                                                                     <  1:129417‑M1/50‑1 (MQ=255)
                                                                                                             cATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCAGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCaa                                 >  4:49956/1‑149 (MQ=255)
                                                                                                                      tggaaatataggggcaaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggtt                        <  4:77787‑M1/127‑1 (MQ=255)
                                                                                                                       ggaaatataggggcaaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCagag                                                                                                            <  8:7774‑M1/43‑1 (MQ=255)
                                                                                                                                   gcaaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGa                                                                                          >  4:23470‑M1/10‑70 (MQ=255)
                                                                                                                                   gcaaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGa                                                                                          <  3:23470‑M1/61‑1 (MQ=255)
                                                                                                                                     aaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggtttg                      >  7:46564‑M1/8‑136 (MQ=255)
                                                                                                                                     aaatccaATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggtttg                      <  8:46564‑M1/129‑1 (MQ=255)
                                                                                                                                      ttGATGATGAGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacaca        <  8:62459‑M1/149‑139 (MQ=255)
                                                                                                                                           aATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTTTGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTacgac              >  2:12933‑M1/2‑138 (MQ=255)
                                                                                                                                           aATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTTTGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTacgac              <  1:12933‑M1/137‑1 (MQ=255)
                                                                                                                                           aATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggtttggt                    >  4:48079‑M1/2‑132 (MQ=255)
                                                                                                                                           aATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggtttggt                    <  3:48079‑M1/131‑1 (MQ=255)
                                                                                                                                            atgaGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCa  <  3:77787/149‑1 (MQ=255)
                                                                                                                                               aGtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatacttggttttagactagtcattggagaacagatgattgatgtcttagggccggagaaaagcagaatgcgtacaacacaggaa    >  7:27284‑M1/1‑2 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TGCTGAAAACAAAACCTCGCAAAACAGCGTTTATGGTTACGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCA  >  NC_000913/168488‑168776

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