Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	169,743	A→G	intergenic (+16/‑35)	fhuA → / → fhuC	ferrichrome outer membrane transporter/phage receptor/iron(III) hydroxamate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	169,743	0	A	G	100.0%	91.2 / NA	29	intergenic (+16/‑35)	fhuA/fhuC	ferrichrome outer membrane transporter/phage receptor/iron(III) hydroxamate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/13);  total (16/13)

ACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCAC  >  NC_000913/169627‑169851                                                                                                                     |                                                                                                             acCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCa                                                                                              >  2:19679/1‑133 (MQ=255)     gTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTc                                                                                          >  1:257480/1‑133 (MQ=255)      ttCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTca                                                                                         >  2:40001/1‑133 (MQ=255)      ttCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTca                                                                                         >  1:501752/1‑133 (MQ=255)       tCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGATGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTcac                                                                                        >  1:80085/1‑133 (MQ=255)       tCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGATGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTcac                                                                                        >  1:99599/1‑133 (MQ=255)               gAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGc                                                                                <  1:454297/133‑1 (MQ=255)                 aTACGTCGTCAGCTGTTTAAACTCTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTg                                                                              <  1:249059/133‑1 (MQ=255)                       cGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTg                                                                              <  1:111200/127‑1 (MQ=255)                       cGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTg                                                                              >  2:111200/1‑127 (MQ=255)                       cGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGc                                                                        >  2:312738/1‑133 (MQ=255)                        gCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCa                                                                       >  2:102567/1‑133 (MQ=255)                                          tATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTc                                                     >  1:351548/1‑133 (MQ=255)                                           aTGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTcc                                                    >  1:361780/1‑133 (MQ=255)                                           aTGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGATGGATGTAATGCAGGAATACAAGAATCATTcc                                                    >  2:178572/1‑133 (MQ=255)                                                tgctTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATAc                                               <  2:50489/133‑1 (MQ=255)                                                tgctTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATAc                                               >  1:303550/1‑133 (MQ=255)                                                   ttCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCAc                                            >  2:313145/1‑133 (MQ=255)                                                     cTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGAAGGAATCGACGCATAATTCCGGTACCACtt                                          >  1:141801/1‑133 (MQ=255)                                                        gggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTg                                       <  2:501752/133‑1 (MQ=255)                                                         ggCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTTTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGc                                      <  2:229537/133‑1 (MQ=255)                                                          gcgcAGAACATCAGGCCGTTGCAACCGCCACCCTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCa                                     <  2:235124/133‑1 (MQ=255)                                                           cgcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCAc                                    <  1:390995/133‑1 (MQ=255)                                                           cgcAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCAc                                    <  1:313145/133‑1 (MQ=255)                                                             cAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTg                                  >  2:400084/1‑133 (MQ=255)                                                             cAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTg                                  <  2:430573/133‑1 (MQ=255)                                                                   gTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAat                            <  1:19679/133‑1 (MQ=255)                                                                     cAGGTCGTTGAAACCGCAACCTTCCGTTTCTATTTTCTCTTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAatat                          <  2:186664/133‑1 (MQ=255)                                                                                             cGTTTCTATTTTCTCGTTTGGGGCGCGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCAc  <  1:178572/133‑1 (MQ=255)                                                                                                                     |                                                                                                             ACCTGTTCGATCGTGAATACGTCGCCAGCTGCTTTAACACTTATGGCTGCTTCTGGGGCGCAGAACGTCAGGTCGTTGCAACCGCAACCTTCCGTTTCTAATTTCTCTTTTGGGGCACGGATTTCCGTGCCCATTTCACAAGTTGGCTGTTATGCAGGAATACACGAATCATTCCGATACCACTTTTGCACTGCGTAATATCTCCTTTCGTGTGCCCGGGCGCAC  >  NC_000913/169627‑169851

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