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

CAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATA  >  NZ_CP009273/585301‑585462
                                                                                     |                                                                            
cAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCggggggggg                                                                          >  2:578699/1‑89 (MQ=255)
    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)
                                                                                     |                                                                            
CAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATA  >  NZ_CP009273/585301‑585462

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: