Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,659,0620CA61.5% 0.5 / 17.1 13A37S (GCG→TCG) gadXacid resistance transcriptional activator GadX
Reads supporting (aligned to +/- strand):  ref base C (5/0);  new base A (0/8);  total (5/8)
Fisher's exact test for biased strand distribution p-value = 7.77e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAC  >  NZ_CP009273/3658978‑3659148
                                                                                     |                                                                                       
tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                                    >  2:13159/1‑90 (MQ=255)
     agaAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAgg                                                                               >  1:122043/1‑90 (MQ=255)
           aatAAAAGCGTGTCCCTTGACACAAAAACAAAATTTTCAACACCCTTATCTACTCAAATTTGGCTTGCATCCGAAAAAACCAGGTCACCg                                                                        <  2:49169/90‑1 (MQ=255)
                      gtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATa                                                             >  1:128847/1‑90 (MQ=255)
                                      taaaaaTTTTCCAACCCACTTATCTACTCGATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCAtt                                              <  1:224346/88‑1 (MQ=255)
                                               tCAACACATTTATCTACTCAATTTTGGCTTGCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                    <  2:37433/90‑1 (MQ=255)
                                                ccAACACTTTTATCTCCTCGATTTTGGCTTCCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                    <  2:122043/89‑1 (MQ=255)
                                                     ccacTTATCTACTCAATTTTGGTTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat                               <  2:138337/89‑1 (MQ=255)
                                                       atttATCTCCTCGATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatat                             <  1:13159/88‑1 (MQ=255)
                                                                cTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTc                    >  1:343456/1‑90 (MQ=255)
                                                                        tGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc            <  1:45249/90‑1 (MQ=255)
                                                                        tGGCTTCCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc            <  1:147026/90‑1 (MQ=255)
                                                                                  ccGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc  >  2:328173/1‑90 (MQ=255)
                                                                                   cccAAAAACCCGGTTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc  <  2:239832/88‑1 (MQ=255)
                                                                                     |                                                                                       
TCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAC  >  NZ_CP009273/3658978‑3659148

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