Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	169,645	A→G	Y721C (TAC→TGC)	fhuA →	ferrichrome outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	169,645	0	A	G	100.0%	185.1 / NA	52	Y721C (TAC→TGC)	fhuA	ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/28);  total (24/28)

TGGTGATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTG  >  NC_000913/169511‑169775                                                                                                                                       |                                                                                                                                   tGGTGATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTc                                                                                                                                >  1:1257223/1‑139 (MQ=255)       tCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAgc                                                                                                                          <  2:1267230/139‑1 (MQ=255)          ggCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAgctgc                                                                                                                       <  1:50880/139‑1 (MQ=255)              aaCTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                   <  2:150337/139‑1 (MQ=255)                         aGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                        <  2:1257223/139‑1 (MQ=255)                           tGGGAAGTTATACGGTCATGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACTTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgctg                                                                                                      >  1:782920/1‑139 (MQ=255)                                         gTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAg                                                                                        <  1:36066/139‑1 (MQ=255)                                          tCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGa                                                                                       >  1:9720/1‑139 (MQ=255)                                            gTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                        <  2:1103896/120‑1 (MQ=255)                                            gTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                        >  1:1103896/1‑120 (MQ=255)                                            gTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAAc                                                                                     >  1:916322/1‑139 (MQ=255)                                              ggATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGt                                                                                   <  2:1172444/139‑1 (MQ=255)                                                 tGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAg                                                                                <  2:916322/139‑1 (MQ=255)                                                      tAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCt                                                                                                 <  2:808173/117‑1 (MQ=255)                                                      tAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCt                                                                                                 >  1:808173/1‑117 (MQ=255)                                                              tATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccg                                                                   <  1:730946/139‑1 (MQ=255)                                                                      ggCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTc                                                           >  2:396855/1‑139 (MQ=255)                                                                       gcgcGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTcc                                                          >  1:129773/1‑139 (MQ=255)                                                                              tCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGttt                                                      >  2:695806/1‑136 (MQ=255)                                                                              tCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGttt                                                      <  1:695806/136‑1 (MQ=255)                                                                                        tggctCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCttt                                         <  2:129773/139‑1 (MQ=255)                                                                                                aCGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCAc                                 >  1:1042940/1‑139 (MQ=255)                                                                                                 cGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCtt                                                                                                   <  1:113495/72‑1 (MQ=255)                                                                                                 cGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCtt                                                                                                   >  2:113495/1‑72 (MQ=255)                                                                                                     gcgcTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                   >  1:867659/1‑52 (MQ=255)                                                                                                     gcgcTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                   <  2:867659/52‑1 (MQ=255)                                                                                                          gCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt                       <  2:782920/139‑1 (MQ=255)                                                                                                                tgttaacCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGgcgc                                                                                          <  1:1169750/62‑1 (MQ=255)                                                                                                                tgttaacCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGgcgc                                                                                          >  2:1169750/5‑66 (MQ=255)                                                                                                                tAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                   <  1:204499/41‑1 (MQ=255)                                                                                                                tAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                   >  2:204499/1‑41 (MQ=255)                                                                                                                tAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                        >  1:559444/1‑52 (MQ=255)                                                                                                                tAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                        <  2:559444/52‑1 (MQ=255)                                                                                                                tAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCAt                 <  2:9720/139‑1 (MQ=255)                                                                                                                   caacCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGTGCACGGATTTCCGTGCCCATTTc              >  2:276687/1‑139 (MQ=255)                                                                                                                     acCTGTTCGATCGTGAATGCGTCGCCAGCTGCttt                                                                                                                    >  1:384843/1‑35 (MQ=255)                                                                                                                     acCTGTTCGATCGTGAATGCGTCGCCAGCTGCttt                                                                                                                    <  2:384843/35‑1 (MQ=255)                                                                                                                     acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                   <  1:1265242/36‑1 (MQ=255)                                                                                                                     acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                   >  2:1265242/1‑36 (MQ=255)                                                                                                                     acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                    >  1:145034/1‑131 (MQ=255)                                                                                                                     acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                    <  2:145034/131‑1 (MQ=255)                                                                                                                         gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                        <  1:435874/43‑1 (MQ=255)                                                                                                                         gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                        >  2:435874/1‑43 (MQ=255)                                                                                                                         gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGt        >  1:109467/1‑139 (MQ=255)                                                                                                                          ttCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTg                                                                                                >  1:1065628/1‑50 (MQ=255)                                                                                                                          ttCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTg                                                                                                <  2:1065628/50‑1 (MQ=255)                                                                                                                            cGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTgg     <  2:1042940/139‑1 (MQ=255)                                                                                                                              aTCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                    <  2:1117710/122‑1 (MQ=255)                                                                                                                              aTCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                    >  1:1117710/1‑122 (MQ=255)                                                                                                                               tCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTg  <  1:396855/139‑1 (MQ=255)                                                                                                                                 gTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccg                                                                   >  2:407338/1‑72 (MQ=255)                                                                                                                                 gTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTgcaaccg                                                                   <  1:407338/72‑1 (MQ=255)                                                                                                                                       |                                                                                                                                   TGGTGATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTG  >  NC_000913/169511‑169775

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