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_0077793,472,5580GC100.0% 128.6 / NA 31A162G (GCA→GGA) btuBvitamin 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: