Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,547,660	G→T	S87I (AGC→ATC)	fdnG →	formate dehydrogenase N subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,547,660	0	G	T	100.0%	148.3 / NA	43	S87I (AGC→ATC)	fdnG	formate dehydrogenase N subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (25/18);  total (25/18)

CGTAACACCTGCACATACTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTC  >  NC_000913/1547539‑1547792                                                                                                                          |                                                                                                                                     cGTAACACCTGCACATACTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCt                                                                                                                          <  2:554116/134‑1 (MQ=255) cGTAACACCTGCACATACTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCt                                                                                                                          >  1:554116/1‑134 (MQ=255)    aaCACCTGCACATACTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGaa                                                                                                                 <  2:344938/140‑1 (MQ=255)                 aCTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTg                                                                                                    >  1:43577/1‑140 (MQ=255)                   tGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCt                                                                                                  >  2:605204/1‑140 (MQ=255)                     ttCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTgg                                                                                                >  1:107454/1‑140 (MQ=255)                     ttCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTgg                                                                                                >  2:455884/1‑140 (MQ=255)                     ttCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTgg                                                                                                >  1:568386/1‑140 (MQ=255)                       ccGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGAt                                                                                              >  2:700010/1‑140 (MQ=255)                          tAGGTTGCAGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCGGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTAc                                                                                           >  2:573017/1‑140 (MQ=255)                               tGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGttt                                                                                                     <  1:13478/125‑1 (MQ=255)                               tGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGttt                                                                                                     >  2:13478/1‑125 (MQ=255)                               tGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTgg                                                                                                >  1:204408/1‑130 (MQ=255)                               tGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTgg                                                                                                <  2:204408/130‑1 (MQ=255)                               tGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTcga                                                                                      >  1:274133/1‑138 (MQ=255)                                gCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAAc                                                                                     <  1:605204/140‑1 (MQ=255)                                gCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAAc                                                                                     <  1:652035/140‑1 (MQ=255)                                       aTTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCGACAGTGaaa                                                                              <  2:274133/140‑1 (MQ=255)                                       aTTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGaaa                                                                              >  1:48123/1‑140 (MQ=255)                                                taGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTgc                                                                     >  1:556136/1‑140 (MQ=255)                                                 aGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTgcg                                                                    <  2:107454/140‑1 (MQ=255)                                                  gCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTgcgc                                                                   >  2:357016/1‑140 (MQ=255)                                                  gCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTgcgc                                                                   >  1:168048/1‑140 (MQ=255)                                                   ccTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCt                                                                  <  2:43577/140‑1 (MQ=255)                                                   ccTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCt                                                                  <  2:27512/140‑1 (MQ=255)                                                   ccTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCt                                                                  <  1:455884/140‑1 (MQ=255)                                                     tGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCt                                                                  >  2:662313/1‑138 (MQ=255)                                                     tGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCt                                                                  <  1:662313/138‑1 (MQ=255)                                                       ggTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTAccc                                                              >  1:454957/1‑140 (MQ=255)                                                        gTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCg                                                             >  1:186287/1‑140 (MQ=255)                                                                     aaaCGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTgcgc                                                <  2:168048/140‑1 (MQ=255)                                                                      aaCGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGcc                                               <  2:568386/140‑1 (MQ=255)                                                                                            acaTTGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAgcgc                         >  1:347519/1‑140 (MQ=255)                                                                                                tGAAGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTa                     >  1:251259/1‑140 (MQ=255)                                                                                                  aaGGTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGc                   >  2:50029/1‑140 (MQ=255)                                                                                                    ggTGACCCGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTg                 >  1:322789/1‑140 (MQ=255)                                                                                                         cccGGATCATCCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGgaag            <  2:454957/140‑1 (MQ=255)                                                                                                              atcatcCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCa       <  2:251259/140‑1 (MQ=255)                                                                                                              atcatcCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCa       >  1:726630/1‑140 (MQ=255)                                                                                                              atcatcCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCa       <  1:50029/140‑1 (MQ=255)                                                                                                                  tcCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATtct   >  2:97795/1‑140 (MQ=255)                                                                                                                  tcCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATtct   >  2:504788/1‑140 (MQ=255)                                                                                                                  tcCGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATtct   >  2:31698/1‑140 (MQ=255)                                                                                                                   ccGGTAATCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATtctc  <  2:347519/140‑1 (MQ=255)                                                                                                                          |                                                                                                                                     CGTAACACCTGCACATACTGTTCCGTAGGTTGCGGGCTATTGATGTATAGCCTGGGTGATGGCGCAAAAAACGCCAGAGAAGCGATTTATCACATTGAAGGTGACCCGGATCATCCGGTAAGCCGTGGTGCGCTGTGCCCGAAAGGGGCCGGTTTGCTGGATTACGTCAACAGTGAAAACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTC  >  NC_000913/1547539‑1547792

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