Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273322,5020CG54.5% 2.7 / 13.5 11G153R (GGC→CGC) betAcholine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base G (0/6);  total (4/7)
Fisher's exact test for biased strand distribution p-value = 1.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGC  >  NZ_CP009273/322413‑322582
                                                                                         |                                                                                
ttCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGcc                                                                                  >  1:406185/1‑90 (MQ=255)
ttCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGcc                                                                                  <  2:263882/90‑1 (MQ=255)
ttCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGcc                                                                                  <  2:44672/90‑1 (MQ=255)
ttCCTGCTAATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGcc                                                                                  >  2:166669/1‑90 (MQ=255)
     gctgATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGtt                                                                             >  1:348885/1‑90 (MQ=255)
                       cGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCa                                                           >  1:106254/1‑90 (MQ=255)
                         tCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCAcc                                                         <  1:11415/90‑1 (MQ=255)
                         tCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCAcc                                                         >  1:241379/1‑90 (MQ=255)
                                             cTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATa                                     >  2:366320/1‑90 (MQ=255)
                                                    cccGCTTCAACCCCCGTTCCAAACCGCGGTTTGCGGCGGGGTTTGGGGGTAGTGCCCCCCCCCGGCCCCTCCCCGCCGTGATAGTCGttt                              <  2:209622/90‑1 (MQ=255)
                                                                   cctccAAACAGCGGTTTGACGCGGGGTTTGGAGGTAGTGCCGCCCCCCGGGCCATCTCCGCCGTGATAGTCGTTTTCACCCATATCGCGa               <  2:348885/86‑1 (MQ=255)
                                                                        aaaCAGCGGTTTGCCGCGGGGTTTGGGGGTAGTGACGCCCCCCGGGCCACCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGtctc          <  2:241379/90‑1 (MQ=255)
                                                                              cGGATTGACGCGGGGTTTGGAGGTAGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCtt    <  1:533865/90‑1 (MQ=255)
                                                                                gATTGGCGCGGGGTTGGGGGGTGGTGCCGCCCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc  <  2:579438/90‑1 (MQ=255)
                                                                                gATTGACGCGGGGTTGGGGGGTGGTGACGCTCCCCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGc  <  2:39596/90‑1 (MQ=255)
                                                                                         |                                                                                
TTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGTGACGCTCACCGGGCCATCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGC  >  NZ_CP009273/322413‑322582

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