Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,522,713	T→C	Q13Q (CAA→CAG)	hofM ←	DNA utilization protein HofM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,522,713	0	T	C	100.0%	147.1 / NA	42	Q13Q (CAA→CAG)	hofM	DNA utilization protein HofM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/29);  total (13/29)

ACGTTTTAGCCAGCTGCTGCGCATCAACAATCCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAA  >  NC_000913/3522577‑3522847                                                                                                                                         |                                                                                                                                       acGTTTTAGCCAGCTGCTGCGCATCAACAATCCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCtgctgt                                                                                                                                     >  2:442458/1‑140 (MQ=255)         gCCAGCTGCTGCGCATCAACAATCCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCaa                                                                                                                             <  1:698047/140‑1 (MQ=255)         gCCAGCTGCTGCGCATCAACAATCCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCaa                                                                                                                             <  1:670075/140‑1 (MQ=255)                    cgcTTCAAAAACCCGCCCCGTTTTAAAAAGGCCGTACTCAGGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGcc                                                                                                                  <  2:338979/140‑1 (MQ=255)                      cATCAACAATCCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAg                                                                                                                >  1:594500/1‑140 (MQ=255)                                cctcccATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCGACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATg                                                                                                      <  2:186587/137‑1 (MQ=255)                                ccGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATg                                                                                                      <  1:442458/140‑1 (MQ=255)                                 cGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGc                                                                                                     <  2:42449/140‑1 (MQ=255)                                 cGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGc                                                                                                     <  1:243805/140‑1 (MQ=255)                                   cccATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCa                                                                                                   >  1:297383/1‑140 (MQ=255)                                   cccATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCa                                                                                                   >  1:716508/1‑140 (MQ=255)                                      aTCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATtg                                                                                                <  1:612224/140‑1 (MQ=255)                                         tttGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGtga                                                                                             <  1:564581/140‑1 (MQ=255)                                             aTAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTgatgat                                                                                         <  2:83292/140‑1 (MQ=255)                                             aTAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTgatgat                                                                                         <  2:230039/140‑1 (MQ=255)                                             aTAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTgatgat                                                                                         <  2:594500/140‑1 (MQ=255)                                              tAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGAtc                                                                                        <  2:364099/140‑1 (MQ=255)                                              tAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGAtc                                                                                        <  2:457239/140‑1 (MQ=255)                                              tAATGTCGTTCGCCGGCGGCAACCGCCACCAGCGGTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGAtc                                                                                        <  1:783299/140‑1 (MQ=255)                                                  gTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCt                                                                                    >  1:760421/1‑140 (MQ=255)                                                     gTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTcc                                                                                     <  2:459131/136‑1 (MQ=255)                                                     gTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTcc                                                                                     >  1:459131/1‑136 (MQ=255)                                                             gCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTc                                                                         >  2:250097/1‑140 (MQ=255)                                                              cGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGAtc                                                                                        <  2:495141/124‑1 (MQ=255)                                                              cGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGAtc                                                                                        >  1:495141/1‑124 (MQ=255)                                                                   aCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTc                                                                         <  1:765956/134‑1 (MQ=255)                                                                   aCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTc                                                                         >  2:765956/1‑134 (MQ=255)                                                                   aCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTg                                                                   >  2:344429/1‑140 (MQ=255)                                                                   aCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTg                                                                   <  1:344429/140‑1 (MQ=255)                                                                   aCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTg                                                                   <  1:110104/140‑1 (MQ=255)                                                                       ccaccaGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACgg                                                               <  2:527979/140‑1 (MQ=255)                                                                         accaGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGt                                                             <  1:580244/140‑1 (MQ=255)                                                                                ttGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCaa                                                      <  1:17969/140‑1 (MQ=255)                                                                                ttGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCaa                                                      <  1:250097/140‑1 (MQ=255)                                                                                           aTTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCt                                           <  2:716508/140‑1 (MQ=255)                                                                                              cTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGc                                        <  2:454035/140‑1 (MQ=255)                                                                                                          gCACGATCGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTAcc                            >  2:307062/1‑140 (MQ=255)                                                                                                                 cGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtt                     <  1:798153/140‑1 (MQ=255)                                                                                                                 cGCAACCGCTACCGCTTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtt                     >  2:499927/1‑140 (MQ=255)                                                                                                                      ccgctaccgctTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtt                     >  2:505649/1‑135 (MQ=255)                                                                                                                      ccgctaccgctTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtt                     <  1:505649/135‑1 (MQ=255)                                                                                                                          taccgctTCTTGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATaa            >  1:588111/1‑140 (MQ=255)                                                                                                                                   ttGCTGCTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCaa   <  1:499927/140‑1 (MQ=255)                                                                                                                                    tgctgcTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCaaa  <  2:588111/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       ACGTTTTAGCCAGCTGCTGCGCATCAACAATCCGCCCATCTTTGATAATGTCGTTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAA  >  NC_000913/3522577‑3522847

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: ‑