Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	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	2,679,525	0	G	T	100.0%	105.0 / NA	33	A1074D (GCT→GAT)	yphG	DUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (19/14);  total (19/14)

GGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCC  >  NC_000913/2679460‑2679587                                                                  |                                                               gggTTTAGTGAACGCGGTTAAAAATAAAAGGCATAACGGTGGTGAATAATGACGCTTTCGGCCAATcc                                                              >  1:295399/1‑68 (MQ=255)   gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCACTCCgg                                                            >  1:758114/1‑68 (MQ=255)        gTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTaa                                                       >  2:1105125/1‑68 (MQ=255)        gTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTaa                                                       >  2:382022/1‑68 (MQ=255)         tGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                      >  1:613162/1‑68 (MQ=255)         tGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                      >  2:692412/1‑68 (MQ=255)         tGAACGCGGATAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                      >  2:669860/1‑68 (MQ=255)          gAACGCGGTTAAAAATTAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATc                                                     >  2:939349/1‑68 (MQ=255)            aCGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGc                                                   <  2:250234/68‑1 (MQ=255)                  ttAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCaa                                             <  2:602747/68‑1 (MQ=255)                       aaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAg                                        <  1:930828/68‑1 (MQ=255)                        aTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAgc                                       >  2:699542/1‑68 (MQ=255)                        aTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAgc                                       >  1:670684/1‑68 (MQ=255)                             agtcATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCTCCGTCAATTCAGCGCGAg                                  <  1:303731/65‑1 (MQ=255)                             aGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCAGTCAATTCAGCGCGAg                                  >  1:1026526/1‑68 (MQ=255)                                   cACGGTGTTTAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                            <  2:972090/68‑1 (MQ=255)                                   cACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                            <  1:156666/68‑1 (MQ=255)                                   cACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                            <  2:945887/68‑1 (MQ=255)                                   cACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                            <  2:173437/68‑1 (MQ=255)                                   cACGGTCGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                            <  1:171856/68‑1 (MQ=255)                                     cGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCa                          <  1:654474/68‑1 (MQ=255)                                      ggtggtGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCat                         >  1:984146/1‑68 (MQ=255)                                      ggtggtGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCat                         >  1:568040/1‑68 (MQ=255)                                              ataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCa                 >  1:288471/1‑68 (MQ=255)                                                aatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCagc               >  2:894640/1‑66 (MQ=255)                                                   gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                              >  1:955799/1‑50 (MQ=255)                                                   gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                              <  2:955796/50‑1 (MQ=255)                                                   gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc            >  1:797613/1‑68 (MQ=255)                                                   gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc            <  1:338761/68‑1 (MQ=255)                                                    ccGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcc           <  1:721913/68‑1 (MQ=255)                                                     cGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTAAGCGCGAGCAGACTCATATTGCGAAACCTccc          >  2:397366/1‑68 (MQ=255)                                                          tCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGt     <  2:444080/68‑1 (MQ=255)                                                             gCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTccc  >  1:369014/1‑68 (MQ=255)                                                                  |                                                               GGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCC  >  NC_000913/2679460‑2679587

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