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: ‑