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%	45.7 / NA	17	A1074D (GCT→GAT)	yphG	DUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (9/8);  total (9/8)

GGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAA  >  NC_000913_3_pae_tpiA/2679382‑2679658                                                                                                                                                |                                                                                                                                      tggttttATTATCAATTGCTTTTTTGATGTATTTTGGGATGAGGGGGAGTTGGTGCCGATAAAATATTTCATTTAATAGGGTTTAGTGAACGCCGTTAAAAATAAAATGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGa                                                                                                                                  <  2:273596/146‑1 (MQ=255)  gTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGAt                                                                                                                                 >  2:317566/1‑149 (MQ=255)        attaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTCGTGAACCCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCg                                                                                                                           <  1:118783/149‑1 (MQ=255)            acaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGt                                                                                                                       >  1:161997/1‑149 (MQ=255)                                     tgatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGAc                                                                                                    <  1:531386/143‑1 (MQ=255)                                     tgatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGAc                                                                                                    >  2:531386/1‑143 (MQ=255)                                      gatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGATAATAGGGTTTAGTGAACGCGGTTAAAAAAAAAAGGCATCACGGAGGTGAATAATGCCGCTTTCAGCCAATCCGGATTAATCGCCGTCAATTCGGCGCAAGCAGACTCATAtt                                                                                             >  2:297717/1‑149 (MQ=255)                                           gCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTCCGCa                                                                                        >  2:371741/1‑149 (MQ=255)                                                   gcggtgcGCTCATGTATTAGAGTTAATAGGTTTTGGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcccc                                                                                <  2:130569/143‑1 (MQ=255)                                                                      cAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGc                                                             <  1:365989/149‑1 (MQ=255)                                                                         ttAATAGGGTTTAGTGAACACTGTTAAAAATAAAAAGCATCACGGTGGTTAATAATGCCGCTTTCGGCCAATCCAGATTAATCGAAATAAATTCAACACGAGCATACTCAACATGCGAAACTTCCCCCATTCCCGCGGACGCCAACATc                                                          <  2:367599/149‑1 (MQ=255)                                                                              aGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATAAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGccac                                                     >  2:59306/1‑149 (MQ=255)                                                                                               gTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCtt                                    >  2:170766/1‑149 (MQ=255)                                                                                                ttAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCttt                                   <  1:314629/149‑1 (MQ=255)                                                                                                                        gtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGAAGCCAGCATCGCAACAATCAGGCATTTTTCTTGATCCTGCTGTTGTAATTccgcgta           >  2:328945/1‑143 (MQ=255)                                                                                                                         tgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTCCTGTTGTAAATCGCCATAc          >  2:615978/1‑149 (MQ=255)                                                                                                                                 gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGacaa  <  1:317566/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                      GGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAA  >  NC_000913_3_pae_tpiA/2679382‑2679658

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