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

GGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCC  >  NC_000913/2679460‑2679579                                                                  |                                                       gggTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATcc                                                      >  2:778999/1‑68 (MQ=255) gggTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATcc                                                      >  2:1202460/1‑68 (MQ=255) gggATTAGTCAACGCGGTTAAAAATAAAAGGAATCACGGTGGAGATTAATTCCGCTTTCGGCCAATcc                                                      >  2:819859/1‑68 (MQ=255)   gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCgg                                                    >  1:1084891/1‑68 (MQ=255)    tttAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGa                                                   <  1:322481/68‑1 (MQ=255)        gTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTaa                                               >  2:274123/1‑68 (MQ=255)        gTGAACGCGGTTAAAAATAAAAGGCATCAAGGTGGTGAATAATGAAGCTTTCGGCCAATCCGGATTaa                                               >  1:374430/1‑68 (MQ=255)         tGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                              >  1:1391420/1‑68 (MQ=255)         tGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                              >  2:534298/1‑68 (MQ=255)             cgcgGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCACTCCGGATTAATCGcc                                          >  2:322239/1‑68 (MQ=255)              gcgGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCg                                         <  1:1267036/68‑1 (MQ=255)              gcgGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCg                                         <  1:756211/68‑1 (MQ=255)                ggTTAAACATAAAAGGCTTCTCGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTc                                       <  2:1105915/68‑1 (MQ=255)                        aTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCgg                                                    >  1:1251356/1‑47 (MQ=255)                        aTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCgg                                                    <  2:1251352/47‑1 (MQ=255)                        aTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAgc                               >  1:1079820/1‑68 (MQ=255)                        aTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAgc                               >  2:1063390/1‑68 (MQ=255)                        aTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAgc                               >  2:916194/1‑68 (MQ=255)                             aGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAg                          >  1:1043292/1‑68 (MQ=255)                             aGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAg                          >  1:1254414/1‑68 (MQ=255)                             aGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAg                          >  1:1434281/1‑68 (MQ=255)                                 aTCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                      <  2:265809/68‑1 (MQ=255)                                 aTCACGGCGGTGAATAAAGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                      <  1:862154/68‑1 (MQ=255)                                   cACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                    <  2:1084887/68‑1 (MQ=255)                                   cACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                    <  1:1456381/68‑1 (MQ=255)                                    aCGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGACGTCAATTCAGCGCGAGCAGACTc                   >  2:1219380/1‑68 (MQ=255)                                     cGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCa                  <  1:183787/68‑1 (MQ=255)                                         ggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATAtt              >  1:923947/1‑68 (MQ=255)                                         ggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATAtt              >  2:1399411/1‑68 (MQ=255)                                              ataatGCCGCTTTCGGCCAATCCGGATTACTCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCa         >  2:1124520/1‑68 (MQ=255)                                              ataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCa         >  2:696666/1‑68 (MQ=255)                                                aatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAc       >  2:396716/1‑68 (MQ=255)                                                 atGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAGcc      >  2:1426336/1‑68 (MQ=255)                                                  tGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCt     <  1:534301/68‑1 (MQ=255)                                                   gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAGCCTc    >  1:805261/1‑68 (MQ=255)                                                   gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc    >  1:995585/1‑68 (MQ=255)                                                   gCCGCTTTCGGCAAATCCGGATTAATCGCCGTCAATACAGCGAGAGCAGACTAATATTGCGCACCCTc    >  1:1066274/1‑68 (MQ=255)                                                     cGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTccc  <  2:1434276/68‑1 (MQ=255)                                                     cGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTccc  <  2:1043287/68‑1 (MQ=255)                                                                  |                                                       GGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCC  >  NC_000913/2679460‑2679579

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