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

ACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACC  >  NZ_CP009273/3012833‑3012994
                                                                                         |                                                                         
aCAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                                                                           <  1:19228/90‑1 (MQ=255)
aCAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                                                                           <  2:486361/90‑1 (MQ=255)
  aGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCg                                                                         <  1:474721/90‑1 (MQ=255)
   gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                                                                           <  1:150958/87‑1 (MQ=255)
   gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                                                                           >  2:150958/1‑87 (MQ=255)
   gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGgcgc                                                                        >  2:552010/1‑90 (MQ=255)
   gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCg                                                                         <  2:329975/89‑1 (MQ=255)
   gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCgg                                                                        >  1:219747/1‑89 (MQ=255)
   gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCgg                                                                        >  1:268491/1‑89 (MQ=255)
   gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCgg                                                                        >  1:649381/1‑89 (MQ=255)
   gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCg                                                                         >  1:329975/1‑89 (MQ=255)
     ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTccgtg                                                                      >  1:502085/1‑86 (MQ=255)
     ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGgtc                                                                      >  1:403583/1‑90 (MQ=255)
     ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGctc                                                                      >  2:83732/1‑90 (MQ=255)
     ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGctc                                                                      >  1:94208/1‑90 (MQ=255)
       aCTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGCCGTATATCGGCGGGGCGCtctc                                                                    >  2:652469/1‑90 (MQ=255)
                                  gagCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGt                                        <  1:637838/90‑1 (MQ=255)
                                                 ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGGGCGCTCCCCCCTTCGCGGCGCCCTTTTTTAAAATGGGCTGggcgcggg                         >  1:387309/1‑88 (MQ=255)
                                                 ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCATTTCGCGACGCCTTTTTTAAAATTTGGCTGcgcgcggg                         >  2:157527/1‑88 (MQ=255)
                                                 ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCAATTCGCGGCGCCTTTTTTAAAATTTGGCcgcgcgcggg                         >  2:404397/1‑88 (MQ=255)
                                                 ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCTCCCCATTTCGCGGCCCCTTTTTTAAAATTTGGCTGcgcgcggg                         >  2:162835/1‑88 (MQ=255)
                                                     aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGAcc                     <  2:543817/90‑1 (MQ=255)
                                                       cGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCTCCCCCTTTCGCGCCGCATTTTTTAAAATTGGGCTGCGCGCGCGGAcccg                   >  2:440814/1‑88 (MQ=255)
                                                           cACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGtttt               <  2:387591/90‑1 (MQ=255)
                                                              cATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGaa            <  2:231645/90‑1 (MQ=255)
                                                                   ctctCCGGCGTATATCGGCGGGTCGCTCCCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAactac       <  2:104998/90‑1 (MQ=255)
                                                                        cGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGAcc  <  1:399765/90‑1 (MQ=255)
                                                                                         |                                                                         
ACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACC  >  NZ_CP009273/3012833‑3012994

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