Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,100,923	A→G	S116P (TCA→CCA)	gnd ←	6‑phosphogluconate dehydrogenase, decarboxylating

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,100,923	0	A	G	100.0%	146.4 / NA	42	S116P (TCA→CCA)	gnd	6‑phosphogluconate dehydrogenase, decarboxylating
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/22);  total (20/22)

CGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGAC  >  NC_000913/2100785‑2101066                                                                                                                                           |                                                                                                                                                cggcggCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCAc                                                                                                                                         <  3:105228/147‑1 (MQ=255)   gcggcgATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGAt                                                                                                                                      >  1:116281/1‑148 (MQ=255)                   ggATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTcc                                                                                                                     <  1:117071/149‑1 (MQ=255)                    gATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCt                                                                                                                    <  5:68996/149‑1 (MQ=255)                         gTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTgg                                                                                                                  >  6:164423/1‑146 (MQ=255)                         gTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTgg                                                                                                                  >  4:41942/1‑146 (MQ=255)                                   ttCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTac                                                                                                     <  2:116281/149‑1 (MQ=255)                                    tCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTAAAAGCCCTCTGCTGGAAGCTCACGATTacga                                                                                                                                 >  4:27439/1‑120 (MQ=255)                                       tAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTaccacc                                                                                                 <  1:120589/149‑1 (MQ=255)                                          gCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTacc                                                                                                    >  1:186030/1‑143 (MQ=255)                                          gCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTacc                                                                                                    <  2:186030/143‑1 (MQ=255)                                                       caccaGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAgg                                                                                                           <  1:183620/123‑1 (MQ=255)                                                       caccaGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAgg                                                                                                           >  2:183620/1‑123 (MQ=255)                                                           aGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCtt                                                                             >  2:197189/1‑149 (MQ=255)                                                            ggCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTc                                                                                                                      <  8:104281/107‑1 (MQ=255)                                                            ggCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTc                                                                                                                      >  7:104281/1‑107 (MQ=255)                                                             gCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGgtg                                                                                                         >  1:18447/1‑119 (MQ=255)                                                             gCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGgtg                                                                                                         <  2:18447/119‑1 (MQ=255)                                                                taataGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCg                                                                        <  1:197189/149‑1 (MQ=255)                                                                   taGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTTCCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCgaga                                                                     <  3:41942/149‑1 (MQ=255)                                                                   taGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTTTTACCACCATCAATGATGATGTCTCCTTTATCgaga                                                                     <  5:164423/149‑1 (MQ=255)                                                                               tCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGatat                                                                  <  1:159608/140‑1 (MQ=255)                                                                               tCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGatat                                                                  >  2:159608/1‑140 (MQ=255)                                                                                cAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAgg                                                         <  5:138015/148‑1 (MQ=255)                                                                                  gcgcCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAatgatg                                                                                       <  1:145731/116‑1 (MQ=255)                                                                                  gcgcCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAatgatg                                                                                       >  2:145731/1‑116 (MQ=255)                                                                                   cgcCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTcc                                                                               <  8:1590/123‑1 (MQ=255)                                                                                   cgcCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTcc                                                                               >  7:1590/1‑123 (MQ=255)                                                                                            ttCACCGCCAGAAACACCGGTACCGATGAACTTAACGCCCTCTGCTGGAAGCTCACCATTACGACGAATAGTGTACTGGAAGAAGGTGTTACCACCATCAATGATGAAGTCTCCTATATCAAGATATGGATTGAGAGAAACAATagaag                                            >  8:154112/1‑149 (MQ=255)                                                                                                  gCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAatgatg                                                                                       >  3:47807/1‑100 (MQ=255)                                                                                                  gCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAatgatg                                                                                       <  4:47807/100‑1 (MQ=255)                                                                                                   ccAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACAACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCATCCGt                                     >  1:170671/1‑149 (MQ=255)                                                                                                         agacCTGTACCGATGTAGTTAACGCGGTCTGCTGGAAGCTCACGATTACGACGAATAGTGTTCTGGAAGAAGGTGGTACCACCATCAATGATGTTGGCTCCTTTTTCGATATATGGTTTGAGGGAATCAATAGCAGCTTGCGTGCCTGc                               <  7:154112/147‑1 (MQ=255)                                                                                                                  ccGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCATCCGTGCCTGCACCTGCttt                      >  7:51/1‑149 (MQ=255)                                                                                                                   cGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAAt                                                     >  2:180926/1‑117 (MQ=255)                                                                                                                   cGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAAt                                                     <  1:180926/117‑1 (MQ=255)                                                                                                                   cGATGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCATCCGTGCCTGCACCTGCttt                      >  1:38335/1‑148 (MQ=255)                                                                                                                     aTGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAgg                                                         >  3:155006/1‑111 (MQ=255)                                                                                                                     aTGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAgg                                                         <  4:155006/111‑1 (MQ=255)                                                                                                                     aTGAAGTTAAAGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCATCCGTGCCTGCACCTGCTTTCAc                   <  8:51/149‑1 (MQ=255)                                                                                                                              aaGCCCTCTGCTGGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCATTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAg          <  2:170671/149‑1 (MQ=255)                                                                                                                                        ctgGAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGcgac  >  5:32811/1‑147 (MQ=255)                                                                                                                                           |                                                                                                                                                CGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGAC  >  NC_000913/2100785‑2101066

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