Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,136,644 T→C T153T (ACT→ACC pflC → [formate‑C‑acetyltransferase]‑activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,136,6440TC95.2% 64.6 / ‑3.2 21T153T (ACT→ACCpflC[formate‑C‑acetyltransferase]‑activating enzyme
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (15/5);  total (16/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 = 9.58e-01

GCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGCGATGAAGTGTTGTTCGATTTAAAAATTATGGAC  >  NZ_CP009273/4136562‑4136725
                                                                                  |                                                                                 
gcgAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACCGCCGgag                                                                            >  1:174145/1‑90 (MQ=255)
               aGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCa                                                             >  1:222195/1‑90 (MQ=255)
               aGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCa                                                             >  1:122780/1‑90 (MQ=255)
               aGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGGGTGTCATGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCa                                                             >  2:362457/1‑90 (MQ=255)
                            tACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCt                                                <  1:14184/90‑1 (MQ=255)
                              cccGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTgg                                              >  2:311632/1‑90 (MQ=255)
                              cccGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTgg                                              >  2:75179/1‑90 (MQ=255)
                                        aCAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATtgtg                                    <  2:77075/90‑1 (MQ=255)
                                                                gtCATGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCAAg                                                           <  1:291512/43‑1 (MQ=255)
                                                                gtCATGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCAAg                                                           >  2:291512/1‑43 (MQ=255)
                                                                 tCATGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGCGATGAAGTGTTGTTCGATTTaaaa           >  1:175325/1‑90 (MQ=255)
                                                                   aTGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGCGATGAAGTGTTGTTCGATTTAAAAAt         >  1:382044/1‑90 (MQ=255)
                                                                   aTGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGCGATGAAGTGTTGTTCGATTTAAAAAt         >  2:517491/1‑90 (MQ=255)
                                                                   aTGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGCGATGAAGTGTTGTTCGATTTAAAAAt         >  1:136895/1‑90 (MQ=255)
                                                                    tGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGc                                   >  2:252496/1‑63 (MQ=255)
                                                                    tGCGCCATTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGc                                   <  1:252496/63‑1 (MQ=255)
                                                                     gcgcCATTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGCGATGAAGTGTTGTTCGATTTAAAAATTa       >  2:270684/1‑90 (MQ=255)
                                                                        ccATTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGCGATGAAGTGTTGTTCGATTTAAAAATTATgg    >  2:481996/1‑90 (MQ=255)
                                                                          aTTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGCGATGAAGTGTTGTTCGATTTAAAAATTATGGAc  <  1:362457/90‑1 (MQ=255)
                                                                          aTTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGCGATGAAGTGTTGTTCGATTTAAAAATTATGGAc  >  2:387642/1‑90 (MQ=255)
                                                                          aTTGAAACCGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGCGATGAAGTGTTGTTCGATTTAAAAATTATGGAc  >  1:201898/1‑90 (MQ=255)
                                                                                  |                                                                                 
GCGAAGTGTTAATGCAGGCGGAGTTTGCTACCCGTTTTTTACAGCGACTGCGGCTGTGGGGTGTGTCATGCGCCATTGAAACTGCCGGAGACGCACCAGCCAGCAAGCTATTACCGCTGGCGAAATTGTGCGATGAAGTGTTGTTCGATTTAAAAATTATGGAC  >  NZ_CP009273/4136562‑4136725

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