Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,058,9580GA54.5% 6.2 / 25.6 22L57L (CTG→TTG) yihOMFS transporter
Reads supporting (aligned to +/- strand):  ref base G (7/3);  new base A (0/12);  total (7/15)
Fisher's exact test for biased strand distribution p-value = 7.04e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGT  >  NZ_CP009273/4058870‑4059027
                                                                                        |                                                                     
aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa                                                                      <  1:599589/90‑1 (MQ=255)
             cccaaaTTTTTTCCCGGGAGTCGAAAAAAAACCCGGTGAGCATATCGGAAAACGCGGAAAAAACCTCCCCCCCAAAAAAAATAACCCCGc                                                         <  2:146719/86‑1 (MQ=255)
                      ttttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAg                                                >  1:456580/1‑90 (MQ=255)
                       tttCGGCGAGTCGAAAAAAAACCGGGTGACCATACGGGAAAACGCGGAAAAAACCTCCCCCACAAAAAAGATAACCCCGCCATAGTAAGc                                               <  2:910227/90‑1 (MQ=255)
                       tttCCGCGAGCCAAAAAAAAATCGGGTGACCATACGGGAAAACGCGAAAAAAACCTCCCCCACCAAAAAAATACTCCCGCCATAGTAAGc                                               <  2:703764/90‑1 (MQ=255)
                           cgcgAGTCGAGAAAAAATCGGGGGAGCATCTGGGAAAACGGGGTAAAAACCTCCGCCACCAGAAAAAAAATCCCGCCATAGTAAGCAGGc                                           <  1:896963/90‑1 (MQ=255)
                           cgcgAGTCGAGAAAAAACCCGGTGAGCATATCGGTAAACGCGGAAAAAAACTTCCCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGc                                           <  1:295553/90‑1 (MQ=255)
                            gcgAGTGGAAAAAAAACCCGGGGAGCATATCGGAAAACGCGGAAAAAAACTCCGCCACAAAAAAGATAACCCCGCCATAGTAAGCAGGCa                                          <  1:619548/90‑1 (MQ=255)
                            gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                          >  1:520729/1‑90 (MQ=255)
                            gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                          >  1:530998/1‑90 (MQ=255)
                            gcgAGTAGAGAAAAAACCGGGGGAGCATATCGGAAAACGGGAAAAAAACCTTCCCCACCAAAAAAAAAACCCCGCCATAGTAAGCAGGCa                                          <  2:987909/90‑1 (MQ=255)
                                         aaaTCGGGTGAGAATATCGGAAAACGCGAAAAAAACCTCCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                             <  2:530998/90‑1 (MQ=255)
                                         aaaTCCGGTGAGCATATCGGTAAACGCGGAAAAAAACTTCCCCCCCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                             <  1:328850/90‑1 (MQ=255)
                                          aaTCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg                            >  1:299593/1‑90 (MQ=255)
                                               ggTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata                       >  2:614409/1‑90 (MQ=255)
                                               ggTGAGCATATCGGAAACCGGGGAAAAAAACTTCCCCCCAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata                       <  2:388625/90‑1 (MQ=255)
                                                       tagcGGAAAACGGGATAAAAACCTCCGCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga               <  1:497686/87‑1 (MQ=255)
                                                            gTAACCGCGAAAAAAACCTTCCCCACCAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga          <  1:471761/90‑1 (MQ=255)
                                                              aaaCGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATa        >  2:593118/1‑90 (MQ=255)
                                                               aaCGCGAAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc       <  2:913185/90‑1 (MQ=255)
                                                                  gcgGAAAAAAACTTCGCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga    <  1:342689/90‑1 (MQ=255)
                                                                    ggTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt  >  1:913185/1‑90 (MQ=255)
                                                                                        |                                                                     
AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGT  >  NZ_CP009273/4058870‑4059027

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