Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273612,0120TG53.3% 10.7 / 21.4 15G800G (GGT→GGGentFenterobactin non‑ribosomal peptide synthetase EntF
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base G (8/0);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04
Rejected as consensus: Frequency below/above cutoff threshold.

TGTCAGCTGGTATCCGGCTTTTGGCGAGGAACTGGCACAGGTGCGCGGCAGCAGTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCT  >  NZ_CP009273/611874‑612142
                                                                                                                                          |                                                                                                                                  
tGTCAGCTGGTATCCGGCTTTTGGCGAGGAACTGGCACAGGTGCGCGGCAGCAGTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGt                                                                                                                                    <  2:222531/139‑1 (MQ=255)
  tCAGCTGGTATCCGGCTTTTGGCGAGGAACTGGCACAGGTGCGCGGCAGCAGTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTg                                                                                                                                   <  1:347664/138‑1 (MQ=255)
                                             cgGCAGCAGTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGGGGTATTTATCACACCGGCTTTATACCGGCGCAGGGATCTCCCggg                                                                                       >  2:131236/1‑138 (MQ=255)
                                                        tgCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGGGGTCTTTTTCTCACTGGCATTTAAATGGCGCCGGGGCATTTCGGCCGcccccgcct                                                                            >  2:57390/1‑134 (MQ=255)
                                                        tgCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGGGGGTTCTTTTTCTCCTGCCTTTAACTTGCGCCGGGGGATCTCGGGCGCCCCgcact                                                                            >  1:16636/1‑135 (MQ=255)
                                                             aTTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTg                                                                                                                                   <  1:88625/79‑1 (MQ=255)
                                                             aTTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTg                                                                                                                                   >  2:88625/1‑79 (MQ=255)
                                                                                    gggTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGGGATTTCTTTTTCACTGGCATTCAACTGGGGGAGGGGTATCTCGGGCGCCCCCCTCTGACCCCCAGCCGCCTTATTTCCGATCCt                                                >  1:136815/1‑139 (MQ=255)
                                                                                    gggTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGGGATCTCTTTCTCCCTGGCCTTCCACTGGCGGCGGGGGCTCTCGGGCGCCCCGGTCTGACCGCCCGCCGCCTTTTTTCCCATCCt                                                >  1:419672/1‑139 (MQ=255)
                                                                                     ggTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCtt                                               <  1:353812/139‑1 (MQ=255)
                                                                                      gTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGGGATCTCTTTCTCTCTCGCATTCAACTGGCGCAGGGGTGTCTCGGACGCCCCCCTCTTTCCGCCCGCCGCCTTTTTTCCCATCCttt                                              >  1:492788/1‑139 (MQ=255)
                                                                                               ttCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAg                                      <  2:16636/138‑1 (MQ=255)
                                                                                                          atgatgCATCCGGTGCCGCCGGGTGTGGCGGGGGATCTCTCTCTCACTGGCGTTCAACTGGCGGCGGGGGGTCTCGGGCGCCCCCGCCTGTCCGCCCGCCCCCTTTTTGCCGCGCCTTTTTCCCCCGGGGGACGGgtgt                          >  2:46583/1‑139 (MQ=255)
                                                                                                           tgatgCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATgta                         <  1:57390/139‑1 (MQ=255)
                                                                                                                     ggTGCCGCCGGGTGTGGCGGGGGATATCTATCTCACCCGCATTCAACCGGCGCCGGGGTATCTCGGGCGCCCCGCTCTGACCGCCAGCCGCCTTTTTTCCGATTCTTTTTTCCCAGGGGGACGGGCgtaccctacccgg               >  2:228920/1‑136 (MQ=255)
                                                                                                                                   tgGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCgct  <  1:52492/138‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                  
TGTCAGCTGGTATCCGGCTTTTGGCGAGGAACTGGCACAGGTGCGCGGCAGCAGTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCT  >  NZ_CP009273/611874‑612142

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