Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,197,168	C→A	A89E (GCG→GAG)	nudC →	NADH pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,197,168	0	C	A	100.0%	144.7 / NA	43	A89E (GCG→GAG)	nudC	NADH pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/20);  total (23/20)

GGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGCCCCCTGCATTATTGT  >  NC_000913/4197041‑4197300                                                                                                                                |                                                                                                                                     ggCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCg                                                                                                                          <  1:288619/140‑1 (MQ=255)   cGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGAt                                                                                                                        <  1:21140/140‑1 (MQ=255)   cGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGAt                                                                                                                        >  2:752566/1‑140 (MQ=255)   cGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGAt                                                                                                                        >  1:319534/1‑140 (MQ=255)   cGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGAt                                                                                                                        <  2:627175/140‑1 (MQ=255)   cGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGAt                                                                                                                        <  1:39317/140‑1 (MQ=255)     aaTGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCg                                                                                                                      <  1:752566/140‑1 (MQ=255)     aaTGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTAAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGAt                                                                                                                        <  1:12581/138‑1 (MQ=255)     aaTGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTAAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGAt                                                                                                                        >  2:12581/1‑138 (MQ=255)        ggCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCAt                                                                                                                   <  2:275228/140‑1 (MQ=255)             ggggAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATa                                                                                                              <  1:508582/140‑1 (MQ=255)                          ggTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGc                                                                                                 <  2:319534/140‑1 (MQ=255)                           gTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCg                                                                                                <  2:71376/140‑1 (MQ=255)                            ttAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCgg                                                                                               >  1:486803/1‑140 (MQ=255)                                     cagcagCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAAt                                                                                      >  1:535028/1‑140 (MQ=255)                                         agcGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTAt                                                                                  >  1:594959/1‑140 (MQ=255)                                                             tCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGc                                                              >  2:509135/1‑140 (MQ=255)                                                              cGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATggg                                                               >  2:634680/1‑138 (MQ=255)                                                              cGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATggg                                                               <  1:634680/138‑1 (MQ=255)                                                               ggTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGa                                                            >  1:532366/1‑140 (MQ=255)                                                                gTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGGGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGAt                                                           >  2:535058/1‑140 (MQ=255)                                                                gTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGAt                                                           >  1:413275/1‑140 (MQ=255)                                                                gTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGAt                                                           >  2:98643/1‑140 (MQ=255)                                                                gTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGAATCCGAGCAAAACCGAATGGGCGAt                                                           >  1:149404/1‑140 (MQ=255)                                                                    gTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCtg                                                       >  2:761213/1‑140 (MQ=255)                                                                       aGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGc                                                    <  1:509135/140‑1 (MQ=255)                                                                             ttGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCAt                                              <  1:761213/140‑1 (MQ=255)                                                                                      aTGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAg                                     >  2:537962/1‑140 (MQ=255)                                                                                          tGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGtt                                 >  1:574554/1‑140 (MQ=255)                                                                                           gggCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTta                                >  1:473211/1‑140 (MQ=255)                                                                                           gggCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTta                                >  1:394337/1‑140 (MQ=255)                                                                                               tGTTTCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTtactac                            <  1:98643/140‑1 (MQ=255)                                                                                                 tttCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTAccc                          >  1:361025/1‑140 (MQ=255)                                                                                                   tCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGc                        <  2:473211/140‑1 (MQ=255)                                                                                                   tCAACTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGc                        <  1:52161/140‑1 (MQ=255)                                                                                                      aCTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCaaa                     >  2:260579/1‑140 (MQ=255)                                                                                                      aCTGGCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCaaa                     >  2:697734/1‑140 (MQ=255)                                                                                                          gCCGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGc                 <  2:149404/140‑1 (MQ=255)                                                                                                           ccGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGcc                >  1:229316/1‑140 (MQ=255)                                                                                                           ccGGACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGcc                >  2:180937/1‑140 (MQ=255)                                                                                                             ggACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGcccc              <  1:535058/140‑1 (MQ=255)                                                                                                             ggACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGcccc              <  1:260579/140‑1 (MQ=255)                                                                                                             ggACGAGGCGTACAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGcccc              <  1:654286/140‑1 (MQ=255)                                                                                                                         cAACTGGAGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGCCCCCTGCATTAttgt  >  1:81275/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                     GGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGCGTCACGATATGGGGTCGGTACGTCAGGTCATTGATCTCGATGTTGGGCTGTTTCAACTGGCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGCCCCCTGCATTATTGT  >  NC_000913/4197041‑4197300

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑