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

ATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATT  >  NC_000913/516532‑516726                                                                                                    |                                                                                                aTAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGc                                                                                                 <  1:1091883/100‑1 (MQ=255)  tAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCt                                                                                                >  1:1764313‑M2/1‑99 (MQ=255)   aaaGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtg                                                                                               >  1:1524808‑M2/1‑98 (MQ=255)   aaaGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtg                                                                                               >  2:1600815‑M2/1‑98 (MQ=255)     aGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggt                                                                                             >  1:1487650‑M2/1‑96 (MQ=255)     aGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggt                                                                                             <  2:1115333‑M2/100‑5 (MQ=255)      gTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtg                                                                                            >  1:1115323‑M2/1‑95 (MQ=255)      gTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggt                                                                                             <  1:541309‑M2/99‑5 (MQ=255)       tGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtgg                                                                                           >  1:105523‑M2/1‑94 (MQ=255)        gATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACTAATCATTAGCACTGGCATTAATGCtg                                                                                               <  1:56473‑M2/95‑3 (MQ=255)           tACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggc                                                                                       <  1:1451240‑M2/100‑11 (MQ=255)            acacTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAATTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtgg                                                                                        <  1:608154‑M2/98‑10 (MQ=255)            acacTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggca                                                                                      <  2:490209‑M2/100‑12 (MQ=255)               cTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatc                                                                                   <  2:820043‑M2/100‑15 (MQ=255)               cTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatc                                                                                   >  1:302146‑M2/1‑86 (MQ=255)               cTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatc                                                                                   >  2:296584‑M2/1‑86 (MQ=255)               cTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGAACTGGCATTAATGCtggtggtggcaa                                                                                     <  1:1548983‑M2/98‑13 (MQ=255)                 gCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggcggtggcaatctt                                                                                 <  1:1680754‑M2/100‑17 (MQ=255)                  cAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatctta                                                                                >  1:260710‑M2/1‑83 (MQ=255)                     gtgcaTGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaat                                                                              >  1:882199‑M2/3‑80 (MQ=255)                     ccGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaatt                                                                             >  1:554650‑M2/1‑80 (MQ=255)                      cGCATGCCGGATAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaatta                                                                            >  2:90612‑M2/1‑79 (MQ=255)                      cGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaatta                                                                            <  2:1308883‑M2/100‑22 (MQ=255)                         aTGCCTGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagc                                                                          >  2:856619‑M2/1‑76 (MQ=255)                         aTGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagcc                                                                         >  1:1257052‑M2/1‑76 (MQ=255)                              ggAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaa                                                                    <  2:716689‑M2/100‑30 (MQ=255)                              ggAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaa                                                                    <  2:848881‑M2/100‑30 (MQ=255)                                         ggAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggc                                                         <  1:128422‑M2/100‑41 (MQ=255)                                             gCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtgttggtggcaatcttaattagccataaagaaaaactgacgct                                                      <  2:642246‑M2/99‑44 (MQ=255)                                             gCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctg                                                     <  1:227074‑M2/100‑45 (MQ=255)                                                 gCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggaga                                                 >  2:938875‑M2/1‑52 (MQ=255)                                                       aaCTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcactggagaaagata                                           <  2:207231‑M2/100‑55 (MQ=255)                                                        aCTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatat                                          >  2:1765630‑M2/1‑45 (MQ=255)                                                         cTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatatt                                         <  2:766385‑M2/100‑57 (MQ=255)                                                           cGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattct                                       <  1:1726101‑M2/100‑59 (MQ=255)                                                            gCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctc                                      <  1:1443090‑M2/100‑60 (MQ=255)                                                              ataataTTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctg                                    >  2:227075‑M2/1‑39 (MQ=255)                                                              ataataTTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctct                                     >  2:1340778‑M2/1‑39 (MQ=255)                                                              acaatattACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctg                                    >  1:1101370‑M2/3‑39 (MQ=255)                                                               taataTTACTAACGAATCCTTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctgg                                   >  2:1806890‑M2/1‑38 (MQ=255)                                                               taataTTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctgg                                   >  2:1778946‑M2/1‑38 (MQ=255)                                                                 atatTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggag                                 <  1:477860‑M2/100‑65 (MQ=255)                                                                   atTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcg                               <  2:779704‑M2/100‑67 (MQ=255)                                                                   atTACTAACGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcg                               >  1:1238900‑M2/1‑34 (MQ=255)                                                                         aaCGAATCATTAGCACTGGCATTAATGCtggtggtggctatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtc                             <  1:1117111‑M2/96‑69 (MQ=255)                                                                         aaCGAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggc                         <  2:1714679‑M2/100‑73 (MQ=255)                                                                            gAATCATTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgag                      <  1:1472767‑M2/100‑76 (MQ=255)                                                                               tCATTAGCACTGGTATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaatatattctctggagcgtcgggcgagc                     >  1:656552‑M2/1‑22 (MQ=255)                                                                                tatTAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgat                  <  2:1285122‑M2/99‑80 (MQ=255)                                                                                  ttAGCACTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgataa                <  1:729988‑M2/100‑82 (MQ=255)                                                                                   tAGCACTGGCATTAATGCtggtggtggcaatcttaaatagccataaagaaaaactggcgctggagaaagatattctctggaacgtcgggcgagcgataa                <  1:807968‑M2/99‑82 (MQ=255)                                                                                       aCTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgataattcag           >  2:814092‑M2/1‑14 (MQ=255)                                                                                        cTGGCATTAATGCtggtggtggcaatcttacttagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgataattcagc          <  2:1476121‑M2/100‑88 (MQ=255)                                                                                        cTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgataattcagc          >  2:1303228‑M2/1‑13 (MQ=255)                                                                                        cTGGCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgataattcagc          <  1:1190151‑M2/100‑88 (MQ=255)                                                                                          ggCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgataattcagctg        >  2:37454‑M2/1‑11 (MQ=255)                                                                                          ggCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgataattcagctg        >  1:490206‑M2/1‑11 (MQ=255)                                                                                          ggCATTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgacgggcgagcgataattca            <  2:1304535‑M2/96‑86 (MQ=255)                                                                                           gCATTAATGAgtgtggtggaaatcttaattagcaataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgataattcagct         >  1:663066‑M2/1‑10 (MQ=255)                                                                                             aTTAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgataattcagctgat      >  2:1641573‑M2/1‑8 (MQ=255)                                                                                               tAATGCtggttgtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagtgataattcagctgattat   >  1:848875‑M2/1‑6 (MQ=255)                                                                                               tAATGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgataattcagctgattat   <  1:370820‑M2/100‑95 (MQ=255)                                                                                                aaTGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcgataattcagctgattatt  <  1:1610979‑M2/100‑96 (MQ=255)                                                                                                aaTGCtggtggtggcaatcttaattagccataaagaaaaactggcgctggagaaagatattctctggagcgtcgggcgagcaataattcagctgattatt  <  1:10175‑M2/100‑96 (MQ=255)                                                                                                    |                                                                                                ATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATT  >  NC_000913/516532‑516726

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑