Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	2,679,525	G→T	A1074D (GCT→GAT)	yphG ←	DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	2,679,525	0	G	T	100.0%	79.8 / NA	26	A1074D (GCT→GAT)	yphG	DUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (10/16);  total (10/16)

TTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAA  >  NC_000913_3_pae_tpiA/2679384‑2679661                                                                                                                                              |                                                                                                                                          tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGAtt                                                                                                                                    >  1:34547/1‑149 (MQ=255)      attaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCg                                                                                                                               >  2:496191/1‑149 (MQ=255)       ttaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGc                                                                                                                              >  1:65761/1‑149 (MQ=255)            aacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCa                                                                                                                         >  1:131155/1‑149 (MQ=255)             acGGCTTATTTTAATTGTTTTGTGTGGAGCGGCAGCGCGTGCCGCTGAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGGAGCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATAGCCGTCaa                                                                                                                        <  2:113568/149‑1 (MQ=255)                 cTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCa                                                                                                                   <  1:29125/149‑1 (MQ=255)                                        aGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTgcgc                                                                                            >  1:1082/1‑149 (MQ=255)                                           ggCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAAc                                                                                         <  2:419461/149‑1 (MQ=255)                                                  cgTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATCTCAG‑GCTAGCAGACTCATATTGCGCAACCTCCCCg                                                                                  <  1:569290/149‑1 (MQ=255)                                                                                tttAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCACT‑TCCGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGCAGCCAGCATCGCCACCATc                                                    >  2:349816/1‑149 (MQ=255)                                                                                tgtAGTTAACGGTGGTAAAAATAAAAGGCATCATGGTGGTGAATAATGTCGCTTTCGGCCAATCCGGATTAATCGCCGTTAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATc                                                    <  2:282323/147‑1 (MQ=255)                                                                                 ttAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCTCGGACGCCAGCATCGCCACCATCa                                                   <  1:496191/149‑1 (MQ=255)                                                                                      gAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCGACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCAt                                              <  1:454074/149‑1 (MQ=255)                                                                                                aaaaaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAt                                    <  1:536000/149‑1 (MQ=255)                                                                                                  aaaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAtgc                                  >  1:424868/1‑149 (MQ=255)                                                                                                  aaaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAtgc                                  >  1:316055/1‑149 (MQ=255)                                                                                                                  ggtggtGAATAAGTCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCACCCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCACCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGc                  <  1:414755/149‑1 (MQ=255)                                                                                                                      gtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATa              >  1:207209/1‑149 (MQ=255)                                                                                                                         aataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAg           <  2:316055/149‑1 (MQ=255)                                                                                                                         aataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAg           <  2:424868/149‑1 (MQ=255)                                                                                                                            aaggcagCTTTCTGGCAATCCGGATTAATCGCCGCCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGACAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGTGa        <  2:183043/143‑1 (MQ=255)                                                                                                                             atGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGac       >  2:87217/1‑149 (MQ=255)                                                                                                                             atGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGac       <  1:196071/149‑1 (MQ=255)                                                                                                                               gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGacaa     <  2:574561/149‑1 (MQ=255)                                                                                                                                  gCTTTCGGCCAATCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGAcaacaa  <  1:461972/149‑1 (MQ=255)                                                                                                                                  gCTTTACGCCAATCCGGATTAATCGGCTTCAAT‑TCAGCGCGAGTAGACTCATATTGCGCAAACTCCCCGTGTTTCGATGAGTTGAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGAcaacaa  <  2:208140/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                          TTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAAT‑TCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAA  >  NC_000913_3_pae_tpiA/2679384‑2679661

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑