Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273585,3860AG60.9% 7.1 / 24.3 23L342P (CTC→CCC) nfrAbacteriophage adsorption protein NfrA
Reads supporting (aligned to +/- strand):  ref base A (0/9);  new base G (14/0);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 1.22e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAAC  >  NZ_CP009273/585300‑585459
                                                                                      |                                                                         
gCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCgggggggg                                                                        >  1:437175/1‑90 (MQ=255)
   aTCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGgggggggggg                                                                     >  1:136725/1‑90 (MQ=255)
                   agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTgcgc                                                     <  2:486956/90‑1 (MQ=255)
                   agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTgcgc                                                     <  1:544880/90‑1 (MQ=255)
                   agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTgcgc                                                     <  2:634766/90‑1 (MQ=255)
                       tCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGGAGGTTTTGTGCcgcg                                                 >  2:769742/1‑90 (MQ=255)
                              tGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTCGCTGCGTCGGACt                                          >  1:757363/1‑90 (MQ=255)
                              tGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGGTTTTCTCTGcgtcggcct                                          >  2:616209/1‑86 (MQ=255)
                                  aCGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGAcc                                      <  1:51865/90‑1 (MQ=255)
                                  aCGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGAcc                                      <  2:757063/90‑1 (MQ=255)
                                        cGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTCGCCGCGTCGTTCTTACCTTCttt                                >  1:400875/1‑90 (MQ=255)
                                         gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTTGCTGCGGCGGACTTACCTTCtttt                               >  2:175948/1‑90 (MQ=255)
                                             cGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACa                           <  2:268371/90‑1 (MQ=255)
                                                tGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCa                        <  2:136725/90‑1 (MQ=255)
                                                   cAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGTGGTTTTTGGGCTGCCGCGGATTTGCCTTTTTTTTATAGCACCg                     >  1:822848/1‑90 (MQ=255)
                                                               cctcAAGCATTTCATTGGCGGGGGGGGTGGGGGGCAGATTTTTCTCTGCGGCGGTGGTACCTTTTTTTTAAAGCACCCGGAGGGgcgcgc         >  1:617553/1‑90 (MQ=255)
                                                               cctcAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGTTTTTTTGCTGCGTCGTACTGACCCTTTTTTTAAAACACCCGGAGCGGCccgc         >  1:632116/1‑90 (MQ=255)
                                                               cctcAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTCGCGGCGTGGTACTGACCTTCTTTTTAAAACACCGGGAGCGgcgcgc         >  1:91241/1‑90 (MQ=255)
                                                               cctcAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGTTTTTTCGCTGCGGCGTACTGACCTTCTTTTTAAAACACCCGTAGCGgcgcgc         >  1:797984/1‑90 (MQ=255)
                                                               cctcAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTTGCTGCGTCGGACTGACCTTTTTTTTACAACACCGGGAGCGgcgcgc         >  1:452672/1‑90 (MQ=255)
                                                               cctcAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAGGTTTTTCGCTGCGTCGTACTGACCTTTTTTTTAAAGCACCGGGAGCGgcgcgc         >  1:221219/1‑90 (MQ=255)
                                                                ctcAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGcc        <  2:632116/90‑1 (MQ=255)
                                                                  cAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTCGCTGCGCCGGACTGACCTTTTTTTTAAAACACCGGGAGGGGCGGGCcaa      >  2:59308/1‑90 (MQ=255)
                                                                    aGCATTTCATTGGCGGGGGGGGGGGCGGGCAGTTTTTTCGCTGCGTCGTACTTACCTTCTTTTTACAGCACCGGGAGCGGCGCGCccaca    >  1:284168/1‑90 (MQ=255)
                                                                      cATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCaacaac  <  1:59308/90‑1 (MQ=255)
                                                                                      |                                                                         
GCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAAC  >  NZ_CP009273/585300‑585459

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