Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,153,816 T→C S84P (TCA→CCA)  btuB → TonB‑dependent vitamin B12 receptor BtuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,153,8160TC100.0% 54.6 / NA 18S84P (TCA→CCA) btuBTonB‑dependent vitamin B12 receptor BtuB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (5/13);  total (5/13)

GACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGG  >  NZ_CP009273/4153729‑4153894
                                                                                       |                                                                              
gACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCa                                                                              <  2:295797/90‑1 (MQ=255)
gACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGGTCAGGTCAGCTCCCa                                                                              >  2:155965/1‑90 (MQ=255)
             cGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTc                                                                 <  2:392083/90‑1 (MQ=255)
               aCCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                               <  1:108281/90‑1 (MQ=255)
                         atgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                               <  1:27369/80‑1 (MQ=255)
                         atgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                               >  2:27369/1‑80 (MQ=255)
                          tgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGc                                                    <  1:491641/90‑1 (MQ=255)
                           caggtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                   <  1:85132/87‑1 (MQ=255)
                                 ctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCAt                                             <  1:398415/90‑1 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCAtgt                                           >  1:403374/1‑90 (MQ=255)
                                       cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                       <  1:473405/90‑1 (MQ=255)
                                        gTCTTCCGGGCGTCGATATCACCCAAAACGGCGGGTCAGGTCACGCCCCCCCTTTTTTTTTTCGCGGTACAAAAGCCACGCATGTGtttg                                      >  2:85132/1‑88 (MQ=255)
                                                 gCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAAAGCCAGTCATGTGTTGGGGGTAATTg                             >  1:448635/1‑90 (MQ=255)
                                                   gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                           <  2:109213/90‑1 (MQ=255)
                                                   gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                           <  1:370498/90‑1 (MQ=255)
                                                         atCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTa                     <  1:186483/90‑1 (MQ=255)
                                                                         gTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCgg     <  1:160979/90‑1 (MQ=255)
                                                                            cAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCggggg  <  2:168868/90‑1 (MQ=255)
                                                                                       |                                                                              
GACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGG  >  NZ_CP009273/4153729‑4153894

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: