Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,368,6190AC57.1% 12.7 / 26.5 21V48G (GTA→GGA) ampCBlaEC family class C beta‑lactamase
Reads supporting (aligned to +/- strand):  ref base A (9/0);  new base C (0/12);  total (9/12)
Fisher's exact test for biased strand distribution p-value = 3.40e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-05
Rejected as consensus: Frequency below/above cutoff threshold.

AAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCA  >  NZ_CP009273/4368534‑4368708
                                                                                        |                                                                                         
aaaCAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgc                                                                                       >  1:175778/1‑90 (MQ=255)
aaaCAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgc                                                                                       >  2:232686/1‑90 (MQ=255)
          tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATa                                                                             >  2:675304/1‑90 (MQ=255)
               tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCgg                                                                        >  1:352683/1‑90 (MQ=255)
                gtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCggg                                                                       >  2:134536/1‑90 (MQ=255)
                   aCGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATc                                                                    >  1:460193/1‑90 (MQ=255)
                     gggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCtt                                                                  >  2:135236/1‑90 (MQ=255)
                           ttttttGGCGATCCCCCCATACCCCCGGGAAAAGTAATAAGGTTTCCCCTAATAAATTCCCGCCACCGCCATACCCGGGATCTTttgttg                                                            <  2:12952/89‑1 (MQ=255)
                            tttttGGCGGTCCCCCCATCCCCCGGGGAAAAATAATATGTTTTCCCCTAAAAAATTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGc                                                           <  1:351397/90‑1 (MQ=255)
                            ttttgGGCGATCCCCCCATACCCCAAGAAAAAAAAAAATGTTTTCCCTTAAAAATTTCCCCCCCCCGCCATACCCGGGAACTTTTGTTGc                                                           <  1:445406/90‑1 (MQ=255)
                                     atccccGCATCGCCCAAGGAAAAGAAATATTTTTTCCCTTAAAAAATCCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc                                                  <  1:350909/86‑1 (MQ=255)
                                     aTGCCCGCATACCCCCAGGAAAAAAAATAAGTTTTCCCCTAAAAAATCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc                                                  <  1:102746/90‑1 (MQ=255)
                                               ccccccAGGAAAAGAAATAGGGTTTCCCCTAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                        <  2:175778/88‑1 (MQ=255)
                                                gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgt                                       >  1:449849/1‑90 (MQ=255)
                                                  cccGGGAAAAGTAATAAGTTTTCCCCTAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc                                     <  2:90582/90‑1 (MQ=255)
                                                   ccAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                    >  2:704705/1‑90 (MQ=255)
                                                   ccAGGAAAAATAAAAAAGTTTTCCCTTAAAAAATTCCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc                                     <  2:319498/90‑1 (MQ=255)
                                                    cgggTAAAGAAATAATTTTTCCCTTGAAAATTTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGa                                   <  2:352683/88‑1 (MQ=255)
                                                     gagTAAAATAATTATGTTTTCCCTTAAAAATTCCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGa                                   <  1:696384/88‑1 (MQ=255)
                                                         aaaGAAATAAGTTTTCCCCTAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                              <  2:283099/90‑1 (MQ=255)
                                                                      ttACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                 >  1:580158/1‑90 (MQ=255)
                                                                       tccccTAAAAAATCCCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                <  2:101464/88‑1 (MQ=255)
                                                                         cccTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt              >  2:517441/1‑90 (MQ=255)
                                                                                        accgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca  >  1:3032/1‑90 (MQ=255)
                                                                                        |                                                                                         
AAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCA  >  NZ_CP009273/4368534‑4368708

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