Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007779	3,472,558	G→C	A162G (GCA→GGA)	btuB ←	vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	3,472,558	0	G	C	100.0%	128.6 / NA	31	A162G (GCA→GGA)	btuB	vitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (10/21);  total (10/21)

ATGGGCATAATCGCCCAACAGCGTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTGCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGA  >  NC_007779/3472461‑3472654                                                                                                  |                                                                                                 atcagCATAATCGCCCAACAGCGTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCctcctg                                                                                               >  2:60482/5‑101 (MQ=255) aTGGGCATAATCGCCCAACAGCGTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCctcctg                                                                                               >  1:428935/1‑101 (MQ=255) aTGGGCATAATCGCCCAACAGCGTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCctcctg                                                                                               >  1:839871/1‑101 (MQ=255)    ggCATAATCGCCCAACAGCGTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGaaa                                                                                            >  3:717732/1‑101 (MQ=255)      cATAATCGCCCAACAGCGTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAAtt                                                                                          <  1:69640/101‑1 (MQ=255)          aTCGCCCAACAGCGTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATttccg                                                                                      <  3:672626/101‑1 (MQ=255)             gCCCAACAGCGTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATttccgttc                                                                                   <  1:101667/101‑1 (MQ=255)             gCCCAACAGCGTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATttccgttc                                                                                   <  1:852045/101‑1 (MQ=255)                 aaCAGCGTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCggg                                                                               <  1:199120/101‑1 (MQ=255)                       gTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATc                                                                         <  3:464686/101‑1 (MQ=255)                       gTTACCCGTGTCTTATCCCCCAGTTATTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATc                                                                         >  3:1584754/1‑101 (MQ=255)                        ttACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATcg                                                                        >  3:84072/1‑101 (MQ=255)                         tACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATcgc                                                                       <  1:1373641/101‑1 (MQ=255)                          aCCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATcgcg                                                                      >  3:973349/1‑101 (MQ=255)                           cccGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATcgcgc                                                                     <  3:1550596/101‑1 (MQ=255)                                 tCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGtg                                                               <  1:43006/101‑1 (MQ=255)                                  cTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGtga                                                              >  3:992335/1‑101 (MQ=255)                                   ttATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGtgat                                                             <  3:895591/101‑1 (MQ=255)                                      tCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTgatgat                                                          <  1:312516/101‑1 (MQ=255)                                        ccccAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGatat                                                        <  3:1256966/101‑1 (MQ=255)                                                gttgCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCaccacc                                                >  1:415849/1‑101 (MQ=255)                                                 ttgCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCAccc                                               >  1:1332047/1‑101 (MQ=255)                                                  tgctgcGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCCCCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCAcccc                                              <  1:511135/101‑1 (MQ=255)                                                      gcGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCt                                          >  1:563998/1‑101 (MQ=255)                                                                    tAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAAcc                            <  1:118064/101‑1 (MQ=255)                                                                    tAGTTCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAAcc                            <  3:1261052/101‑1 (MQ=255)                                                                       ttCTGATAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATa                         >  1:936610/1‑101 (MQ=255)                                                                             tAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGc                   >  3:706405/1‑101 (MQ=255)                                                                             tAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGc                   <  1:165190/101‑1 (MQ=255)                                                                             tAACTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGc                   <  2:82911/101‑1 (MQ=255)                                                                               aCTATTGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCCCCCCGCCTATTGCATCGGAACCATAAACAGCgg                 <  3:688498/101‑1 (MQ=255)                                                                                   ttGCTTCCCCACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAgcg             <  1:982740/101‑1 (MQ=255)                                                                                            cACCCTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACg    <  3:1060064/101‑1 (MQ=255)                                                                                              cccTCCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGa  <  1:91981/101‑1 (MQ=255)                                                                                                  |                                                                                                 ATGGGCATAATCGCCCAACAGCGTTACCCGTGTCTTATCCCCCAGTTGTTGCTGCGTAGAGACATCATAGTTCTGATAACTATTGCTTCCCCACCCTGCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGA  >  NC_007779/3472461‑3472654

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