Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	813,338	A→G	intergenic (+391/‑188)	bioD → / → uvrB	dethiobiotin synthetase/excision nuclease subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	813,338	0	A	G	100.0%	171.5 / NA	50	intergenic (+391/‑188)	bioD/uvrB	dethiobiotin synthetase/excision nuclease subunit B
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/27);  total (23/27)

AAAGAAATATCTTTTATTTTTTAACTGGTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAACCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGT  >  NC_000913/813210‑813459                                                                                                                                 |                                                                                                                          aaaGAAATATCTTTTATTTTTTAACTGGTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGc                                                                                                                >  2:293752/1‑140 (MQ=255)   aGAAATATCTTTTATTTTTTAACTGGTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCg                                                                                                              <  1:803159/140‑1 (MQ=255)        tatCTTTTATTTTTTAACTGGTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGtt                                                                                                         >  2:420462/1‑140 (MQ=255)           cTTTTATTTTTTAACTGGTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTa                                                                                                      <  2:263542/140‑1 (MQ=255)             tttATTTTTTAACTGGTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTAtc                                                                                                    >  1:642899/1‑140 (MQ=255)             tttATTTTTTAACTGGTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTAtc                                                                                                    >  1:228015/1‑140 (MQ=255)               tATTTTTTAACTGGTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTAtctc                                                                                                  >  1:617969/1‑140 (MQ=255)                    tttAACTGGTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCcac                                                                                             >  1:174136/1‑140 (MQ=255)                            gTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGcc                                                                                               >  2:526907/1‑130 (MQ=255)                            gTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGcc                                                                                               <  1:526907/130‑1 (MQ=255)                                  tAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGc                                                                               <  2:210139/140‑1 (MQ=255)                                  tAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGc                                                                               <  2:16627/140‑1 (MQ=255)                                      tGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCAt                                                                           >  2:50709/1‑140 (MQ=255)                                      tGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCAt                                                                           >  1:436357/1‑140 (MQ=255)                                      tGCAATGGCAGCCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCACCACAGAGTAAATTTTGCTCAt                                                                           >  2:204326/1‑140 (MQ=255)                                           tGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTg                                                                      >  1:698707/1‑140 (MQ=255)                                           tGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTg                                                                      >  1:588212/1‑140 (MQ=255)                                            ggCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGa                                                                     <  1:496823/140‑1 (MQ=255)                                                gTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGAtt                                                                       <  1:584483/134‑1 (MQ=255)                                                gTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGAtt                                                                       >  2:584483/1‑134 (MQ=255)                                                gTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGc                                                                 <  1:579417/140‑1 (MQ=255)                                                gTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGc                                                                 <  1:175327/140‑1 (MQ=255)                                                 tCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGc                                                                 <  1:442494/139‑1 (MQ=255)                                                 tCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGc                                                                 >  2:442494/1‑139 (MQ=255)                                                     tGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGt                                                            >  2:818599/1‑140 (MQ=255)                                                       aaCAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGc                                                                 <  1:599446/133‑1 (MQ=255)                                                       aaCAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGc                                                                 >  2:599446/1‑133 (MQ=255)                                                       aaCAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGttt                                                          >  2:385668/1‑140 (MQ=255)                                                        aCAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTa                                                         <  2:588212/140‑1 (MQ=255)                                                          aGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACg                                                       <  2:698707/140‑1 (MQ=255)                                                          aGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACg                                                       <  2:228015/140‑1 (MQ=255)                                                           ggCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGc                                                      <  1:424322/140‑1 (MQ=255)                                                               tctcTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGt                                                   <  1:493997/139‑1 (MQ=255)                                                               tctcTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGt                                                   >  2:493997/1‑139 (MQ=255)                                                                ctctTGCCATAAAACTGTCATCACTCACCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTAt                                                 <  1:818599/140‑1 (MQ=255)                                                                   ttGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAg                                              <  1:112187/140‑1 (MQ=255)                                                                        aTAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAata                                         >  1:153252/1‑140 (MQ=255)                                                                                        tCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTg                         <  1:253480/140‑1 (MQ=255)                                                                                        tCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTg                         <  1:113261/140‑1 (MQ=255)                                                                                        tCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTg                         <  2:617969/140‑1 (MQ=255)                                                                                           tCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTg                         >  2:646015/1‑137 (MQ=255)                                                                                           tCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTg                         <  1:646015/137‑1 (MQ=255)                                                                                                  aaaTGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATc               <  1:156108/140‑1 (MQ=255)                                                                                                   aaTGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATcc              <  1:305312/140‑1 (MQ=255)                                                                                                     tGTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAg            >  2:353142/1‑140 (MQ=255)                                                                                                      gTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGt           <  1:293752/140‑1 (MQ=255)                                                                                                      gTTAAAAAAGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGt           >  2:691572/1‑140 (MQ=255)                                                                                                             aaGCCGTTGCTTTGGGGATAGCCCGGTAAGGCCGGAGTTTTATCTCACCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAAttt    <  1:204326/140‑1 (MQ=255)                                                                                                                           gggATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATttgt  >  2:461185/1‑128 (MQ=255)                                                                                                                           gggATAGCCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATttgt  <  1:461185/128‑1 (MQ=255)                                                                                                                                 |                                                                                                                          AAAGAAATATCTTTTATTTTTTAACTGGTTAGATAAATGCAATGGCAGTCACTGAACAGGCATCTCTTGCCATAAAACTGTCATCACTCATCTTGACAAATGTTAAAAAAGCCGTTGCTTTGGGGATAACCCGGTAAGGCCGGAGTTTTATCTCGCCACAGAGTAAATTTTGCTCATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGT  >  NC_000913/813210‑813459

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