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

GGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAA  >  NC_000913_3_pae_tpiA/2679381‑2679661                                                                                                                                                 |                                                                                                                                         tgggttttATTATCAACGGCTTATTCTAATTATTTTGTGATGGGGGGCAGGGTGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCgg                                                                                                                                      <  2:518199/145‑1 (MQ=255) gggTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCgg                                                                                                                                      <  1:398/149‑1 (MQ=255)   gTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGAt                                                                                                                                    >  1:455347/1‑149 (MQ=255)       ttattaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                                                                                                <  1:18812/149‑1 (MQ=255)        tattaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATc                                                                                                                               <  1:218815/149‑1 (MQ=255)            aacaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCg                                                                                                                           >  1:502581/1‑149 (MQ=255)                    cTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAAGAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCa                                                                                                                   <  1:72493/149‑1 (MQ=255)                          cTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCCATTGGTGTCAATCCGGATTAATCGCCGTCAATTCAGCGCGa                                                                                                             <  1:316522/149‑1 (MQ=255)                               tattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                                                                                                        <  1:106591/149‑1 (MQ=255)                                       gatgaGCGGCAGCGCGTGCCACTCAAATATTACAGTTAACAGCGCGTCGCGACCGCGATTAACAATCCAAGGCATTTCGACGTTATATAATGCCGCTTTCGGCCAATCCGGATTAATCCCCCTCAATTCACCGCGCGCAGACTCATCtt                                                                                                >  2:162462/1‑149 (MQ=255)                                                    gcgTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcccc                                                                                   <  1:186016/149‑1 (MQ=255)                                                         ccGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTc                                                                              >  2:359106/1‑149 (MQ=255)                                                           gCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCATGTCATCGCCAGCAGACTCATATTGCGCAACCTCCCCGAGTccc                                                                            >  1:177554/1‑149 (MQ=255)                                                             tCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCCAGTccccc                                                                          >  1:585105/1‑147 (MQ=255)                                                                    ttACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGcc                                                                   >  1:114585/1‑149 (MQ=255)                                                                     tACAGTTCATATTGTTTAGTGATAGAGGTTAAAAATCAAAGGCGATACGGTAGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCa                                                                  <  2:297477/149‑1 (MQ=255)                                                                             atatGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcccc                                                                                   <  2:613626/124‑1 (MQ=255)                                                                             aTAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcccc                                                                                   >  1:613626/1‑124 (MQ=255)                                                                              tAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATGAGGGTGGTGAAGAATGCCGCTTTCGGCCAATCCGGCTTAATCGCCGTCAATTCAGGGCGCGCAGCCTGATATTGCGCAACCTCGCGAAGTCCCGCGGGCGCCTGCATCGcca                                                         >  2:339615/1‑149 (MQ=255)                                                                                    ttAGTGACCGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCa                                                   <  1:74100/149‑1 (MQ=255)                                                                                        tGAACGCGGTTAAAAATAGAAGGAATCACGGTGGTGAATAATGCCGCTTTCGGCCGATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCAGATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGct                                               >  1:593540/1‑148 (MQ=255)                                                                                          aaCGCGGTTAAAAATAAAAGTCATCACGGTGGTGAATAATGCTGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCATCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCAtt                                             <  2:472190/149‑1 (MQ=255)                                                                                                 ttAAAAATAAAAGGCATCACGGTGGTGAATAATGCAGCTTTCGGCCAGTCCGGATTAATCGACGCCGATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCttt                                      >  1:604598/1‑149 (MQ=255)                                                                                                    aaaaTACAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCACTCCGGATTAATCGCCGTCACGTCAGCGCAAGCCGGCTCGGAGGTGGTGACCTTGTCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAtg                                   >  1:128228/1‑149 (MQ=255)                                                                                                    aaaaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAtg                                   <  2:114585/149‑1 (MQ=255)                                                                                                        tAAAAGGCATCACGGTGGTGAATCTTGAAGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGCGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCCCCATCAGGCATTTTTCTTTAtcctgc                               >  1:479082/1‑149 (MQ=255)                                                                                                                aTCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCAGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCGTTATGCTGCTGTTGTaa                       >  1:337040/1‑149 (MQ=255)                                                                                                                                 tGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATa              <  1:367430/141‑1 (MQ=255)                                                                                                                                 tGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATa              >  2:367430/1‑141 (MQ=255)                                                                                                                                     gCTTTCGGCCAATCCTGATTAATCGCCGCCAATTCAGCGCGAGCAGACTCATCTTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGAcaacaa  <  2:632623/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                         GGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAA  >  NC_000913_3_pae_tpiA/2679381‑2679661

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: ‑