Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,173,777	T→G	N232K (AAT→AAG)	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,777	0	T	G	100.0%	183.8 / NA	49	N232K (AAT→AAG)	birA	bifunctional biotin‑[acetylCoA carboxylase] holoenzyme synthetase/ DNA‑binding transcriptional repressor, bio‑5'‑AMP‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/27);  total (22/27)

GGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTT  >  NC_000913/4173633‑4173915                                                                                                                                                 |                                                                                                                                           ggcaggcaTTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGAt                                                                                                                                      <  1:790808/151‑1 (MQ=255)  gcaagcaTTCTGGTGGAGCTGACTGGCAAAAGTGGCGATGCGGAGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGCATCACGCTGCAGGAAGCGGGGATCAAGCGCGATc                                                                                                                                     <  1:366349/151‑1 (MQ=255)       caTTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAAt                                                                                                                                <  1:57633/151‑1 (MQ=255)          tCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATAc                                                                                                                              <  1:442595/150‑1 (MQ=255)             ggtggAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACg                                                                                                                             <  2:163113/148‑1 (MQ=255)                 gAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTggcgg                                                                                                                      <  2:78876/151‑1 (MQ=255)                 gAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTggcgg                                                                                                                      >  1:595035/1‑151 (MQ=255)                     tGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCAt                                                                                                                  >  1:463663/1‑151 (MQ=255)                      gACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATg                                                                                                                 <  1:469023/151‑1 (MQ=255)                          ggCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCt                                                                                                               <  1:340944/149‑1 (MQ=255)                            cAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAAt                                                                                                            >  1:233776/1‑150 (MQ=255)                            cAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAAt                                                                                                            >  2:566842/1‑150 (MQ=255)                             aaaaCTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGTGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATAc                                                                                                          >  1:182666/1‑151 (MQ=255)                             aaaaCTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATAc                                                                                                          <  2:565753/151‑1 (MQ=255)                               aaCTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACg                                                                                                         >  1:699732/1‑150 (MQ=255)                                           ggcgCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGtgct                                                                                            >  2:23760/1‑151 (MQ=255)                                              gcAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCg                                                                                         <  1:770836/151‑1 (MQ=255)                                               cAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGt                                                                                        >  2:799333/1‑151 (MQ=255)                                                 aaTAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTg                                                                                      <  2:167633/151‑1 (MQ=255)                                                  aTAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTgg                                                                                     >  2:780681/1‑151 (MQ=255)                                                        aTTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAActct                                                                               >  2:403656/1‑151 (MQ=255)                                                           ggAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTc                                                                             <  2:248823/150‑1 (MQ=255)                                                                cGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACaa                                                                       <  2:544116/151‑1 (MQ=255)                                                                cGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACa                                                                        >  2:931743/1‑150 (MQ=255)                                                                      cAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAgg                                                                  >  2:634181/1‑150 (MQ=255)                                                                         cATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTg                                                              <  2:283227/151‑1 (MQ=255)                                                                         cATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTg                                                              >  2:348015/1‑151 (MQ=255)                                                                          aTGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTgg                                                             <  2:483637/151‑1 (MQ=255)                                                                                    gcCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCt                                                   <  1:178792/151‑1 (MQ=255)                                                                                       gtgtTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTc                                                <  2:63379/151‑1 (MQ=255)                                                                                          ttGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTcgcg                                             >  2:708070/1‑151 (MQ=255)                                                                                           tGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTcgcg                                             <  1:544312/150‑1 (MQ=255)                                                                                            gAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCt                                           >  2:782257/1‑151 (MQ=255)                                                                                            gAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCt                                           <  1:721108/151‑1 (MQ=255)                                                                                               gagagTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTggg                                        >  2:785965/1‑151 (MQ=255)                                                                                                agagTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGa                                       <  1:172167/151‑1 (MQ=255)                                                                                                   gtgtCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGaaaa                                    <  2:145913/151‑1 (MQ=255)                                                                                                   gtgtCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGaaaa                                    <  2:202836/151‑1 (MQ=255)                                                                                                    tgtCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGaaa                                     >  2:436606/1‑149 (MQ=255)                                                                                                     gtCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGc                                  >  1:97637/1‑151 (MQ=255)                                                                                                            aTCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATa                            >  2:120685/1‑150 (MQ=255)                                                                                                              cAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAAtt                         >  2:162755/1‑151 (MQ=255)                                                                                                                ggggTGGATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAAtttt                       <  2:626232/151‑1 (MQ=255)                                                                                                                     ggATCACGCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTaa                  <  1:114821/151‑1 (MQ=255)                                                                                                                            gCTGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCa           <  2:727002/151‑1 (MQ=255)                                                                                                                              tGCAGGAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGt         >  1:572434/1‑151 (MQ=255)                                                                                                                                  ggAAGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGaaa     <  2:905892/151‑1 (MQ=255)                                                                                                                                     aGCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACtt  <  1:365978/151‑1 (MQ=255)                                                                                                                                      gCGGGGATCAAGCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACt   >  1:328631/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                           GGCAGGCATTCTGGTGGAGCTGACTGGCAAAACTGGCGATGCGGCGCAAATAGTCATTGGAGCCGGGATCAACATGGCAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTT  >  NC_000913/4173633‑4173915

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