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

GGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCC  >  NC_000913/2679460‑2679587                                                                  |                                                               gggTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATcc                                                              >  2:92898/1‑68 (MQ=255) gggTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATcc                                                              >  1:923881/1‑68 (MQ=255) gggTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATcc                                                              >  1:801811/1‑68 (MQ=255) gggTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATcc                                                              >  1:758104/1‑68 (MQ=255)   gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCgg                                                            >  2:771297/1‑68 (MQ=255)   gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCgg                                                            >  2:121849/1‑68 (MQ=255)   gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCgg                                                            >  2:8971/1‑68 (MQ=255)        gTGAGCGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTaa                                                       >  2:131831/1‑68 (MQ=255)        gTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTaa                                                       >  2:73261/1‑68 (MQ=255)         tGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                      >  2:395946/1‑68 (MQ=255)         tGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                      >  1:814318/1‑68 (MQ=255)              gcgGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCg                                                 <  1:791429/68‑1 (MQ=255)              gcgGTTAAAAATAAAAGGCATCACGGCGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCg                                                 <  2:460141/68‑1 (MQ=255)                      aaaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCAGTCATTTCa                                         >  1:943312/1‑68 (MQ=255)                        aTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                      <  2:511372/53‑1 (MQ=255)                        aTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAgc                                       >  1:899194/1‑68 (MQ=255)                        aTAAAAGGCATCACGGTGGAGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                      >  1:511372/1‑53 (MQ=255)                            aaGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGa                                   <  1:48326/68‑1 (MQ=255)                            aaGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGa                                   >  2:733621/1‑68 (MQ=255)                               gCATCAAGGTGGTGAATAATGCCGCTTTCGGCAAATCCGGATTAATCGCCGTAAATTAAGAGCGAgaa                                >  1:292583/1‑66 (MQ=255)                                 aTCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                              <  1:8971/68‑1 (MQ=255)                                 aTCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                              <  1:549878/68‑1 (MQ=255)                                     cGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCa                          >  1:595080/1‑68 (MQ=255)                                     cGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCa                          >  1:588185/1‑68 (MQ=255)                                      ggtggtGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCat                         >  2:20869/1‑68 (MQ=255)                                      ggtcgtTAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCATTTCAGCGCGAGCAGACTCat                         <  1:184290/68‑1 (MQ=255)                                         ggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATAtt                      >  2:61724/1‑68 (MQ=255)                                          gtgAATAATGCCTCTTTCGGCCAATCCGGATTTATCTCCGTCATTTCAGCTCGAGCAGACTCATATTg                     <  2:292583/68‑1 (MQ=255)                                             aataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCGATTCAGCGCGAGCAGACTCATATTgcgc                  >  2:807425/1‑68 (MQ=255)                                              ataatGCCGCTTTCGGCCTATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGCCTCATATTGCGCa                 >  2:892254/1‑68 (MQ=255)                                              ataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgc                    <  1:962269/65‑1 (MQ=255)                                              ataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgc                    >  2:962266/1‑65 (MQ=255)                                              ataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCa                 >  2:465245/1‑68 (MQ=255)                                              ataatGCAGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCa                 >  1:109152/1‑68 (MQ=255)                                                   gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc            >  2:299867/1‑68 (MQ=255)                                                   gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc            >  2:15799/1‑68 (MQ=255)                                                   gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc            >  2:963660/1‑68 (MQ=255)                                                    ccGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcc           <  2:899191/68‑1 (MQ=255)                                                     cGCTTTCGGCCAATCCTTAATAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTccc          <  1:892257/68‑1 (MQ=255)                                                        tttCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCgc       >  2:779369/1‑67 (MQ=255)                                                             gCCAATCCGGATTAATCGCCGTCAATTCAGAGCGAGCAGACTCATATTGAGCAACCTCCACGAGTcca  >  1:611115/1‑67 (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 < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑