Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,102,2550TC53.8% 1.5 / 13.6 13E87G (GAG→GGG) speCornithine decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (6/0);  new base C (0/7);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 5.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGTTGCCCATCTCAACGTACTGCGTCAGCGTGTCATAAAACGGTGGCAGCAAATTCTCTTCATACTGACAGGCTGCGGATTCCAGCTCCAGCCACTGCTGCTCGTTGCCGTTGATTACCGCCGTAACGCCCGCAGGTAATTCAACAGCATGTTCGGAATACAAAAACACCGGTA  >  NZ_CP009273/3102167‑3102342
                                                                                        |                                                                                       
gCTGTTGCCCATCTCAACGTACTGCGTCAGCGTGTCATAAAACGGTGGCAGCAAATTCTCTTCATACTGACAGGCTGCGGATtccagctc                                                                                        <  1:7778/90‑1 (MQ=255)
                     cTGCGTCAGCGTGTCATAAAACGGTGGCAGCAAATTCTCTTCATACTGACAGGCTGCGGATTCCAGCTCCAGCCACTGCTGCTCGTTGcc                                                                   >  1:286171/1‑90 (MQ=255)
                      tGCGTCAGCGTGTCATAAAACGGTGGCAGCAAATTCTCTTCATACTGACAGGCTGCGGATTCCAGCTCCAGCCACTGCTGCTCGTTGCCg                                                                  >  1:101893/1‑90 (MQ=255)
                      tGCGTCAGCGTGTCATAAAACGGTGGCAGCAAATTCTCTTCATACTGACAGGCTGCGGATTCCAGCTCCAGCCACTGCTGCTCGTTGCCg                                                                  >  1:66109/1‑90 (MQ=255)
                            aGCGTGTCATAAAACGGTGGCAGCAAATTCTCTTCATACTGACAGGCTGCGGATTCCAGCTCCAGCCACTGCTGCTCGTTGCCGTTGAtt                                                            >  1:69181/1‑90 (MQ=255)
                                        aaCGGTGGCAGCAATTTCTCTTCATCCTGACAGGCTGCGGATTCCAGCCCCAGCCACTGCTGCTCGTTGCCGTTGATTACCGCCGTAACg                                                <  2:202051/90‑1 (MQ=255)
                                            gTGGCAGCAAATTCTCTTCATACTGACAGGCTGCGGATTCCAGCTCCAGCCACTGCTGCTCGTTGCCGTTGATTACCGCCGTAACGCCCg                                            >  1:384029/1‑90 (MQ=255)
                                               gcagaaATTTCTTTTCATACTGACAGGCGGGGGATCCCAGCCCCAGCCACTGCTGCTCGTTGCCGTTGATTACCGCCGTAACGCCCGCAg                                         <  2:69181/90‑1 (MQ=255)
                                                 agaaaaTTCTTTTCATACTGACAGGCTGGGGATTCCAGCCCCAGCCACTGCTGCTCGTTGCCGTTGATTACCGCCGTAACGCCCGCAGGt                                       <  2:384029/87‑1 (MQ=255)
                                                        tctcTTCATACTGACGGGCGGCGGATTCCAGCCCCAGCCACTGCTGCTCGTTGCCGTTGATTACCGCCGTAACGCCCGCAGGTAATTCaa                                <  2:466103/90‑1 (MQ=255)
                                                            ttCATACTGACAGGCTGCGGATTCCAGCTCCAGCCACTGCTGCTCGTTGCCGTTGATTACCGCCGTAACGCCCGCAGGTAATTCAACAGc                            >  2:333540/1‑90 (MQ=255)
                                                                  cTGACGGGCTGCGGATCCCAGCCCCAGCCCCTGCTGCTCGTTGCCGTTGATTACCGCCGTAACGCCCGCAGGTAATTCAACAGCATGTTc                      <  2:66109/90‑1 (MQ=255)
                                                                  cTGACGGGCGGGGGTTTCCAGCCCCAGCCACTGCTGCTCGTTGCCGTTGTTTACCGCCGTAACGCCCGCAGGTAATTCAACAGCATGTTc                      <  2:101893/90‑1 (MQ=255)
                                                                      cAGGCGGGGGTTCCCAGCCCCAGCCACTGCTGCTCGTTGCCGTTGATTACCGCCGTAACGCCGCAGGTAATTCAACAGCATGTTCGGAAt                 <  2:286171/90‑1 (MQ=255)
                                                                                      gctccagcCACTGCTGCTCGTTGCCGTTGATTACCGCCGTAACGCCCGCAGGTAATTCAACAGCATGTTCGGAATACAAAAACACCGGTa  <  1:333540/90‑1 (MQ=255)
                                                                                        |                                                                                       
GCTGTTGCCCATCTCAACGTACTGCGTCAGCGTGTCATAAAACGGTGGCAGCAAATTCTCTTCATACTGACAGGCTGCGGATTCCAGCTCCAGCCACTGCTGCTCGTTGCCGTTGATTACCGCCGTAACGCCCGCAGGTAATTCAACAGCATGTTCGGAATACAAAAACACCGGTA  >  NZ_CP009273/3102167‑3102342

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