Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	413,307	T→C	D816G (GAC→GGC)	sbcC ←	exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	413,307	0	T	C	100.0%	167.8 / NA	47	D816G (GAC→GGC)	sbcC	exonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/14);  total (33/14)

TTGTTGTTGCTGACGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGTCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCC  >  NC_000913/413162‑413453                                                                                                                                                  |                                                                                                                                                   ttgttgttgCTGACGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGt                                                                                                                                               >  2:641940/1‑151 (MQ=255)     tgttgCTGACGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCgg                                                                                                                                            >  1:419904/1‑150 (MQ=255)        tgCTGACGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGtgtt                                                                                                                                        >  1:433950/1‑151 (MQ=255)         gCTGACGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGtgttg                                                                                                                                       >  2:628596/1‑151 (MQ=255)         gCTGACGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGtgttg                                                                                                                                       >  2:479546/1‑151 (MQ=255)          cTGACGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGtgttgt                                                                                                                                      >  2:493753/1‑151 (MQ=255)            gACGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGttgttg                                                                                                                                    >  2:799209/1‑151 (MQ=255)            gACGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGttgttg                                                                                                                                    >  2:135280/1‑151 (MQ=255)             aCGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGttgttg                                                                                                                                    >  1:627218/1‑150 (MQ=255)              cGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGAt                                                                                                                                  >  1:279478/1‑151 (MQ=255)                   aTCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCtgt                                                                                                                             >  2:97938/1‑151 (MQ=255)                      tGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGcc                                                                                                                          >  1:744174/1‑151 (MQ=255)                       gCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCGCCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGc                                                                                                                           <  1:369198/149‑1 (MQ=255)                         aTCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGc                                                                                                                       <  2:268140/151‑1 (MQ=255)                          tCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGc                                                                                                                       <  2:559805/150‑1 (MQ=255)                                  cAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTctgct                                                                                                              <  2:343408/151‑1 (MQ=255)                                  cAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTctgct                                                                                                              <  2:301449/151‑1 (MQ=255)                                             gAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGAc                                                                                                   >  1:456786/1‑151 (MQ=255)                                             gAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGAc                                                                                                   >  2:380032/1‑151 (MQ=255)                                                    gCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGtt                                                                                            <  1:706247/151‑1 (MQ=255)                                                    gCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGtt                                                                                            <  2:136901/151‑1 (MQ=255)                                                        tGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTgc                                                                                         <  2:816706/150‑1 (MQ=255)                                                        tGACTCGTGGGGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTgcg                                                                                        >  1:217891/1‑151 (MQ=255)                                                         gACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTgcgc                                                                                       >  2:208887/1‑151 (MQ=255)                                                           cTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCt                                                                                     >  1:79002/1‑151 (MQ=255)                                                                ggtgTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACgg                                                                                >  1:643151/1‑151 (MQ=255)                                                                   gtTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGgcgc                                                                             >  2:74038/1‑151 (MQ=255)                                                                           gCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTc                                                                     <  1:599826/151‑1 (MQ=255)                                                                                       gAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCtt                                                         >  1:635690/1‑151 (MQ=255)                                                                                            tGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCt                                                    >  1:497236/1‑151 (MQ=255)                                                                                               gcTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGtt                                                 <  1:813836/151‑1 (MQ=255)                                                                                                cTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTc                                                >  2:596635/1‑151 (MQ=255)                                                                                                 tAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTcc                                               >  1:451484/1‑151 (MQ=255)                                                                                                        tGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCt                                           <  2:195969/148‑1 (MQ=255)                                                                                                           tGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGtt                                      >  1:642744/1‑150 (MQ=255)                                                                                                              aTCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAgtg                                  >  2:44126/1‑151 (MQ=255)                                                                                                                  gCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGtgtttg                              >  1:645488/1‑151 (MQ=255)                                                                                                                  gCTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGtgtttg                              <  2:39605/151‑1 (MQ=255)                                                                                                                   cTCCACCGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGtgtttgt                             >  1:803051/1‑151 (MQ=255)                                                                                                                     ccACCGTCACAGTGAGAGCCAACCCGTCACCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTc                           >  1:790912/1‑151 (MQ=255)                                                                                                                        ccGTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGATTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATc                        >  2:671831/1‑151 (MQ=255)                                                                                                                          gTCACAGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATcc                       >  1:605979/1‑150 (MQ=255)                                                                                                                              caGTGAGAGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCAtt                    >  2:173928/1‑149 (MQ=255)                                                                                                                                  gagaGCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAgcgc               <  2:113310/150‑1 (MQ=255)                                                                                                                                   agagCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGcgccg             <  1:326198/151‑1 (MQ=255)                                                                                                                                   agagCCAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGcgccg             >  2:819188/1‑151 (MQ=255)                                                                                                                                        cAACCCGTCGCCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAgg        >  1:249297/1‑151 (MQ=255)                                                                                                                                              gtcgccAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGcc  <  2:707028/151‑1 (MQ=255)                                                                                                                                              gtcgccAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGcc  >  2:312353/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TTGTTGTTGCTGACGGTTATCTGCATCCTGCTTCAGCTGCTGGCGAATCTCGCCTTGACTCGTGGTGTTTTCACGCAACTTTTGGTGAGTTTGCGCTAACTCTTGCTGAATCTGCTCCACCGTCACAGTGAGAGCCAACCCGTCGTCAGGTCGGTGTTGTTGATGCTGTGCCAGCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCC  >  NC_000913/413162‑413453

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