Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 167,520 IS2 (–) +5 bp coding (37‑41/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 167520 =0 (0.000)35 (0.870) 34/252 0.2 100% coding (37/2244 nt) fhuA ferrichrome outer membrane transporter
?NC_000913 381260 = NA (NA)noncoding (1/1331 nt) IS2 repeat region
* ? NC_000913 = 1675240 (0.000)40 (1.000) 33/252 0.2 100% coding (41/2244 nt) fhuA ferrichrome outer membrane transporter
?NC_000913 = 382590 NA (NA)noncoding (1331/1331 nt) IS2 repeat region

AGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAAC  >  NC_000913/382446‑382737
                                                                                                                                                |                                                                                                                                                   
aGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTa                                                                                                                                                 >  2:79241/1‑149 (MQ=31)
aGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTa                                                                                                                                                 >  2:82208/1‑149 (MQ=31)
 gAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTAc                                                                                                                                                <  5:8745/149‑1 (MQ=32)
          agagGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAgg                                                                                                                                        >  5:56485/1‑148 (MQ=38)
                       cATTATAACGAATGGCATCCGCATAGTGAGCTGGGTTATCTCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAAAATAGGGGCAAATCCAATAACCTATAAGGCAGTTTGCGCCAAt                                                                                                                          >  3:9408/1‑149 (MQ=14)
                        aTTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                                                          <  4:92683/148‑1 (MQ=255)
                        aTTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATc                                                                                                                         <  8:27221/149‑1 (MQ=255)
                            taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                                                               >  6:98970/1‑139 (MQ=255)
                            taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                                                               <  5:98970/139‑1 (MQ=255)
                               cGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCg                                                                                                                  >  8:3007/1‑149 (MQ=255)
                                                gTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGc                                                                                                  <  7:3007/148‑1 (MQ=255)
                                                 tGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAc                                                                                                <  8:149698/149‑1 (MQ=255)
                                                       gggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGa                                                                                           >  7:144757/1‑148 (MQ=255)
                                                         ggtATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaa                                                                                        <  6:56485/147‑1 (MQ=255)
                                                         gTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaa                                                                                        <  1:79241/149‑1 (MQ=255)
                                                         gTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaa                                                                                        <  1:82208/149‑1 (MQ=255)
                                                          ttATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACaaa                                                                                       <  8:144757/149‑1 (MQ=255)
                                                               gctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACgg                                                                                  <  2:143144/149‑1 (MQ=255)
                                                                 tcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACg                                                                                   <  5:182532/146‑1 (MQ=255)
                                                                 tcgcCACGGGAAAATCTGCGGCAGCGGGCTAGTAATGGGTTAAGTGATAACAGAAGACTCGAAAAATACGCGCAAATCCAAATACCTATCAGCCAGATTCCCCCAATCACAACCGGACCCACAAAAACGCCACCAAAACaaaaaac                                                                                   >  6:182532/1‑142 (MQ=12)
                                                                        gggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAgg                                                                                                                                        >  1:140208/1‑86 (MQ=255)
                                                                        gggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAgg                                                                                                                                        <  2:140208/86‑1 (MQ=255)
                                                                          gAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTgcgc                                                                                                                              <  8:11265/94‑1 (MQ=255)
                                                                          gAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTgcgc                                                                                                                              >  7:11265/1‑94 (MQ=255)
                                                                           aaTATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCa                                                                                                                                          <  5:140444/81‑1 (MQ=255)
                                                                           aaTATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCa                                                                                                                                          >  6:140444/1‑81 (MQ=255)
                                                                                                   aTGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAAc                                                                                                     <  3:45607/94‑1 (MQ=255)
                                                                                                   aTGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAAc                                                                                                     >  4:45607/1‑94 (MQ=255)
                                                                                                       gTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAAc                                                                                                     >  3:169493/1‑90 (MQ=255)
                                                                                                       gTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAAc                                                                                                     <  4:169493/90‑1 (MQ=255)
                                                                                                           aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtc                                               <  6:83792/140‑1 (MQ=255)
                                                                                                           aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtc                                               >  5:83792/1‑140 (MQ=255)
                                                                                                                 aaCAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTg                                >  4:3803/1‑149 (MQ=255)
                                                                                                                                     ggggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGaa                                                                                                         >  3:48269/1‑56 (MQ=255)
                                                                                                                                     ggggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGaa                                                                                                         <  4:48269/56‑1 (MQ=255)
                                                                                                                                     ggggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATAttt             <  3:3803/148‑1 (MQ=255)
                                                                                                                                        gCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAGCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAAt                       <  4:17071/135‑1 (MQ=255)
                                                                                                                                              ccAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaa                                                                  >  7:116055/1‑86 (MQ=255)
                                                                                                                                              ccAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaa                                                                  <  8:116055/86‑1 (MQ=255)
                                                                                                                                                aaTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAAc  <  6:28882/148‑1 (MQ=255)
                                                                                                                                                aaTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAAc  >  5:28882/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAAC  >  NC_000913/382446‑382737

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.