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

ACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAAC  >  NZ_CP009273/3012833‑3013002
                                                                                         |                                                                                  
aCAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                                                                                    <  1:646280/90‑1 (MQ=255)
  aGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCg                                                                                  <  2:264892/90‑1 (MQ=255)
   gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCgg                                                                                 >  1:588250/1‑89 (MQ=255)
     ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt                                                                                    <  1:534309/85‑1 (MQ=255)
     ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGGGctc                                                                               >  1:483134/1‑90 (MQ=255)
                               tAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTaaa                                                   <  1:216095/90‑1 (MQ=255)
                                            aCATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTgcgc                                      <  1:497667/90‑1 (MQ=255)
                                            aCATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTgcgc                                      <  1:460054/90‑1 (MQ=255)
                                               tCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCATTGCGCGGCGCATTTTTTAAAATTGGGCTGcgcgcg                                   >  1:224445/1‑90 (MQ=255)
                                                 ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGGGCCCTCCCCCCTTCGCGGCGCCTTTTTTTAAAATGGGCTGcgcgcggt                                 >  2:242006/1‑88 (MQ=255)
                                                 ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCATTTTGCGACCCCTTTTTTAAAAATTGGCTGcgggcggg                                 >  1:613487/1‑88 (MQ=255)
                                                  gTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGGGCCCTTCCCATTTCGCGGCGCATTTTTTAAAATTGGCCTGcgcgcgggg                                >  2:497667/1‑87 (MQ=255)
                                                    tATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCCATTCGCGACGCCCTTTTTAAAAATTGGCCGCGCGCGGTGAc                              >  1:82475/1‑90 (MQ=255)
                                                     aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGAcc                             <  2:700885/90‑1 (MQ=255)
                                                     aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGAcc                             <  1:600541/90‑1 (MQ=255)
                                                      tCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCt                            <  2:82475/90‑1 (MQ=255)
                                                      tCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCt                            <  1:795983/90‑1 (MQ=255)
                                                         aTCACCATGCCTCTCCGGCGTATATCGGCGGGGCGGGCTCCCCCTTTCGCCGCGCCTTTTTTTAAAATTGGCCGGGCGGGGTGACCCGtt                         >  2:598096/1‑90 (MQ=255)
                                                         aTCACCATGCCTCTCCGGCGTATATCGGCGGGGCGGCCTCCCCCTTTCGCGGCGCCTTTTTTTAAAATTGGCTGCGCGCGGGGGCCTGtt                         >  2:148338/1‑90 (MQ=255)
                                                             ccATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGa                     <  2:592498/90‑1 (MQ=255)
                                                                     ctcCGGCGTATATCGGCGGGGCGCTCTCCCCATTGCGCGACGCCTTTTTTAAAATTGGCCCGCGCGCGGTGACCCGTTTTTTAACTACTg             >  1:52988/1‑90 (MQ=255)
                                                                       ccGGCGTATATCGGCGGGGCGCTCTCCCCCTTTCGCGGCGCCTTTTTTTAAATTTGGCCGCGCGCGGGGGCCCGTTTTTTAAATACTTAc           >  2:195196/1‑90 (MQ=255)
                                                                       ccGGCGTATATCGGCGGGGCGCCCTCCCCCATTCGCGGCGCCTTTTTTTAAAATTGCCCGCGCGCGGGGGCCCGTTTTTTAAATAATGAc           >  2:178066/1‑90 (MQ=255)
                                                                                atCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTaacaac  <  1:148338/90‑1 (MQ=255)
                                                                                atCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTaacaac  <  1:598096/90‑1 (MQ=255)
                                                                                         |                                                                                  
ACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAAC  >  NZ_CP009273/3012833‑3013002

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: