Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,791,751	A→G	R267R (CGT→CGC)	kbl ←	2‑amino‑3‑ketobutyrate CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,791,751	0	A	G	100.0%	97.1 / NA	28	R267R (CGT→CGC)	kbl	2‑amino‑3‑ketobutyrate CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/13);  total (15/13)

TCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCA  >  NC_000913/3791620‑3791883                                                                                                                                    |                                                                                                                                     tCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGgcg                                                                                                                              <  1:362872/140‑1 (MQ=255)     ggAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGc                                                                                                                          <  2:290824/140‑1 (MQ=255)        aCTGACGCGCGTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCAc                                                                                                                       <  1:127379/140‑1 (MQ=255)        aCTGACGCGCGTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCAc                                                                                                                       <  1:531123/140‑1 (MQ=255)                             gACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGcgcc                                                                                                  <  2:254698/140‑1 (MQ=255)                                 gTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCgg                                                                                              >  1:352073/1‑140 (MQ=255)                                  tCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGt                                                                                             >  1:544884/1‑140 (MQ=255)                                             gCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAg                                                                                  <  1:276226/140‑1 (MQ=255)                                                        cGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGccc                                                                       <  2:94140/140‑1 (MQ=255)                                                            cATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGTCCCGCCCAgcg                                                                   <  2:544884/140‑1 (MQ=255)                                                                    gTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCa                                                           >  1:126979/1‑140 (MQ=255)                                                                          tGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTa                                                     <  1:61163/140‑1 (MQ=255)                                                                          tGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTa                                                     >  2:401354/1‑140 (MQ=255)                                                                          tGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTa                                                     >  2:163886/1‑140 (MQ=255)                                                                          tGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTa                                                     >  2:362287/1‑140 (MQ=255)                                                                                 cgcgGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGa                                              >  1:31928/1‑140 (MQ=255)                                                                                 cgcgGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGCACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGa                                              >  2:519415/1‑140 (MQ=255)                                                                                   cgGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGata                                            >  1:201222/1‑140 (MQ=255)                                                                                            tGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGAcc                                   <  1:290467/140‑1 (MQ=255)                                                                                                  gTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCgg                                <  1:248391/137‑1 (MQ=255)                                                                                                  gTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCgg                                >  2:248391/1‑137 (MQ=255)                                                                                                   tGCCAGCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCa                            >  2:541095/1‑140 (MQ=255)                                                                                                        gCGAGTTGGAGAACAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATACCGCCCCGGCCCATcaca                       <  2:278936/140‑1 (MQ=255)                                                                                                                     cAGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCAt          >  2:399019/1‑140 (MQ=255)                                                                                                                      aGGTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATg         <  2:126979/140‑1 (MQ=255)                                                                                                                        gTACGGACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATggg       >  1:300121/1‑140 (MQ=255)                                                                                                                            ggACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTATCAAGCGTACCGGTGATAATATCGACCCGGTCCATCACATCGCAGTATTCCTGGGAAcc   >  1:272334/1‑140 (MQ=255)                                                                                                                             gACGAGAGCGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCa  >  1:511581/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                     TCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCA  >  NC_000913/3791620‑3791883

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