Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273585,3900TG60.0% ‑6.5 / 14.4 15T341P (ACC→CCC) nfrAbacteriophage adsorption protein NfrA
Reads supporting (aligned to +/- strand):  ref base T (0/6);  new base G (9/0);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 2.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAA  >  NZ_CP009273/585321‑585479
                                                                     |                                                                                         
cTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTg                                                                       <  2:285736/90‑1 (MQ=255)
  tCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTCGCTGcg                                                                     >  2:18971/1‑90 (MQ=255)
         tGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTCGCTGCGGCGTACt                                                              >  1:526069/1‑90 (MQ=255)
         tGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGGGGGGAGGTTTTTTGCTTCGTCGGACt                                                              >  1:974000/1‑90 (MQ=255)
             aCGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGAcc                                                          <  2:638341/90‑1 (MQ=255)
               gggTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGTTTTTTCGCTGCGTCGGACTGACctt                                                        >  1:243863/1‑90 (MQ=255)
               gggTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTTGCTGCGGCGTACTGACctt                                                        >  1:633799/1‑90 (MQ=255)
                  tCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCGGGTTTTTCGCTGCGTCGTACTGACcttctt                                                     >  2:525087/1‑90 (MQ=255)
                   cGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAGGTTTTTCGGCGCGGCGTACTTGCCTTCttt                                                    >  1:141841/1‑90 (MQ=255)
                    gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTGGTGGCGGCGGCCTGACCTtttttt                                                   >  2:912462/1‑90 (MQ=255)
                                       gCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTcg                                <  1:525087/90‑1 (MQ=255)
                                          cctcAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAGTTTTTTCGCTGCGTCGTACTGACCTTCTTTTTACAGCACCGGGAGCGgcgcgc                             >  1:279381/1‑90 (MQ=255)
                                           ctcAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGcc                            <  2:500335/90‑1 (MQ=255)
                                            tcAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCca                           <  2:924805/90‑1 (MQ=255)
                                                                  gggTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCgg     <  1:912462/90‑1 (MQ=255)
                                                                     tGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCaa  <  1:954441/90‑1 (MQ=255)
                                                                     |                                                                                         
CTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAA  >  NZ_CP009273/585321‑585479

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