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

CGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGG  >  NZ_CP009273/4153742‑4153894
                                                                          |                                                                              
cGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTc                                                                 <  2:388497/90‑1 (MQ=255)
    cTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCgg                                                             <  1:436895/90‑1 (MQ=255)
      cGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCCTCTTTTTTTTTTTGCggga                                                           >  2:190005/1‑88 (MQ=255)
            atgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTTAGGGCAGCTCCCCTCTATTTTTTTTTGCGGGACAAATg                                                     >  2:373529/1‑90 (MQ=255)
                 gtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCCGCTCCCCTCTTTTTTTTTTTGCGGTACAAAAGCCAGt                                                >  2:404001/1‑90 (MQ=255)
                     tgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCAtg                                            <  1:404001/90‑1 (MQ=255)
                      gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTTGCGGGACAAATGCCAGTCAtgt                                           >  1:355448/1‑90 (MQ=255)
                         ccGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGtt                                        <  1:67138/90‑1 (MQ=255)
                          cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                       <  2:359474/90‑1 (MQ=255)
                                     cGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTTTTCGCGGGAAAAAAGCCAGTCATGTGTTGGTGGTTATTGa                            >  1:478835/1‑90 (MQ=255)
                                      gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                           <  2:345616/90‑1 (MQ=255)
                                                cccAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGcc                 >  1:13934/1‑90 (MQ=255)
                                                        ggcggTTCAGGTCAGCTCCCCTCTATTTTTTTTCGCGGGACAAAAGCCCGTCATGTGTTGGGGTTAATTTATGGCGTACGCCTGAATCTg         >  2:423899/1‑90 (MQ=255)
                                                            gTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCgg     <  1:423899/90‑1 (MQ=255)
                                                               cAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCggggg  <  2:478835/90‑1 (MQ=255)
                                                                      gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAAt            <  2:145954/73‑1 (MQ=255)
                                                                      gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAAt            >  1:145954/1‑73 (MQ=255)
                                                                      gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCggg    <  2:13520/81‑1 (MQ=255)
                                                                      gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCggg    >  1:13520/1‑81 (MQ=255)
                                                                          |                                                                              
CGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGG  >  NZ_CP009273/4153742‑4153894

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