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%	139.4 / NA	44	*226* (TAG→TAA)	bioD	dethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/23);  total (21/23)

ATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAGCCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAT  >  NC_000913/812834‑813072                                                                                                                   |                                                                                                                               atGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTgcg                                                                                                      >  1:1245786/1‑139 (MQ=255)   gACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTgcgag                                                                                                    <  2:477364/139‑1 (MQ=255)      cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTgcgagcga                                                                                                 >  2:872874/1‑139 (MQ=255)          cTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATc                                                                                             <  2:1245786/139‑1 (MQ=255)             aCCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATcgcg                                                                                          >  1:1006341/1‑139 (MQ=255)             aCCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATcgcg                                                                                          >  2:811830/1‑139 (MQ=255)                   aTGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCg                                                                                    >  2:1374658/1‑139 (MQ=255)                      aTTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAtt                                                                                 <  1:1374658/139‑1 (MQ=255)                          ccGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACt                                                                             <  1:886029/139‑1 (MQ=255)                           cgcgCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTg                                                                            <  1:1169747/139‑1 (MQ=255)                               ccgcTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGGTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              <  2:482627/134‑1 (MQ=255)                               ccgcTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGGTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              >  1:482627/1‑134 (MQ=255)                                         gagagaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAAtt                                                              <  1:890891/139‑1 (MQ=255)                                          agagaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTa                                                             <  1:714254/139‑1 (MQ=255)                                          agagaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTa                                                             <  1:811830/139‑1 (MQ=255)                                             gaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTa                                                             <  2:339782/136‑1 (MQ=255)                                             gaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTa                                                             >  1:339782/1‑136 (MQ=255)                                              aTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTgcgagc                                                                                                   <  2:1652729/97‑1 (MQ=255)                                              aTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTgcgagc                                                                                                   >  1:1652729/1‑97 (MQ=255)                                                  ccTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCATTTAAATCCTaa                                                     <  2:1006341/139‑1 (MQ=255)                                                       cTTGCAGAAAATCCAGGAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAaataa                                                >  2:356206/1‑139 (MQ=255)                                                          gCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTgcga                                                                                                     >  1:1732812/1‑83 (MQ=255)                                                          gCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTgcga                                                                                                     <  2:1732812/83‑1 (MQ=255)                                                          gCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTaaa                                                           >  1:1290069/1‑125 (MQ=255)                                                          gCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTaaa                                                           <  2:1290069/125‑1 (MQ=255)                                                           cAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCAc                                            >  2:374710/1‑139 (MQ=255)                                                                aaTCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACC‑TTCTGTTATTTGGTTAATTTGCGAGCAAGAT‑GCGGTTTTGAATTACTTCCATTTTAACAATTAAATTTTTAAAATATTACGaaaaa                                      >  2:622102/1‑139 (MQ=255)                                                                        aaaTGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTa                               >  2:11772/1‑139 (MQ=255)                                                                            gCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTgcgag                                                                                                    >  2:1516044/1‑66 (MQ=255)                                                                            gCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTgcgag                                                                                                    <  1:1516044/66‑1 (MQ=255)                                                                              ggCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              >  1:379030/1‑86 (MQ=255)                                                                              ggCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              <  2:379030/86‑1 (MQ=255)                                                                                 aaCCGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCt                      <  1:11772/139‑1 (MQ=255)                                                                                   ccGGAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCa                    >  1:900017/1‑139 (MQ=255)                                                                                      gAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCa                                             <  2:1381520/111‑1 (MQ=255)                                                                                      gAAAGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCa                                             >  1:1381520/1‑111 (MQ=255)                                                                                       aaaGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcg                <  2:877470/139‑1 (MQ=255)                                                                                        aaGTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcgg               <  1:872874/139‑1 (MQ=255)                                                                                          gTACATAAACCTTGCCTT‑GTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggcg             <  1:1005146/139‑1 (MQ=255)                                                                                                     ttgcctt‑gTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              >  1:1299360/1‑63 (MQ=255)                                                                                                     ttgcctt‑gTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              <  2:1299360/63‑1 (MQ=255)                                                                                                     ttgcctt‑gTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt  >  2:2312/1‑139 (MQ=255)                                                                                                       gcctt‑gTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAtt                                                                                 <  1:1523642/58‑1 (MQ=255)                                                                                                       gcctt‑gTTGTAACCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAtt                                                                                 >  2:1523642/1‑58 (MQ=255)                                                                                                                   |                                                                                                                               ATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTT‑GTTGTAGCCATTCTG‑TATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAT  >  NC_000913/812834‑813072

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