Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	168,660	C→T	Q393* (CAG→TAG)	fhuA →	ferrichrome outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	168,660	0	C	T	100.0%	121.5 / NA	36	Q393* (CAG→TAG)	fhuA	ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/18);  total (18/18)

TTTATGGTTACGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAA  >  NC_000913/168517‑168793                                                                                                                                                |                                                                                                                                      tttATGGTTACGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTg                                                                                                                                  >  5:75496/1‑149 (MQ=255)           cGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAg                                                                                                                         <  6:75496/148‑1 (MQ=255)             gCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGttt                                                                                                                      <  5:103745/149‑1 (MQ=255)                       cGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGc                                                                                                             >  5:112767/1‑148 (MQ=255)                        gATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGa                                                                                                           >  4:86107/1‑149 (MQ=255)                              gCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGa                                                                                                     <  7:76249/149‑1 (MQ=255)                               cGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGAc                                                                                                    <  6:40874/149‑1 (MQ=255)                               cGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGAc                                                                                                    <  6:112767/149‑1 (MQ=255)                                              aaCAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGAcc                                                                                     >  8:5227/1‑149 (MQ=255)                                                            ttAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTtat                                                                       >  1:121976/1‑149 (MQ=255)                                                                            aaGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGAc                                                       <  7:2461/149‑1 (MQ=255)                                                                             aGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTaa                                                           <  3:67939/144‑1 (MQ=255)                                                                             aGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTaa                                                           >  4:67939/1‑144 (MQ=255)                                                                               gCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACAt                                                     >  3:63264/1‑148 (MQ=255)                                                                               gCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACAt                                                     >  1:77489/1‑148 (MQ=255)                                                                                  aTTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAAc                                                 >  1:72308/1‑149 (MQ=255)                                                                                   ttATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACg                                                >  1:51231/1‑149 (MQ=255)                                                                                       cTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCt                                             <  2:72308/148‑1 (MQ=255)                                                                                        tGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtgg                                           >  2:58971/1‑149 (MQ=255)                                                                                           cACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggttt                                        <  4:63264/149‑1 (MQ=255)                                                                                            aCGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggtttg                                       <  3:86107/149‑1 (MQ=255)                                                                                                 aaTACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTa                                   <  2:7768/148‑1 (MQ=255)                                                                                                    aCGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTacgac                               >  8:68244/1‑149 (MQ=255)                                                                                                         gtTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGActc                          >  7:20495/1‑149 (MQ=255)                                                                                                          ttGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggtttg                                       >  3:82101/1‑135 (MQ=255)                                                                                                          ttGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCtggtttg                                       <  4:82101/135‑1 (MQ=255)                                                                                                             atgatgaGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTtatgcgtatgcg                                                              >  6:71601/1‑109 (MQ=255)                                                                                                             atgatgaGAAGCTGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTtatgcgtatgcg                                                              <  5:71601/109‑1 (MQ=255)                                                                                                                         tGCAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAAt          >  6:71684/1‑149 (MQ=255)                                                                                                                           cAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAAt          <  6:123381/147‑1 (MQ=255)                                                                                                                           cAAAACTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCt        <  2:121976/149‑1 (MQ=255)                                                                                                                              aaCTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTa     <  8:20495/149‑1 (MQ=255)                                                                                                                               aCTTCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTAc    <  7:68244/149‑1 (MQ=255)                                                                                                                                 ttCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTtat                                                                       >  3:57315/1‑80 (MQ=255)                                                                                                                                 ttCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTtat                                                                       <  4:57315/80‑1 (MQ=255)                                                                                                                                 ttCTCCGTTGATACCTAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACaa  >  5:137510/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                      TTTATGGTTACGGCGTCTGCTCCGATCCGGCGAATGCTTACAGCAAACAGTGTGCGGCATTAGCGCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCGTATGCGTAATGACATCAACGCCTGGTTTGGTTACGACGACTCTGTGCCACTGCTCAATCTGTACAA  >  NC_000913/168517‑168793

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