Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	812,947	G→A	*226* (TAG→TAA)	bioD →	dethiobiotin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	812,947	0	G	A	100.0%	117.3 / NA	37	*226* (TAG→TAA)	bioD	dethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/17);  total (20/17)

GAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAATTGCCCAGG  >  NC_000913/812814‑813080                                                                                                                                        |                                                                                                                                        gAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCAtt                                                                                                                                    <  2:776038/139‑1 (MQ=255)        tCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTAtt                                                                                                                             <  1:868272/139‑1 (MQ=255)                     atGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcg                                                                                                                >  2:275988/1‑139 (MQ=255)                          cACGCTCACCCGCATGATTGCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgagc‑‑ga                                                                                                         >  1:702609/1‑139 (MQ=255)                          cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgagc‑‑ga                                                                                                         >  1:802719/1‑139 (MQ=255)                            cGCTCACCCGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCAtt                                                                                                                                    <  2:624106/112‑1 (MQ=255)                            cGCTCACCCGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCAtt                                                                                                                                    >  1:624106/1‑112 (MQ=255)                              cTCACCCGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATc                                                                                                     <  2:186783/139‑1 (MQ=255)                              cTCACCCGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATc                                                                                                     >  2:535659/1‑139 (MQ=255)                                ctcccGCATGCTTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTGTCAGAAAATCCAGAAAATGCGTCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATcgc                                                                                                   <  2:979249/137‑1 (MQ=255)                                 aCCCGCATGATTCACGCACCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATcgcg                                                                                                  >  1:660075/1‑139 (MQ=255)                                  cccGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGt                                                                                                 <  1:597539/139‑1 (MQ=255)                                  cccGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGt                                                                                                 >  1:483719/1‑139 (MQ=255)                                       aTGATTCCCGCGCCGCTGCTGGGAGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCtt                                                                                              >  2:960867/1‑139 (MQ=255)                                             cccGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGAc                                                                                      <  1:435263/139‑1 (MQ=255)                                              ccGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACt                                                                                     <  2:483719/139‑1 (MQ=255)                                                 cgccgcTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCa                                                                                  >  2:699185/1‑139 (MQ=255)                                                               gagagaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAAtt                                                                      <  1:960867/139‑1 (MQ=255)                                                                     tCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTc                                                                >  2:167502/1‑139 (MQ=255)                                                                       cccTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTa                                                              <  1:422781/139‑1 (MQ=255)                                                                          tGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTaaaa                                                           >  1:950140/1‑139 (MQ=255)                                                                                gCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCa                                                     >  1:284312/1‑139 (MQ=255)                                                                                 cAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGAc                                                                                      <  1:756230/105‑1 (MQ=255)                                                                                 cAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAAGTGCGAGC‑‑GAGATCGCGTCTTCGATTGAc                                                                                      >  2:756230/1‑105 (MQ=255)                                                                                      aaTCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGaaaa                                               >  1:46793/1‑139 (MQ=255)                                                                                                               agaaCATAAACCTTGCCTTGTTGTAACCTTTCTGTATTTGTTTAAATTGCGACC‑‑GAGTTCGCGTCTTCGATTGACTGCAATTAAACAAATAAAATTCTAAAATAATCCGGAAAAAAATTTTACTTCCGCCTCATgcggc                      <  1:699185/136‑1 (MQ=255)                                                                                                               aGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggc                      <  2:284312/139‑1 (MQ=255)                                                                                                                gTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggct                     <  2:660075/139‑2 (MQ=255)                                                                                                                           ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt          >  1:283412/1‑139 (MQ=255)                                                                                                                           ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt          >  2:855064/1‑139 (MQ=255)                                                                                                                           ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt          >  1:52555/1‑139 (MQ=255)                                                                                                                           ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt          >  2:355981/1‑139 (MQ=255)                                                                                                                                  gttgtAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCt                                                                                               <  1:495718/47‑1 (MQ=255)                                                                                                                                  gttgtAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCt                                                                                               >  2:495718/1‑47 (MQ=255)                                                                                                                                   ttgtAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAATTGCCCAgg  <  2:912174/139‑1 (MQ=255)                                                                                                                                   ttgtAACCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAATTGCCCAgg  <  1:1040446/139‑1 (MQ=255)                                                                                                                                    tgtAACCATTCTGTATTTGGTTAAATTGCGAGCGAGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAATTGCCCAg   >  2:379609/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        GAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGG‑‑GAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGC‑‑GAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAATTGCCCAGG  >  NC_000913/812814‑813080

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