Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	170,771	G→A	R66H (CGC→CAC)	fhuD →	iron(3+)‑hydroxamate import ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	170,771	0	G	A	100.0%	156.7 / NA	45	R66H (CGC→CAC)	fhuD	iron(3+)‑hydroxamate import ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/26);  total (19/26)

TTCTCCGTTGTTATGGCAGATGAATACCGCCCACGCGGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCGCCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAGGATATGGCCCTTCACCAGAAATGCTGGCTC  >  NC_000913/170628‑170917                                                                                                                                                |                                                                                                                                                   ttCTCCGTTGTTATGGCAGATGAATACCGCCCACGCGGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTgg                                                                                                                                             >  2:243392/1‑151 (MQ=255)   ctcCGTTGTTATGGCAGATGAATACCGCCCACGCGGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGGTACCATCCACTATCACCTGTggg                                                                                                                                            >  1:464628/1‑150 (MQ=255)       gttgttATGGCAGATGAATACCGCCCACGCGGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGc                                                                                                                                       <  2:673290/151‑1 (MQ=255)                     tGAATACCGCCCACGCGGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCCGCGAACCACCCTTGcc                                                                                                                         >  1:526141/1‑151 (MQ=255)                     tGAATACCGCCCACGCGGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCCTTGcc                                                                                                                         >  2:291741/1‑151 (MQ=255)                       aaTACCGCCCACGCGGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCgg                                                                                                                       <  2:464664/151‑1 (MQ=255)                          aCCGCCCACGCGGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGAc                                                                                                                     <  1:465532/150‑1 (MQ=255)                            cGCCCACGCGGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCa                                                                                                                  <  1:431980/151‑1 (MQ=255)                              cccACGCGGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGt                                                                                                                <  2:269756/151‑1 (MQ=255)                                 aCGCGGCGGTTATTGATCCCAATCGTATTGTGGCGATGGAGTTGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCAATGCCGGCCTCAGTGAt                                                                                                             <  1:480345/151‑1 (MQ=255)                                  cgcgGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATc                                                                                                            <  1:517481/151‑1 (MQ=255)                                    cggcggCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGa                                                                                                          <  2:16634/151‑1 (MQ=255)                                                aTCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTgc                                                                                               >  1:421982/1‑150 (MQ=255)                                                  cccAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGca                                                                                            <  2:576696/151‑1 (MQ=255)                                                   ccAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGcac                                                                                           <  2:377248/151‑1 (MQ=255)                                                     aaTCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTgcgc                                                                                             >  2:420446/1‑147 (MQ=255)                                                      aTCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGa                                                                                        <  1:733689/151‑1 (MQ=255)                                                        cGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGaac                                                                                      <  1:840660/151‑1 (MQ=255)                                                        cGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGGCAGCGAACCACCATTGCCTGACTCAGTGATCGACGTCGGGTTTCGCACAGaa                                                                                       >  2:611043/1‑150 (MQ=255)                                                                   gcTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGaa                                                                           <  2:211029/151‑1 (MQ=255)                                                                      ggAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAActg                                                                        >  2:676095/1‑151 (MQ=255)                                                                         gTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAActgctg                                                                     <  1:274680/151‑1 (MQ=255)                                                                         gTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAActgctg                                                                     <  1:655834/151‑1 (MQ=255)                                                                                  ggTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATg                                                            >  2:113306/1‑151 (MQ=255)                                                                                   gTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATg                                                            >  2:353088/1‑150 (MQ=255)                                                                                      gAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGaaa                                                         >  1:831363/1‑150 (MQ=255)                                                                                          tACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCAt                                                     <  1:619983/150‑1 (MQ=255)                                                                                                 ggCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTa                                               <  2:831421/149‑1 (MQ=255)                                                                                                       cGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTAtggtctg                                        <  2:839328/150‑1 (MQ=255)                                                                                                        ggCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTAtggtctggt                                      <  1:404314/151‑1 (MQ=255)                                                                                                           aTCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCgg                                   >  2:221533/1‑151 (MQ=255)                                                                                                           aTCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCgg                                   <  2:180295/151‑1 (MQ=255)                                                                                                           aTCGTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCgg                                   >  2:572485/1‑151 (MQ=255)                                                                                                              gTGCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAg                                <  2:293071/151‑1 (MQ=255)                                                                                                                gCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAgg                               >  1:85783/1‑150 (MQ=255)                                                                                                                gCCTTACGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAgg                               >  2:740995/1‑150 (MQ=255)                                                                                                                   ttACGGCCTGGCGGATAACATCAACCATTAACTTACGGGCCGGGAACCACCAATGCAGGACCCAGAGATCCACGACGGTTTGCGCAACGGACCCCAACTTGAAATGCTGACCGAACTGAACCCATCGTTTATGGTCTGGTCGGAAgg                               >  1:684128/1‑147 (MQ=255)                                                                                                                     atggCGTGTCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACGCAGTGATCGACGTCGGTTTGCGCACATAACCTGGTCTTGAACTGCTGGCCGAAATGAAACCATCGTTTATGGTCGGGTCGGCAGGatat                           <  1:564207/147‑1 (MQ=255)                                                                                                                     aCGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAGGATATg                          <  2:174538/150‑1 (MQ=255)                                                                                                                      cGGCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAGGATATgg                         >  1:502509/1‑150 (MQ=255)                                                                                                                       ggCGTGGCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAGGATATGGcc                       >  2:721688/1‑151 (MQ=255)                                                                                                                            ggCGGATACCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAGGATATGGCCCTTCa                  <  2:716016/151‑1 (MQ=255)                                                                                                                                   aCCATCAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAGGATATGGCCCTTCACCAGaa            >  2:120150/1‑150 (MQ=255)                                                                                                                                        cAACTATCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAGGATATGGCCCTTCACCAGAAATGCTg      <  1:238655/151‑1 (MQ=255)                                                                                                                                             aTCACCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAGGATATGGCCCTTCACCAGAAATGCTGGCTc  <  1:501916/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   TTCTCCGTTGTTATGGCAGATGAATACCGCCCACGCGGCGGCTATTGATCCCAATCGTATTGTGGCGCTGGAGTGGTTGCCGGTGGAATTACTGCTGGCGCTCGGCATCGTGCCTTACGGCGTGGCGGATACCATCAACTATCGCCTGTGGGTCAGCGAACCACCATTGCCGGACTCAGTGATCGACGTCGGTTTGCGCACAGAACCTAACCTTGAACTGCTGACCGAAATGAAACCATCGTTTATGGTCTGGTCGGCAGGATATGGCCCTTCACCAGAAATGCTGGCTC  >  NC_000913/170628‑170917

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