Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,659,0580TC61.1% ‑0.5 / 10.1 18D38G (GAT→GGT) gadXacid resistance transcriptional activator GadX
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base C (0/11);  total (4/14)
Fisher's exact test for biased strand distribution p-value = 1.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGC  >  NZ_CP009273/3658977‑3659140
                                                                                 |                                                                                  
aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaa                                                                            >  1:48408/1‑90 (MQ=255)
           aaatGAAAGCGTGCCCCTTGACACAAATAAAAAATTTTCAACACACTTATCTACTCAAATTTGGCTTGCACCCGAAAAAACCAGGTCAcc                                                                 <  1:96927/89‑1 (MQ=255)
                      ttccccTTAACAAAAATAAAAATTTTCCAACCCCCTTATCTTCTAAATTTTGGCTTCCACCCGAAAAACCCAGGCCACCGCCATTAAAAt                                                      <  2:352719/87‑1 (MQ=255)
                         cccTTGACACAAATAAAAAATTTCCAACACACTTATCTCCTCGATTTTGGCTTGCATCCGCAAAACCCAGGTCACCGCCATTAAAATAGc                                                   <  1:18578/90‑1 (MQ=255)
                         cccTTAACACAAATAAAAAATTTCCAACACATTTATCTCCTCGATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGc                                                   <  2:39772/90‑1 (MQ=255)
                                cacGAATACAAAATTTTCAACCCATTTATCTCCTCAATTTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTc                                            <  1:63664/90‑1 (MQ=255)
                                   gAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc                                         >  1:10101/1‑90 (MQ=255)
                                     aTACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGa                                                 >  2:11100/1‑80 (MQ=255)
                                     aTACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGa                                                 <  1:11100/80‑1 (MQ=255)
                                                taaaCACTTTTATTTACTAAATTTTGGCTTCCACCCGAAAAACCCGGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                            <  2:48408/88‑1 (MQ=255)
                                                tCAACACCTTTATCTACTCAAATTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                            <  1:131979/90‑1 (MQ=255)
                                                      cccTTATCTACTAGATTTTGTTTTCCACCCAAAAAAACCGGCTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAata                      <  2:233428/88‑1 (MQ=255)
                                                                  cGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCtt          >  1:244477/1‑90 (MQ=255)
                                                                        tGGTTTGCACCCAAAAAAACCGGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc    <  2:214675/90‑1 (MQ=255)
                                                                        tGGCTTCCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc    <  1:257262/90‑1 (MQ=255)
                                                                        tGGCTTCCACCCAAAAAAACCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc    <  2:44965/90‑1 (MQ=255)
                                                                        tGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc    <  1:25407/90‑1 (MQ=255)
                                                                          gCTTGCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGc  <  1:434377/90‑1 (MQ=255)
                                                                                 |                                                                                  
ATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGC  >  NZ_CP009273/3658977‑3659140

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