Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	178,175	C→A	Q96H (CAG→CAT)	btuF ←	vitamin B12 ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	178,175	0	C	A	100.0%	95.1 / NA	29	Q96H (CAG→CAT)	btuF	vitamin B12 ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (13/16);  total (13/16)

GATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTGCGCTTGTGGAGG  >  NC_000913/178050‑178303                                                                                                                              |                                                                                                                                 gATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCg                                                                                                                    >  2:591189/1‑140 (MQ=255)  aTTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGc                                                                                                                   <  2:534703/140‑1 (MQ=255)  aTTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGc                                                                                                                   <  2:29192/140‑1 (MQ=255)   ttGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCt                                                                                                                  >  2:15544/1‑140 (MQ=255)     gcgcGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCg                                                                                                                >  2:664103/1‑140 (MQ=255)      cgcgGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTACCGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCgg                                                                                                               <  2:304809/140‑1 (MQ=255)      cgcgGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCgg                                                                                                               <  1:135038/140‑1 (MQ=255)         ggCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCAt                                                                                                            <  1:196871/140‑1 (MQ=255)               ttCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTc                                                                                                      <  1:591189/140‑1 (MQ=255)                    ccTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGc                                                                                                 <  1:474896/140‑1 (MQ=255)                         tCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGc                                                                                            >  2:5193/1‑140 (MQ=255)                          cTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCa                                                                                           <  2:463575/140‑1 (MQ=255)                          cTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCa                                                                                           <  2:260103/140‑1 (MQ=255)                               tGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCAc                                                                                      <  2:612662/140‑1 (MQ=255)                                 cGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAACCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACACCAGGCAATCACCa                                                                                    >  2:47435/1‑140 (MQ=255)                                                cAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAgcg                                                                     >  1:272644/1‑140 (MQ=255)                                                cAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAgcg                                                                     >  1:143789/1‑140 (MQ=255)                                                  gTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAGCCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGcgcg                                                                   <  2:653139/140‑1 (MQ=255)                                                  gTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGcgcg                                                                   >  2:493273/1‑140 (MQ=255)                                                      aCGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACaa                                                               <  1:493273/140‑1 (MQ=255)                                                            cGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGtt                                                         >  1:174973/1‑140 (MQ=255)                                                                                      tCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGaa                               <  1:138983/140‑1 (MQ=255)                                                                                          aTCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTg                                  >  2:154089/1‑133 (MQ=255)                                                                                          aTCGACCCACATCACTTTTATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTg                                  <  1:154089/133‑1 (MQ=255)                                                                                                         ttttATTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTATCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTgcg            <  1:664103/140‑1 (MQ=255)                                                                                                             aTTCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTGCGCTtg        >  1:150105/1‑140 (MQ=255)                                                                                                              ttCCCAGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTGCGCTtgt       >  2:47308/1‑140 (MQ=255)                                                                                                                   aGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTGCGCTTGTggagg  >  1:89744/1‑140 (MQ=255)                                                                                                                   aGCGAAGCCAGATGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTGCGCTTGTggagg  <  2:62567/140‑1 (MQ=255)                                                                                                                              |                                                                                                                                 GATTGCGCGGCTTGTTCGGCCTTGTCTGGTTGCGGACTCCAGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTGCGCTTGTGGAGG  >  NC_000913/178050‑178303

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