Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,356,9000TC53.8% 5.3 / 13.8 13G71G (GGA→GGGyfaYnicotinamide mononucleotide deamidase‑related protein YfaY
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.

TCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCAC  >  NZ_CP009273/2356811‑2356986
                                                                                         |                                                                                      
tCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTggtcggt                                                                                        >  1:288273/1‑90 (MQ=255)
                gggCTTCATGCACCCCCGGCCCTCCCCTTTTTGCTGCCGCGGGGGCGAGGGCGCTAAAATCACCGCGGGTCGGCCCCAGCCCGCCGTTaa                                                                        <  2:130129/89‑1 (MQ=255)
                 ggCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAAc                                                                       >  2:212174/1‑90 (MQ=255)
                   ctccATGCACCCCCAGGCCTTCCCTTTTTGCTGTCGCGGCGGGGAGTGCGCTTAAATCATCGCGGGTGGGCCCCAGCCCGCCGTTAACGa                                                                     <  1:154353/87‑1 (MQ=255)
                          cagcaCCAGCCCTTCCCTTTTTGCTGCCGCAGGGGCGAGTGCGTTTAAATCATCGCGGGTGGGCCCCAGCCCGCCGTTAACGATCAGcac                                                              <  1:212174/90‑1 (MQ=255)
                            gcaCCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACAt                                                            >  1:275137/1‑90 (MQ=255)
                                cggcccTTCCCTTTTTGCTGTCGCGGGGGGGGGTGCGTTTAAATCACCGCGGGTCGGCCCCACCCCGCCGTTAACGATCAGCACATCGGc                                                        <  1:375984/86‑1 (MQ=255)
                                 aGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCg                                                       >  1:256015/1‑90 (MQ=255)
                                         cctttttGCTGCCGCGGGGGCGGGTGCGCTTAAACCATCGCTGGTCGGCCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACt                                               <  1:208968/87‑1 (MQ=255)
                                             tttGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGt                                           >  1:176567/1‑90 (MQ=255)
                                                   gTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCga                                     >  1:291234/1‑90 (MQ=255)
                                                         gcggcgAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATg                               >  2:82014/1‑90 (MQ=255)
                                                                           tcatcGCTGGTCGGCCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAgg             <  1:300257/90‑1 (MQ=255)
                                                                            ccgcgcgGGTCCGCCCCCCCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGt            <  2:58168/87‑1 (MQ=255)
                                                                                      cggccCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCAc  <  1:97803/90‑1 (MQ=255)
                                                                                         |                                                                                      
TCCATCTCTTTGAGCCAGGCTTCATGCAGCACCAGGCCTTCACCTTTTGCTGTCGCAGCGGCGAGTGCGCTTAAATCATCGCTGGTCGGTCCCAGCCCGCCGTTAACGATCAGCACATCGGCGTGCTGACTACGTTCGCGAAGAATGGTGACTAAGTCATCAAGGTTATCCCCCAC  >  NZ_CP009273/2356811‑2356986

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