Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	517,047	G→T	Q155H (CAG→CAT)	fetB →	putative iron ABC exporter membrane subunit FetB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	517,047	0	G	T	100.0%	161.9 / NA	46	Q155H (CAG→CAT)	fetB	putative iron ABC exporter membrane subunit FetB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (23/23);  total (23/23)

CTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCAAAAACGGTTGGCTTAGTGAGTTTACCAG  >  NC_000913/516916‑517177                                                                                                                                    |                                                                                                                                   ctctcAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCTATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGga                                                                                                                            <  1:461121/140‑1 (MQ=255) ctctcAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGga                                                                                                                            <  1:678870/140‑1 (MQ=255)   ctcAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGgaga                                                                                                                          >  1:159732/1‑140 (MQ=255)      aGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAAc                                                                                                                       >  2:762393/1‑140 (MQ=255)       gggTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACt                                                                                                                      <  2:481343/140‑1 (MQ=255)       gggTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACt                                                                                                                      <  2:653622/140‑1 (MQ=255)       gggTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACt                                                                                                                      <  1:253271/140‑1 (MQ=255)       gggTCGATTGAATTTATCCCGAGGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACt                                                                                                                      <  1:904544/140‑1 (MQ=255)        ggTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTg                                                                                                                     >  2:748634/1‑140 (MQ=255)             aTTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCt                                                                                                                >  2:550939/1‑140 (MQ=255)                  aTTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTg                                                                                                           >  2:33931/1‑140 (MQ=255)                  aTTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTg                                                                                                           <  2:325085/140‑1 (MQ=255)                  aTTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTg                                                                                                           >  2:183741/1‑140 (MQ=255)                   tttATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGc                                                                                                          >  2:103892/1‑140 (MQ=255)                      aTCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGAc                                                                                                       >  2:316564/1‑140 (MQ=255)                       tCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACg                                                                                                      >  2:568314/1‑140 (MQ=255)                       tCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACg                                                                                                      <  2:159732/140‑1 (MQ=255)                           gATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGa                                                                                                  >  2:300897/1‑140 (MQ=255)                           gATGCAGGTGATCCCTATCGCCGGGACGAATGCCGGTAATCCCATCGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGTAACATATTCAGGAGAAACTGAGTCTTGGTCCGACGCCgc                                                                                                  >  2:201090/1‑139 (MQ=255)                                       cccTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAg                                                                                      <  1:183741/140‑1 (MQ=255)                                              gCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATAtt                                                                               >  2:849844/1‑140 (MQ=255)                                                  ggATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGAtt                                                                           <  2:604343/140‑1 (MQ=255)                                                    atgatTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGatat                                                                                <  1:691022/133‑1 (MQ=255)                                                    atgatTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGatat                                                                                >  2:691022/1‑133 (MQ=255)                                                    atgatTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTcg                                                                         >  1:617939/1‑140 (MQ=255)                                                       atTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGa                                                                      <  1:103892/140‑1 (MQ=255)                                                                  cGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAATGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTcgcg                                                           >  2:925945/1‑140 (MQ=255)                                                                   gCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTcgcgc                                                          <  2:309829/140‑1 (MQ=255)                                                                   gCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTcgcgc                                                          <  1:568314/140‑1 (MQ=255)                                                                   gCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTcgcgc                                                          <  2:617939/140‑1 (MQ=255)                                                                   gCCATGGTAGCGGTAGGGTTGTGTTACAACAATTTAGGGCAATGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTcgcgc                                                          <  1:925945/140‑1 (MQ=255)                                                                      aTGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGc                                                       <  2:647404/140‑1 (MQ=255)                                                                             cGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAAtt                                                >  1:340312/1‑140 (MQ=255)                                                                               gTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTcc                                              >  2:680053/1‑140 (MQ=255)                                                                                    gTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACgg                                         >  2:501920/1‑140 (MQ=255)                                                                                       gtgtTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCg                                      >  2:376074/1‑140 (MQ=255)                                                                                             caacaaTTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAg                                >  1:432888/1‑140 (MQ=255)                                                                                                  aTTTAGGGCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCaaaa                           >  2:901379/1‑140 (MQ=255)                                                                                                         gCAACGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCAAAAACGGTTg                    <  2:432888/140‑1 (MQ=255)                                                                                                           aaCGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCAAAAACGGTTGGc                  >  2:513759/1‑140 (MQ=255)                                                                                                             cGGGTCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCAAAAACGGTTGGCtt                >  1:467008/1‑140 (MQ=255)                                                                                                                 tCATTAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCAAAAACGGTTGGCTTAGTg            <  2:110005/140‑1 (MQ=255)                                                                                                                    ttAGCGAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCAAAAACGGTTGGCTTAGTGAGt         <  1:316564/140‑1 (MQ=255)                                                                                                                         gAACAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCAAAAACGGTTGGCTTAGTGAGTTTAcc    <  1:501920/140‑1 (MQ=255)                                                                                                                          aaCAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCAAAAACGGTTGGCTTAGTGAGTTTACCa   >  1:947816/1‑140 (MQ=255)                                                                                                                          aaCAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCAAAAACGGTTGGCTTAGTGAGTTTACCa   <  1:849844/140‑1 (MQ=255)                                                                                                                           aCAGCAACATATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCAAAAACGGTTGGCTTAGTGAGTTTACCAg  <  1:550939/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                   CTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAGCGATATTGATTCGCGACAGTATTCGCGCGGCTTTAATTCCGACGGTCGATTCAGCAAAAACGGTTGGCTTAGTGAGTTTACCAG  >  NC_000913/516916‑517177

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