Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273585,3900TG57.1% 11.7 / 25.8 21T341P (ACC→CCC) nfrAbacteriophage adsorption protein NfrA
Reads supporting (aligned to +/- strand):  ref base T (0/9);  new base G (12/0);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 3.40e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCA  >  NZ_CP009273/585305‑585478
                                                                                     |                                                                                        
cGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGgggggggggggg                                                                                      >  2:42740/1‑88 (MQ=255)
 gTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGTGGGGGCgg                                                                                     >  1:443356/1‑89 (MQ=255)
    ccAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCa                                                                                  <  1:313161/90‑1 (MQ=255)
              agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTgcgc                                                                        <  1:410330/90‑1 (MQ=255)
                cTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTg                                                                      <  2:698787/90‑1 (MQ=255)
                   cAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGcgt                                                                   <  2:144291/90‑1 (MQ=255)
                          gTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTg                                                            <  1:723946/90‑1 (MQ=255)
                             aCGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGAcc                                                         <  1:42740/90‑1 (MQ=255)
                             aCGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGAcc                                                         <  1:220782/90‑1 (MQ=255)
                               gggTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGCGGGTTTGCGCCGCGTCGGAGTGagcct                                                       >  1:46413/1‑86 (MQ=255)
                                   cGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGGGCCGGTTTTTCGCTGCGGCGTACCGGCCCTCttt                                                   >  2:642512/1‑90 (MQ=255)
                                    gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTTGCTGGGGCGTGCTTACCTTCtttt                                                  >  1:537626/1‑90 (MQ=255)
                                                        ctcctcAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTTGCTGCGTCGTTCTTACCCTCTTTTTACACGACCGGGAGGGTggc                              >  1:554502/1‑90 (MQ=255)
                                                          cctcAAGCATTTCATTGGCGGGGGGGGGGGGGGGGGGGTTTTTCGCTGCGGCGTACTGACCCTCTTTTTAAAACACCCGGAGGGgcgcgc                            >  2:218013/1‑90 (MQ=255)
                                                          cctcAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTCGCTGCGTCGTACTGACCTTTTTTTTAAAACACCGGGAGCGgcgcgc                            >  2:395979/1‑90 (MQ=255)
                                                           ctcAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGTTTTTTCCCTGCGTCGTACTTACCCTCTTTTTTAAACACCCGGAGCGTCGCGcc                           >  1:715335/1‑90 (MQ=255)
                                                           ctcAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTCTCCGCGGCGGTCTTACCTTCTTTTTAAAACACCGGGGGGGGCGCGcc                           >  1:74061/1‑90 (MQ=255)
                                                           ctcAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGcc                           <  1:725205/90‑1 (MQ=255)
                                                                    ttCATTGGCGGGGGGGGGGGGGGGCAGTTTTTTTGCTGCGGCGTTCTGACCTTCTTTTTATAACACCGGGAGCGGCGCGCCCACAACATa                  >  1:32874/1‑90 (MQ=255)
                                                                    ttCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTCGCTGCGGCGTACTGACCTTCTTTTTAAAACACCGGGAGCGGCGCGCCCacaacaaa                  >  2:478109/1‑88 (MQ=255)
                                                                    ttCATTGGCGGGGGGGGGGGCGAGGAGGTTTTTCGCTGCGTCGGGCTGACCTTCTTTTTACAACACCGGGGGCGGCGCGCCCacaacaaa                  >  2:232116/1‑88 (MQ=255)
                                                                                    gTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCa  <  1:427939/90‑1 (MQ=255)
                                                                                     |                                                                                        
CGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCA  >  NZ_CP009273/585305‑585478

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