Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	812,748	C→T	A160V (GCA→GTA)	bioD →	dethiobiotin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	812,748	0	C	T	100.0%	113.3 / NA	31	A160V (GCA→GTA)	bioD	dethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/14);  total (17/14)

TAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAC  >  NC_000913/812673‑812824                                                                            |                                                                             tAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAgg                                                                         >  2:166218/1‑81 (MQ=255)   acacAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGt                                                                        >  2:166219/1‑80 (MQ=255)       aGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATac                                                                   <  2:166178/81‑1 (MQ=255)             aacTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACg                                                             >  1:166215/1‑81 (MQ=255)                tGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCg                                                          >  1:166217/1‑81 (MQ=255)                  ccGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCgg                                                         <  1:166179/80‑1 (MQ=255)                    ggTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGgact                                                      <  2:166180/81‑1 (MQ=255)                     gTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGgact                                                      >  1:166222/1‑80 (MQ=255)                        aTACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGgactgact                                                  >  2:166229/1‑81 (MQ=255)                        aTACTGGTAGTTGGTATGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGgactgact                                                  >  2:166212/1‑81 (MQ=255)                          aCTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGActct                                                >  1:166221/1‑81 (MQ=255)                             ggTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGc                                             >  2:166228/1‑81 (MQ=255)                                aGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCggg                                          >  1:166223/1‑81 (MQ=255)                                 gTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGt                                         <  2:166181/81‑1 (MQ=255)                                  ttGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGtt                                        <  2:166182/81‑1 (MQ=255)                                  ttGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGtt                                        >  2:166216/1‑81 (MQ=255)                                     gtgtGAAACTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTggg                                     <  2:166183/81‑1 (MQ=255)                                          aaaCTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGc                                 <  2:166184/80‑1 (MQ=255)                                            aCTCGGCTGTATTAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGaa                              >  2:166227/1‑81 (MQ=255)                                                       ttAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACg                    <  1:166185/80‑1 (MQ=255)                                                        tAATCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGcc                  >  1:166234/1‑81 (MQ=255)                                                          aTCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTc                <  2:166186/81‑1 (MQ=255)                                                          aTCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTc                >  1:166238/1‑81 (MQ=255)                                                           tCACGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTcc               <  1:166187/81‑1 (MQ=255)                                                             aCGCGATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCg              <  2:166188/80‑1 (MQ=255)                                                               gcgATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCggg            >  2:166235/1‑80 (MQ=255)                                                                 gATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGaa          >  2:166220/1‑80 (MQ=255)                                                                 gATGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGaa          >  1:166225/1‑80 (MQ=255)                                                                  aTGTTGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGaaaa        <  1:166189/81‑1 (MQ=255)                                                                     ttGACTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGt     <  2:166190/81‑1 (MQ=255)                                                                        aCTGTACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc  <  2:166191/81‑1 (MQ=255)                                                                            |                                                                             TAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAC  >  NC_000913/812673‑812824

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