Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 516,641 (TGG)4→3 coding (59‑61/780 nt) fetB → iron export ABC transporter permease; peroxide resistance protein

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 516632 516634 3 64 [0] [1] 64 fetB iron export ABC transporter permease; peroxide resistance protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 516635 =1 (0.010)53 (0.780) 42/180 0.8 99.2% coding (53/780 nt) fetB iron export ABC transporter permease; peroxide resistance protein
?NC_000913 = 516631 0 (0.000)coding (49/780 nt) fetB iron export ABC transporter permease; peroxide resistance protein

AGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATT  >  NC_000913/516536‑516729
                                                                                                   |                                                                                              
taaagtgattacactgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgct                                                                                                >  1:1764313‑M1/100‑100 (MQ=255)
 aaagtgattacactgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgctg                                                                                               >  2:1600815‑M1/99‑100 (MQ=255)
 aaagtgattacactgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgctg                                                                                               >  1:1524808‑M1/99‑100 (MQ=255)
   agtgattacactgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgctggt                                                                                             <  2:1115333‑M1/4‑1 (MQ=255)
   agtgattacactgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgctggt                                                                                             >  1:1487650‑M1/97‑100 (MQ=255)
    gtgattacactgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgctggtg                                                                                            >  1:1115323‑M1/96‑100 (MQ=255)
    gtgattacactgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgctggt                                                                                             <  1:541309‑M1/4‑1 (MQ=255)
     tgattacactgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgctggtgg                                                                                           >  1:105523‑M1/95‑100 (MQ=255)
      gattacactgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaactaatcattagcactggcattaatgctg                                                                                               <  1:56473‑M1/2‑1 (MQ=255)
         tacactgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGc                                                                                       <  1:1451240‑M1/10‑1 (MQ=255)
          acactgcaaccgcatgccggagaaatgcaggaagcacgctatgaattcgcataatattactaacgaatcattagcactggcattaatgctggtggtgg                                                                                        <  1:608154‑M1/9‑1 (MQ=255)
          acactgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCa                                                                                      <  2:490209‑M1/11‑1 (MQ=255)
             ctgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATc                                                                                   >  2:296584‑M1/87‑100 (MQ=255)
             ctgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATc                                                                                   >  1:302146‑M1/87‑100 (MQ=255)
             ctgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATc                                                                                   <  2:820043‑M1/14‑1 (MQ=255)
             ctgcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagaactggcattaatgcTGGTGGTGGCaa                                                                                     <  1:1548983‑M1/12‑1 (MQ=255)
               gcaaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGCGGTGGCAATCtt                                                                                 <  1:1680754‑M1/16‑1 (MQ=255)
                caaccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTa                                                                                >  1:260710‑M1/84‑100 (MQ=255)
                   gtgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAAt                                                                              >  1:882199‑M1/81‑99 (MQ=255)
                   ccgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAAtt                                                                             >  1:554650‑M1/81‑100 (MQ=255)
                    cgcatgccggataaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTa                                                                            >  2:90612‑M1/80‑100 (MQ=255)
                    cgcatgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTa                                                                            <  2:1308883‑M1/21‑1 (MQ=255)
                       atgcctgagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGc                                                                          >  2:856619‑M1/77‑99 (MQ=255)
                       atgccggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGcc                                                                         >  1:1257052‑M1/77‑100 (MQ=255)
                            ggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATaaa                                                                    <  2:848881‑M1/29‑1 (MQ=255)
                            ggagaaatgcaggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATaaa                                                                    <  2:716689‑M1/29‑1 (MQ=255)
                                       ggaagcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGgc                                                         <  1:128422‑M1/40‑1 (MQ=255)
                                           gcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGTTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGACGCt                                                      <  2:642246‑M1/43‑1 (MQ=255)
                                           gcacgctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTg                                                     <  1:227074‑M1/44‑1 (MQ=255)
                                               gctatgaactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGgaga                                                 >  2:938875‑M1/53‑100 (MQ=255)
                                                     aactcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCACTGGAGAAAGata                                           <  2:207231‑M1/54‑1 (MQ=255)
                                                      actcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGatat                                          >  2:1765630‑M1/46‑100 (MQ=255)
                                                       ctcgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATAtt                                         <  2:766385‑M1/56‑1 (MQ=255)
                                                         cgcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATtct                                       <  1:1726101‑M1/58‑1 (MQ=255)
                                                          gcataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATtctc                                      <  1:1443090‑M1/59‑1 (MQ=255)
                                                            ataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTctct                                     >  2:1340778‑M1/40‑99 (MQ=255)
                                                            ataatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTg                                    >  2:227075‑M1/40‑100 (MQ=255)
                                                            acaatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTg                                    >  1:1101370‑M1/40‑100 (MQ=255)
                                                             taatattactaacgaatccttagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTgg                                   >  2:1806890‑M1/39‑100 (MQ=255)
                                                             taatattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTgg                                   >  2:1778946‑M1/39‑100 (MQ=255)
                                                               atattactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAg                                 <  1:477860‑M1/64‑1 (MQ=255)
                                                                 attactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCg                               >  1:1238900‑M1/35‑100 (MQ=255)
                                                                 attactaacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCg                               <  2:779704‑M1/66‑1 (MQ=255)
                                                                       aacgaatcattagcactggcattaatgcTGGTGGTGGCTATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTc                             <  1:1117111‑M1/68‑1 (MQ=255)
                                                                       aacgaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGgc                         <  2:1714679‑M1/72‑1 (MQ=255)
                                                                          gaatcattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGgcgag                      <  1:1472767‑M1/75‑1 (MQ=255)
                                                                             tcattagcactggtattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAATATATTCTCTGGAGCGTCGGgcgagc                     >  1:656552‑M1/23‑98 (MQ=255)
                                                                              tattagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGAt                  <  2:1285122‑M1/79‑1 (MQ=255)
                                                                                ttagcactggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATaa                <  1:729988‑M1/81‑1 (MQ=255)
                                                                                 tagcactggcattaatgcTGGTGGTGGCAATCTTAAATAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAACGTCGGGCGAGCGATaa                <  1:807968‑M1/81‑1 (MQ=255)
                                                                                     actggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAg           >  2:814092‑M1/15‑100 (MQ=255)
                                                                                      ctggcattaatgcTGGTGGTGGCAATCTTACTTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGc          <  2:1476121‑M1/87‑1 (MQ=255)
                                                                                      ctggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGc          <  1:1190151‑M1/87‑1 (MQ=255)
                                                                                      ctggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGc          >  2:1303228‑M1/14‑100 (MQ=255)
                                                                                        ggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTg        >  2:37454‑M1/12‑100 (MQ=255)
                                                                                        ggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTg        >  1:490206‑M1/12‑100 (MQ=255)
                                                                                        ggcattaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGACGGGCGAGCGATAATTCa            <  2:1304535‑M1/85‑1 (MQ=255)
                                                                                         gcattaatgaGTGTGGTGGAAATCTTAATTAGCAATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCt         >  1:663066‑M1/11‑98 (MQ=255)
                                                                                           attactcgtggtggtggCAATCTGAATGAGTCGTAAAGATAAACGGGTGCTGGAGAAAGGGATTCTCCGTAGCGTTGATCGAGCGATAATTCAGCTg        <  1:1329843/90‑1 (MQ=255)
                                                                                           attaatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGat      >  2:1641573‑M1/9‑99 (MQ=255)
                                                                                             taatgcTGGTTGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGattat   >  1:848875‑M1/7‑100 (MQ=255)
                                                                                             taatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGattat   <  1:370820‑M1/94‑1 (MQ=255)
                                                                                              aatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGattatt  <  1:1610979‑M1/95‑1 (MQ=255)
                                                                                              aatgcTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCAATAATTCAGCTGattatt  <  1:10175‑M1/95‑1 (MQ=255)
                                                                                                   |                                                                                              
AGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATT  >  NC_000913/516536‑516729

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