Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	168,708	A→G	T409A (ACC→GCC)	fhuA →	ferrichrome outer membrane transporter/phage receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	168,708	0	A	G	100.0%	165.2 / NA	48	T409A (ACC→GCC)	fhuA	ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/23);  total (25/23)

ATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCG  >  NC_000913/168575‑168843                                                                                                                                      |                                                                                                                                        aTTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGgtg                                                                                                                                   >  2:760641/1‑140 (MQ=255)   tAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTc                                                                                                                                 >  2:218011/1‑140 (MQ=255)       gcCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACt                                                                                                                             <  2:192390/140‑1 (MQ=255)            cGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTtatg                                                                                                                        >  1:211890/1‑140 (MQ=255)               aTAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTtatgcgt                                                                                                                     >  1:549965/1‑140 (MQ=255)                      ccATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCATAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTAtgcgtatgcgta                                                                                                              >  1:418355/1‑140 (MQ=255)                             cTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACAt                                                                                                       >  2:662459/1‑140 (MQ=255)                                    gTAAATACGTTGTTCATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGcc                                                                                                <  2:510774/140‑1 (MQ=255)                                      aaaTACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtg                                                                                              >  2:489613/1‑140 (MQ=255)                                        aTACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggt                                                                                            <  2:211890/140‑1 (MQ=255)                                             tcgtTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggtttggt                                                                                       <  1:463320/140‑1 (MQ=255)                                             tcgtTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggtttggt                                                                                       <  1:218011/140‑1 (MQ=255)                                              cgtTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggtttggtt                                                                                      >  1:244518/1‑140 (MQ=255)                                               gtTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTa                                                                                     >  1:715617/1‑140 (MQ=255)                                                    tgatgaGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTacgacg                                                                                <  1:662459/140‑1 (MQ=255)                                                     gatggGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTacgacgt                                                                               <  2:466166/140‑2 (MQ=255)                                                        gagaAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCTATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGAct                                                                             >  1:366147/1‑139 (MQ=255)                                                        gagaAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCTATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGAct                                                                             <  2:366147/139‑1 (MQ=255)                                                        gagaAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGActc                                                                            >  2:788651/1‑140 (MQ=255)                                                           aaGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCtgt                                                                         >  2:534160/1‑140 (MQ=255)                                                           aaGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCtgt                                                                         >  1:224878/1‑140 (MQ=255)                                                             gCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTAcgacga                                                                               <  1:354995/132‑1 (MQ=255)                                                             gCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTAcgacga                                                                               <  1:743422/132‑1 (MQ=255)                                                             gCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTAcgacga                                                                               >  2:354995/1‑132 (MQ=255)                                                             gCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGc                                                                       <  1:275489/140‑1 (MQ=255)                                                             gCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGGTAcgacga                                                                               >  2:743422/1‑132 (MQ=255)                                                                   aaaCTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGc                                                                 <  2:244518/140‑1 (MQ=255)                                                                   aaaCTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGc                                                                 <  2:224878/140‑1 (MQ=255)                                                                    aaCTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCt                                                                <  1:788651/140‑1 (MQ=255)                                                                              ttGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGcc                                                                      <  2:496097/124‑1 (MQ=255)                                                                              ttGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGcc                                                                      >  1:496097/1‑124 (MQ=255)                                                                                 aTACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAAt                                                   <  2:549965/140‑1 (MQ=255)                                                                                   aCCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATcc                                                 >  2:405147/1‑140 (MQ=255)                                                                                         ttGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGaa                                           >  1:190943/1‑140 (MQ=255)                                                                                           gCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATa                                         >  2:723927/1‑140 (MQ=255)                                                                                           gCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATa                                         >  1:749124/1‑140 (MQ=255)                                                                                                          cACTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACAACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCa                          <  2:418355/140‑1 (MQ=255)                                                                                                           aCTGGCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCaa                         >  1:259070/1‑140 (MQ=255)                                                                                                               gCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTc                                 >  2:627949/1‑128 (MQ=255)                                                                                                               gCGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTc                                 <  1:627949/128‑1 (MQ=255)                                                                                                                cGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCa                    <  1:335953/140‑1 (MQ=255)                                                                                                                cGATATCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCa                    <  2:749124/140‑1 (MQ=255)                                                                                                                    atCGACCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGa                >  1:777647/1‑140 (MQ=255)                                                                                                                        aCCACACCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCg            >  2:441035/1‑140 (MQ=255)                                                                                                                           acacCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCa         <  2:275518/140‑1 (MQ=255)                                                                                                                           acacCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCa         <  2:190943/140‑1 (MQ=255)                                                                                                                             acCCTGCTGGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCaaa       >  1:604863/1‑140 (MQ=255)                                                                                                                                  gctgGCCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCTTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCg  <  1:534160/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        ATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAATCCGGTGAATACCGATTTCGACTTCAATGCCAAAGATCCGGCAAACTCCG  >  NC_000913/168575‑168843

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