Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,173,601	C→T	P174S (CCT→TCT)	birA →	bifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase/ DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,173,601	0	C	T	100.0%	192.3 / NA	50	P174S (CCT→TCT)	birA	bifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase/ DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/19);  total (31/19)

TCGCCTTTTGGCGCAAACTTATATTTGTCGATGTTCTGGCGTCTGGAACAAGGCCCGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGG  >  NC_000913/4173454‑4173747                                                                                                                                                    |                                                                                                                                                   tCGCCTTTTGGCGCAAACTTATATTTGTCGATGTTCTGGCGTCTGGAACAAGGCCCGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCt                                                                                                                                                  <  1:457579/150‑1 (MQ=255)       tttGGCGCAAACTTATATTTGTCGATGTTCTGGCGTCTGGAACAAGGCCCGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGa                                                                                                                                             <  2:169536/149‑1 (MQ=255)          ggCGCAAACTTATATTTGTCGATGTTCTGGCGTCTGGAACAAGGCCCGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACctc                                                                                                                                         <  1:195952/150‑1 (MQ=255)           gcgcAAACTTATATTTGTCGATGTTCTGGCGTCTGGAACAAGGCCCGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACctct                                                                                                                                        >  2:216981/1‑150 (MQ=255)               aaaCTTATATTTGTCGATGTTCTGGCGTCTGGAACAAGGCCCGGCGGCGGCGATTGGTTTAAGTCAGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCt                                                                                                                                    <  2:240532/150‑1 (MQ=255)                aaCTTATATTTGTCGATGTTCTGGCGTCTGGAACAAGGCCCGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGc                                                                                                                                  <  2:577330/151‑1 (MQ=255)                aaCTTATATTTGTCGATGTTCTGGCGTCTGGAACAAGGCCCGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGc                                                                                                                                  >  2:456387/1‑151 (MQ=255)                             cGATGTTCTGGCGTCTGGAACAAGGCCCGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCt                                                                                                                     >  1:621627/1‑151 (MQ=255)                               aTGTTCTGGAGTCTGGAACAAGGCCCGGCGGCGGCGATTGGTTTATGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGATCGTGTTAAATGGTCTAATGAACTCTATCTGCAGGATCGCCAGCTgg                                                                                                                   >  2:402821/1‑151 (MQ=255)                                            tGGAACAAGGCCCGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCtggt                                                                                                      >  1:309863/1‑151 (MQ=255)                                                      cccGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACt                                                                                            <  1:273419/151‑1 (MQ=255)                                                       ccGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTg                                                                                           >  1:558832/1‑151 (MQ=255)                                                        cggcggcggCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTgg                                                                                          <  2:430792/151‑1 (MQ=255)                                                           cggcggcgATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCa                                                                                        >  2:273273/1‑150 (MQ=255)                                                           cggcggcgATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCa                                                                                        >  2:334655/1‑150 (MQ=255)                                                            ggcggcgATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCaaa                                                                                      <  1:428242/151‑1 (MQ=255)                                                                  gATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGc                                                                                <  1:21891/151‑1 (MQ=255)                                                                    ttGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCg                                                                               >  1:64138/1‑150 (MQ=255)                                                                      ggTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATg                                                                            >  1:244864/1‑151 (MQ=255)                                                                      ggTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATg                                                                            <  1:145834/151‑1 (MQ=255)                                                                        tttAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATgcg                                                                          >  2:188559/1‑151 (MQ=255)                                                                        tttAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATgcg                                                                          >  2:140384/1‑151 (MQ=255)                                                                          tAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATgcgg                                                                         <  1:87694/150‑1 (MQ=255)                                                                            aGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGgcgc                                                                      >  1:213214/1‑151 (MQ=255)                                                                                 ggTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGTAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATa                                                                 <  2:659735/151‑1 (MQ=255)                                                                                 ggTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATa                                                                 <  2:359372/151‑1 (MQ=255)                                                                                        ggTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTg                                                          >  1:246515/1‑151 (MQ=255)                                                                                         gTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTgg                                                         <  2:632864/151‑1 (MQ=255)                                                                                           aTCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTg                                                          >  2:432825/1‑148 (MQ=255)                                                                                              gTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCg                                                    >  1:221120/1‑151 (MQ=255)                                                                                                gATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCggg                                                  <  2:185511/151‑1 (MQ=255)                                                                                                gATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCggg                                                  >  1:128459/1‑151 (MQ=255)                                                                                                         aGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATg                                         >  1:99856/1‑151 (MQ=255)                                                                                                         aGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATg                                         <  2:606108/151‑1 (MQ=255)                                                                                                            aTTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGc                                       >  1:295591/1‑150 (MQ=255)                                                                                                             ttACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCa                                      >  2:286397/1‑150 (MQ=255)                                                                                                                cGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATgc                                  >  1:398488/1‑151 (MQ=255)                                                                                                                   aaGCTGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGcc                               >  2:546670/1‑151 (MQ=255)                                                                                                                       tGGGTGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGTATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCgtgt                           >  1:51012/1‑151 (MQ=255)                                                                                                                           tGCAGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGa                        >  1:561982/1‑150 (MQ=255)                                                                                                                              aGATAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAagag                    <  1:103926/151‑1 (MQ=255)                                                                                                                                 tAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAagagag                  >  1:216422/1‑150 (MQ=255)                                                                                                                                 tAAAGTTCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAgt                 >  2:278180/1‑151 (MQ=255)                                                                                                                                       tCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCg                                                    >  1:250614/1‑110 (MQ=255)                                                                                                                                       tCGTGTTAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCg                                                    <  2:250609/110‑1 (MQ=255)                                                                                                                                            ttAAATGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATc       >  1:223911/1‑150 (MQ=255)                                                                                                                                              aaaTGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAgg    >  1:602796/1‑151 (MQ=255)                                                                                                                                              aaaTGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAgg    >  2:547052/1‑151 (MQ=255)                                                                                                                                                aTGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAgggg  >  2:219641/1‑151 (MQ=255)                                                                                                                                                aTGGTCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAgggg  <  1:593844/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TCGCCTTTTGGCGCAAACTTATATTTGTCGATGTTCTGGCGTCTGGAACAAGGCCCGGCGGCGGCGATTGGTTTAAGTCTGGTTATCGGTATCGTGATGGCGGAAGTATTACGCAAGCTGGGTGCAGATAAAGTTCGTGTTAAATGGCCTAATGACCTCTATCTGCAGGATCGCAAGCTGGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGG  >  NC_000913/4173454‑4173747

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