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%	183.5 / NA	52	Y721C (TAC→TGC)	fhuA	ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/24);  total (28/24)

GATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATG  >  NC_000913/169515‑169780                                                                                                                                   |                                                                                                                                             gATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCa                                                                                                                                      <  1:202110/139‑1 (MQ=255)  aTCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                                                     <  1:93595/139‑1 (MQ=255)    ccGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAgct                                                                                                                                   >  2:616898/1‑139 (MQ=255)          aaCTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTa                                                                                                                             <  2:435328/139‑1 (MQ=255)                         ggAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                                                     <  1:1087066/116‑1 (MQ=255)                         ggAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAg                                                                                                                                     >  2:1087066/1‑116 (MQ=255)                         ggAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgctgct                                                                                                              >  1:1653343/1‑139 (MQ=255)                                        gTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAAc                                                                                          >  2:1190704/1‑139 (MQ=255)                                           gATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGgcg                                                                                                     <  1:409794/130‑1 (MQ=255)                                           gATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGgcg                                                                                                     >  2:409794/1‑130 (MQ=255)                                                gTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAgc                                                                                                                                    >  2:412795/1‑94 (MQ=255)                                                gTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTc                                                                                  >  2:853068/1‑139 (MQ=255)                                                gTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAAACTGTTCGATCGTGAATGCGTCGCCAgc                                                                                                                                    <  1:412795/94‑1 (MQ=255)                                                    gTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTg                                                                              <  2:1301075/139‑1 (MQ=255)                                                     tACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTgc                                                                             <  1:805131/139‑1 (MQ=255)                                                      aCGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTgca                                                                            >  2:1616670/1‑139 (MQ=255)                                                       cGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTgcaa                                                                           <  2:1132925/139‑1 (MQ=255)                                                         ttATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTgcaacc                                                                         <  1:1578325/139‑1 (MQ=255)                                                          tATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAac                                                                                                                           <  1:1521144/93‑1 (MQ=255)                                                          tATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAac                                                                                                                           >  2:1521144/1‑93 (MQ=255)                                                            tGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACtt                                                                                                                       <  1:849803/95‑1 (MQ=255)                                                            tGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACtt                                                                                                                       >  2:849803/1‑95 (MQ=255)                                                              aTCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTgcaaccgcaac                                                                    >  2:500131/1‑139 (MQ=255)                                                                  ggCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                                  >  1:796503/1‑94 (MQ=255)                                                                  ggCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGgc                                                                                                                  <  2:796503/94‑1 (MQ=255)                                                                   gcgcGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTcc                                                               >  2:1098219/1‑139 (MQ=255)                                                                   gcgcGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGCCAGGTCGTTGCAACCGCAACCTTcc                                                               >  1:722391/1‑139 (MQ=255)                                                                      cgAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCtt                                                                 >  2:1566573/1‑134 (MQ=255)                                                                      cgAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCtt                                                                 <  1:1566573/134‑1 (MQ=255)                                                                        aGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTc                                                          >  2:188428/1‑139 (MQ=255)                                                                              cATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAAttt                                                    >  2:562165/1‑139 (MQ=255)                                                                                           aaCGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCa                                       <  1:1098219/139‑1 (MQ=255)                                                                                              gTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGAGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACgg                                    >  1:996319/1‑139 (MQ=255)                                                                                                      gCATGTTAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGt                            <  2:1653343/139‑1 (MQ=255)                                                                                                            tAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTACCACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                         >  1:193553/1‑136 (MQ=255)                                                                                                            tAACAACCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                         <  2:193553/136‑1 (MQ=255)                                                                                                              acaacCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                         <  1:616898/139‑1 (MQ=255)                                                                                                                 acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGgcg                                                                                                     <  2:496825/60‑1 (MQ=255)                                                                                                                 acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGgcg                                                                                                     >  1:496825/1‑60 (MQ=255)                                                                                                                 acCTGTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTcac                 >  1:770311/1‑139 (MQ=255)                                                                                                                     gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCt                                                         >  1:66826/1‑95 (MQ=255)                                                                                                                     gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCt                                                         <  2:66826/95‑1 (MQ=255)                                                                                                                     gTTCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGt             >  2:370966/1‑139 (MQ=255)                                                                                                                      ttCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGTCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCACCCTTCCGTCTCTAATTTCTCTTTTGGAGCACGGATTTCCGTGCCCATTTCACAAGtt            >  2:28231/1‑139 (MQ=255)                                                                                                                      ttCGATCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGtt            >  2:1176337/1‑139 (MQ=255)                                                                                                                           tCGTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTg       <  2:770311/139‑1 (MQ=255)                                                                                                                             gTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                         <  2:590539/119‑1 (MQ=255)                                                                                                                             gTGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGcc                         >  1:590539/1‑119 (MQ=255)                                                                                                                              tGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTa    >  2:755952/1‑139 (MQ=255)                                                                                                                              tGAATGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTa    >  1:942612/1‑139 (MQ=255)                                                                                                                                aaTGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATg  <  2:942612/139‑1 (MQ=255)                                                                                                                                aaTGCGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGCCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATg  <  2:722391/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                             GATCCGGCTAACTCCTTTAAAGTGGGAAGTTATACGGTCGTGGATGCGTTAGTACGTTATGATCTGGCGCGAGTCGGCATGGCTGGCTCCAACGTGGCGCTGCATGTTAACAACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCG‑‑‑‑‑CAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATG  >  NC_000913/169515‑169780

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