Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,408,0140AT54.5% 6.1 / 13.2 11N57K (AAT→AAAralRtype I toxin‑antitoxin system endodeoxyribonuclease toxin RalR
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base T (6/0);  total (7/4)
Fisher's exact test for biased strand distribution p-value = 1.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGC  >  NZ_CP009273/1407931‑1408103
                                                                                   |                                                                                         
ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAAAACCCCCCTtttttttttta                                                                                     >  2:192256/1‑87 (MQ=255)
      atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGtt                                                                               <  1:85787/90‑1 (MQ=255)
                          gTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTAATAAAACCCCCCTTATTTTTTTCCCTGGGTTGTTTTTTTCCTCtttttag                                                           >  1:634693/1‑86 (MQ=255)
                          gTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAAAACCTCCCTTTTTTTTTTCCGTGGGTTGTTTTTTTCCTCttttgtg                                                           >  1:225485/1‑86 (MQ=255)
                            gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAAAACCTCCCTTTTTTTTTTCCGTGGTTTGTTTTTTTCATTTTTTGGGTg                                                         >  1:371692/1‑90 (MQ=255)
                            gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAAAACCTCCCTTTTTTTTTTCCGTGGGTTGTTTTTTTCATCTTTTGAGTg                                                         >  1:222929/1‑90 (MQ=255)
                                    ccttTTCAATAGTGGCGGTAATTTTTTACATTAAAACCCCCATTTTTTTTTTCAGTGGTTTGTTTTTTTCCTCTTTTTAGTGGTtttttt                                                 >  1:5539/1‑90 (MQ=255)
                                      ttttCAATAGTGGCGGTAATTTTTTAAATGAAAACCCCCCTTATTTTTTTCATTGGGTTGTTTTTTTCCTCTTTTGAGTGCTTCtttttt                                               >  2:137890/1‑89 (MQ=255)
                                                      tAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCg                               <  1:514558/90‑1 (MQ=255)
                                                            tataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCgg                         <  1:23316/90‑1 (MQ=255)
                                                                    aaCACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCa                 <  2:29636/90‑1 (MQ=255)
                                                                          tCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGtt           >  2:298561/1‑90 (MQ=255)
                                                                             attattattTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTaca        >  1:682197/1‑90 (MQ=255)
                                                                                  tattTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCg   >  2:153224/1‑90 (MQ=255)
                                                                                  tattTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCg   >  2:8533/1‑90 (MQ=255)
                                                                                   attTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGc  <  1:298561/90‑1 (MQ=255)
                                                                                   |                                                                                         
TTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGC  >  NZ_CP009273/1407931‑1408103

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