Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,611,7230GA41.7% 22.9 / 11.2 24K196K (AAG→AAAnikEnickel import ATP‑binding protein NikE
Reads supporting (aligned to +/- strand):  ref base G (3/11);  new base A (10/0);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTAAGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGTGGTGGGAGAG  >  NZ_CP009273/3611588‑3611849
                                                                                                                                       |                                                                                                                              
gTTAAGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGGGTCATTTGCCCGCTGAAAAAACTa                                                                                                                             >  1:66477/1‑139 (MQ=255)
    aGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTacaac                                                                                                                         <  1:388358/139‑1 (MQ=255)
    aGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTacaac                                                                                                                         <  2:123347/139‑1 (MQ=255)
         cggcCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAATCTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGGGTTATGCCTCCTGTAAAAACCacaacaacc                                                                                                                     >  2:402960/1‑137 (MQ=255)
                                      gcgcgcTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCCCCGCCTGCCTGTTCATCACCCa                                                                                       <  1:140958/139‑1 (MQ=255)
                                      gcgcgcTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCa                                                                                       <  1:56693/139‑1 (MQ=255)
                                                 gTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGcc                                                                            <  1:16095/139‑1 (MQ=255)
                                                            actactGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCCTTCGCCCGCTTAAAAAACTACAACAACAACTTGGCACCGCCCGCCTGTTCATCACCCCCGACTTAAGCCCGGTGGAAAGc                                                                 >  2:244790/1‑139 (MQ=255)
                                                                    ttCTGGATGACGCTGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCCCGACTTACGCCTGGTAGAACGCTTTTGCCa                                                         <  2:163249/139‑1 (MQ=255)
                                                                    ttCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCa                                                         <  1:236505/139‑1 (MQ=255)
                                                                     tCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAACTACAACAACAGTTTGGGACCGCCCGCCCGTTCCTTACCCCCCACTTATGCCTGGTAGAAAGCTTTTGcccg                                                        >  2:387662/1‑137 (MQ=255)
                                                                     tCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGGCATTTGCCTCCTGAAAAAAATACAACAACAGATTGGCGCCGCCCGCCTGTTTCTCACCCCCCACCTACCCCCGGGAGGACGCCTTTGcccg                                                        >  2:654568/1‑137 (MQ=255)
                                                                            gagacagTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa                                                 <  1:579851/133‑1 (MQ=255)
                                                                             aaGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa                                                 <  2:43650/138‑1 (MQ=255)
                                                                             aaGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGGGTCATTCGCCTGCTGAAAAAACTACAACAAAAATTTTGCACCGCCCGCCTGTTCATCACCCACCACCTACGCCCGGTAGAAAGCTTTTGCCAGCgggggat                                                >  1:60893/1‑135 (MQ=255)
                                                                                               aTCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGAc                               <  1:654568/138‑1 (MQ=255)
                                                                                               aTCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAAATAAAACAACAGTTTGGGACCGCCCGCCCGTTCATTACCCACGACCTACGCCCGGTAGAACGCTTTTTCCCGCGGGGAAATGTTTTGGGCAACGGAc                               >  2:560259/1‑138 (MQ=255)
                                                                                                   cGTGTTACAGGCGGGGGTCATTCGCCTGCTGAAAAAACTACAAAAAAAGTTTGGCACCGCCTGCCTGTTCATCACCCCCGCCTTACGCCCGGGAGAACGCTTTTGCCAGCGGGGAATGGGTATGGACAACAGACaaaac                          >  1:128592/1‑137 (MQ=255)
                                                                                                     tgtTACAGGCGGGTGTCATTTTCCTGCTGAAAAAACTACAACAACAGTTTGGCACCGCCCGCCTGTTCATCACCCACGACTTACGCCCGGTAGAACGCTTTTTCCAGCGGGGAATGGTTATGGACAAAGGACAAAtcgt                        >  1:140024/1‑139 (MQ=255)
                                                                                                     tgtTACAGGCGGGTGTCATTCGCCCGCTGAAAAAGCTACAACCAAAGATTGGCACCGCCCGCCTGTTCATCAACCCCGACTTACGCCTGGTAGAACGCTTTTTCCAGCGGGGAATGGTTATGGACAACGGACAAAacgg                        >  2:638169/1‑138 (MQ=255)
                                                                                                            ggCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAAcc                 <  2:60893/139‑1 (MQ=255)
                                                                                                                ggTGTTATTCTCCCGCTGGAAAAACAACAAAAAAAATTTGGGACCGCCCGCCTGTTCATCACCCACGACTTTTGCCCGGTAGAAAGCTTTTTCCACCGGGTAATGGGTTTGGACAACGGACAAAACGGCGAAACCCAgg             >  2:496637/1‑139 (MQ=255)
                                                                                                                         cGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGtggtgg       >  1:633832/1‑136 (MQ=255)
                                                                                                                            ctgctgAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGTGGTGGgagag  >  1:379834/1‑138 (MQ=255)
                                                                                                                                       |                                                                                                                              
GTTAAGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGTGGTGGGAGAG  >  NZ_CP009273/3611588‑3611849

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