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% 65.6 / NA 20S84P (TCA→CCA) btuBTonB‑dependent vitamin B12 receptor BtuB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/11);  total (9/11)

GACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG  >  NZ_CP009273/4153729‑4153901
                                                                                       |                                                                                     
gACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCa                                                                                     <  1:521715/90‑1 (MQ=255)
                    ggTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTAc                                                                 <  1:106925/90‑1 (MQ=255)
                      tCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACaa                                                               <  1:433463/90‑1 (MQ=255)
                          tgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGc                                                           <  2:279509/90‑1 (MQ=255)
                           gatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                          <  2:782004/90‑1 (MQ=255)
                           gatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                          <  1:769288/90‑1 (MQ=255)
                             tgtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAg                                                        <  1:407309/90‑1 (MQ=255)
                                  tgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCCTCTTTTTTTATTCGCGGTACAAAAGCCCGTCAtg                                                   >  1:274282/1‑90 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTTAGGGCAGCTCCCCTCTATTTTTTTTTGCGGGACAAAAGCCAGTCAttt                                                  >  2:832273/1‑88 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGGACAAGTGCCAGTCAtgt                                                  >  1:913024/1‑90 (MQ=255)
                                   gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCCGCTCCCCTCTTTTTTTTTTCGCGGGACAAATGCCCGTCAtgt                                                  >  2:801656/1‑90 (MQ=255)
                                       cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                              <  2:2363/90‑1 (MQ=255)
                                                   gTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTTTTCGCGGGACAAATGCCAGCCATGTTTTGGGGTTAATTGAt                                  >  1:52342/1‑90 (MQ=255)
                                                   gTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCCGTCATGTGTTGGGGTTAATTGAt                                  >  1:557101/1‑90 (MQ=255)
                                                         atCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTa                            <  2:199685/90‑1 (MQ=255)
                                                                 aaaCGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGGGTTAATTGATGGCGTACGCCTGaa                    >  1:401586/1‑90 (MQ=255)
                                                                       cggTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGc              <  2:605078/90‑1 (MQ=255)
                                                                                   gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGt                                                       <  2:20636/37‑1 (MQ=38)
                                                                                   gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGt                                                       >  1:20636/1‑37 (MQ=38)
                                                                                   gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTgg  >  1:184459/1‑90 (MQ=255)
                                                                                       |                                                                                     
GACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG  >  NZ_CP009273/4153729‑4153901

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