Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,197,289	(C)5→4	coding (387/774 nt)	nudC →	NADH pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,197,285	0	C	.	100.0%	152.7 / NA	37	coding (383/774 nt)	nudC	NADH pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (19/18);  total (19/18)

GCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGCCCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGGTGTCCATACAGTACTTGCCGGATTCGTCGAAGTGGGCGAAACCCTCGAGCAGGCAGTCG  >  NC_000913/4197146‑4197416                                                                                                                                            |                                                                                                                                    gCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGc                                                                                                                                     <  1:602088/140‑1 (MQ=255)  ccGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGcc                                                                                                                                    >  2:488859/1‑140 (MQ=255)   cGGACGAGGCGTACAACTGGCGTAGTTTTACCGATCGCATAAATACTGCGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGccc                                                                                                                                   <  2:326271/140‑1 (MQ=255)        gAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGcccctgca                                                                                                                              <  1:164110/140‑5 (MQ=255)              tACAACTGGCGGAGTCTTACCGCTCACATCAATACTGAGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTAttg                                                                                                                       <  1:253967/140‑1 (MQ=255)                     ggCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTAttg                                                                                                                       >  2:650997/1‑133 (MQ=255)                     ggCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTAttg                                                                                                                       <  1:650997/133‑1 (MQ=255)                            ttttACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCATTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTcgc                                                                                                         <  1:584486/140‑1 (MQ=255)                             tttACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTcgcg                                                                                                        >  1:414757/1‑140 (MQ=255)                             tttACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTcgcg                                                                                                        <  1:470784/140‑1 (MQ=255)                              ttACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCga                                                                                                       <  1:488859/140‑1 (MQ=255)                              ttACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCga                                                                                                       <  1:698341/140‑1 (MQ=255)                                     tCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCg                                                                                                <  1:743284/140‑1 (MQ=255)                                            aaTACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGAtcctcc                                                                                         <  2:807343/140‑1 (MQ=255)                                                   tgGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCa                                                                                  <  2:604015/140‑1 (MQ=255)                                                   tgGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCa                                                                                  >  1:94836/1‑140 (MQ=255)                                                   tgGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCa                                                                                  >  1:89982/1‑140 (MQ=255)                                                   tgGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCa                                                                                  >  1:776182/1‑140 (MQ=255)                                                   tgGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCa                                                                                  >  2:494615/1‑140 (MQ=255)                                                     gTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGGGATGATTCGATCCTCCTCGCCCAGc                                                                                >  1:403137/1‑140 (MQ=255)                                                     gTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGc                                                                                >  1:411763/1‑140 (MQ=255)                                                        aCTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATa                                                                             >  2:604361/1‑140 (MQ=255)                                                        aCTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTACATTATTGTTGCCATCCGTCGCGACGAATCGATCGTCCTCGCCCAGCATa                                                                             >  1:328801/1‑140 (MQ=255)                                                          tGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATAcc                                                                           >  1:734882/1‑140 (MQ=255)                                                           gCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATAccc                                                                          >  1:813583/1‑140 (MQ=255)                                                             gggCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGc                                                                        <  1:114738/140‑1 (MQ=255)                                                               gCATGAAATGTCTCCAAAAAAACCCGAAGGGGCGACCCTGTACACCCTTTGCCGTCAGCGTTTCTCCCCGCAAA‑CAGCCCCTGACTTATTGTCGCCATCCGTAGCGATGATTCGATCCTCCTCGCCCACCATGCCCGCCa                                                                      <  2:527436/140‑1 (MQ=255)                                                                      aTGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAAc                                                               <  2:94836/140‑1 (MQ=255)                                                                        gTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACgg                                                             <  2:734882/140‑1 (MQ=255)                                                                         tATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTTTTGTTGCCATCCGTCGGGATGATTCGATCCTCCTCGCCCAGCATCCCCGCCATCGTAACGgt                                                            <  2:89982/140‑1 (MQ=255)                                                                         tATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGgt                                                            <  2:411763/140‑1 (MQ=255)                                                                           tCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGgtgt                                                          <  1:604361/140‑1 (MQ=255)                                                                           tCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCACGCCCAGCATACCCGCCATCGTAACGgtgt                                                          <  2:813583/140‑1 (MQ=255)                                                                           tCCGAGCAAAACCGAATCGGCGCTGCTGTGCCGCCATTGACGTGAGCGGTACTCCCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGgtgt                                                          <  2:254815/140‑1 (MQ=255)                                                                               aGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACg                                                              <  1:207298/132‑1 (MQ=255)                                                                               aGCAAAACCGAATGGGCGACGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACg                                                              >  2:207298/1‑132 (MQ=255)                                                                                gCAAAACCGAATGGGCACTGCTTTGCAGCCTCTGCGGTGTGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATTCCCGCCATCGTAACGGTGTCCATa                                                     <  2:403137/140‑1 (MQ=255)                                                                                    aaCCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCAAACCGGCTATCGTAACGGTGTCCATACAGt                                                 >  1:402719/1‑140 (MQ=255)                                                                                             ggCGATGCTGTGCAGCCATTGCCGTGAGCGTGACTACCCTCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGGTGTCCATACAGTACTTGCCgg                                        <  2:325082/140‑1 (MQ=255)                                                                                              gCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGGTGTCCATACAGTACTTGCCGGa                                       >  1:801861/1‑140 (MQ=255)                                                                                                            ccATTGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGGTGTCCAGACAATACTTGCCGAATTCGTCGAAGTggg                         >  1:294279/1‑140 (MQ=255)                                                                                                               ttGCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGGTGTCCATACAGTACTTGCCGGATTCGTCGAAGTGGGCGa                      <  2:801861/140‑1 (MQ=255)                                                                                                                 gCCGTGAGCGTTACTACCCGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGGTGTCCATACAGTACTTGCCGGATTCGTCGAAGTGGGCGaaa                    >  1:344410/1‑140 (MQ=255)                                                                                                                             actccCCGCAAATCG‑CCCCTAACTTATTGTTCCCAGCCGTCGCGATGAGTCGATCCTCCTCGCCCAGCATCCCCGCCATCGTAACGGTGTCCATACAGTACTTGCCGGATTCGTCGAAGTGGGCGAAACCCTCGAgcagg        <  2:328801/140‑1 (MQ=255)                                                                                                                                 cccGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGGTGTCCATACAGTACTTGCCGGATTCGTCGAAGTGGGCGAAACCCTCGAGCAGGCAGt    >  1:624396/1‑140 (MQ=255)                                                                                                                                   cGCAAATCG‑CCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGGTGTCCATACAGTACTTGCCGGATTCGTCGAAGTGGGCGAAACCCTCGAGCAGGCAGTcg  >  1:718106/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                    GCCGGACGAGGCGTACAACTGGCGGAGTTTTACCGATCGCATAAATACTGTGGTTACTGCGGGCATGAAATGTATCCGAGCAAAACCGAATGGGCGATGCTGTGCAGCCATTGCCGTGAGCGTTACTACCCGCAAATCGCCCCCTGCATTATTGTTGCCATCCGTCGCGATGATTCGATCCTCCTCGCCCAGCATACCCGCCATCGTAACGGTGTCCATACAGTACTTGCCGGATTCGTCGAAGTGGGCGAAACCCTCGAGCAGGCAGTCG  >  NC_000913/4197146‑4197416

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