Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,012,9220TG52.6% 18.6 / 23.9 19S135A (TCG→GCG) ssnAputative aminohydrolase SsnA
Reads supporting (aligned to +/- strand):  ref base T (1/8);  new base G (10/0);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.19e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAA  >  NZ_CP009273/3012836‑3013010
                                                                                      |                                                                                         
gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGttcg                                                                                        >  1:191060/1‑87 (MQ=255)
gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGgcgc                                                                                        >  1:38861/1‑90 (MQ=255)
gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCgg                                                                                        >  2:567757/1‑89 (MQ=255)
 cGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGct                                                                                       <  2:190754/90‑1 (MQ=255)
  ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCggcc                                                                                      >  2:505934/1‑87 (MQ=255)
    aCTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCtctc                                                                                    <  1:235788/90‑1 (MQ=255)
                                 gCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                                                                                           >  1:125639/1‑54 (MQ=255)
                                 gCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                                                                                           >  2:313588/1‑54 (MQ=255)
                                 gCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                                                                                           <  2:125639/54‑1 (MQ=255)
                                 gCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                                                                                           <  1:313588/54‑1 (MQ=255)
                                  cGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCCTTTCGCGACGCATTTTTTAAAATTgg                                                     >  1:261883/1‑90 (MQ=255)
                                   ggATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGGGGGCTCCCCACTTCCCGACGCCCTTTTTAAAATTTgc                                                     >  2:445675/1‑89 (MQ=255)
                                            tCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCCCCCCCCTTCGCGGCGCCTTTTTTTAAATTTGGCTTcgcgcg                                           >  2:205380/1‑90 (MQ=255)
                                                 tATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCGTTCCCCAATGCGCGACGCATTTTTTAAAAATGGGCTGCGCCCGCTGag                                      >  1:388096/1‑89 (MQ=255)
                                                  aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGAcc                                     <  2:261883/90‑1 (MQ=255)
                                                        cACCATGCCTCTCCGGCGTATATCGGCGGGGGGGCCTCCCCCTTTCGCGGCCCATTTTTTTAAATTTGTCCGCGGGCGGTGGCCCGtttt                               >  2:505039/1‑90 (MQ=255)
                                                            aTGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCCCATTGCGCGGCGCCTTTTTTAAAATTGGGCTGCGCGCGGTGGCCTGttttttaa                           >  1:128928/1‑86 (MQ=255)
                                                            aTGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCAATTCGCGACGCATTTTTTAAAAATGGCCTGCGCGCGGTGACCCGttttttaa                           >  1:174115/1‑86 (MQ=255)
                                                                ctctCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAactac                       <  2:191060/90‑1 (MQ=255)
                                                                  ctcCGGCGTATATCGGCGGGGCGCGCTCCCCCATGCCCGACGCCCTTTTTAAAATTTGGCTGCGGGCCGGGGCCTGTTTTTATaatactt                     >  2:344675/1‑89 (MQ=255)
                                                                     cGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGAcc                  <  2:175132/90‑1 (MQ=255)
                                                                         gTATATCGGCGGGGGGCTCTCCCCCTTTCGCGGCGCCTTTTTTTAAATTTGGCCGCGCGCGGGGACCCGTTTTTTAAATACTGACCCTaa              >  2:496208/1‑90 (MQ=255)
                                                                             atCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTaacaac          <  2:388096/90‑1 (MQ=255)
                                                                               cggcggGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACgg        <  2:174115/90‑1 (MQ=255)
                                                                                     gTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCaa  <  1:97074/90‑1 (MQ=255)
                                                                                      |                                                                                         
GCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAA  >  NZ_CP009273/3012836‑3013010

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: