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,8160TC88.9% 16.9 / ‑3.3 9S84P (TCA→CCA) btuBTonB‑dependent vitamin B12 receptor BtuB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (4/4);  total (4/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

ACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG  >  NZ_CP009273/4153744‑4153901
                                                                        |                                                                                     
aCCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                                      <  1:213375/90‑1 (MQ=255)
              tgtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCa                                                                                     <  2:169112/61‑1 (MQ=255)
              tgtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAg                                                        <  2:144053/90‑1 (MQ=255)
                 gctgaGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCa                                                     <  2:238635/90‑1 (MQ=255)
                    gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCCGGGCAGCTCCCCTCTATTTTTTTTTTCGGTACAAAAGCCAGTCAtgt                                                  >  1:347011/1‑90 (MQ=255)
                    gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGGACAAATGCCCGTCAtgt                                                  >  2:200673/1‑90 (MQ=255)
                        cGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTg                                              <  1:142279/90‑1 (MQ=255)
                         gTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTATTCGCGGGACAAATGGCAGTCATGTGTTgg                                             >  2:406903/1‑90 (MQ=255)
                                                                    gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTgg  >  2:126255/1‑90 (MQ=255)
                                                                        |                                                                                     
ACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG  >  NZ_CP009273/4153744‑4153901

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