Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273502,4140TC66.7% 20.2 / 22.6 27D42G (GAT→GGT) ybaKCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK
Reads supporting (aligned to +/- strand):  ref base T (6/3);  new base C (14/4);  total (20/7)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCCAGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACG  >  NZ_CP009273/502327‑502498
                                                                                       |                                                                                    
tCCAGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGAccc                                                                                    <  1:77649/90‑1 (MQ=255)
   aGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGa                                                                                 >  2:86835/1‑90 (MQ=255)
       gACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCa                                                                             >  2:402223/1‑90 (MQ=255)
                             cACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTc                                                       >  1:157160/1‑90 (MQ=255)
                             cACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTc                                                       >  2:691752/1‑90 (MQ=255)
                             cACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTc                                                       >  1:630603/1‑90 (MQ=255)
                                         tttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAtt                                           >  2:383094/1‑90 (MQ=255)
                                          ttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTg                                          >  2:429642/1‑90 (MQ=255)
                                             atatCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                                       <  1:402223/90‑1 (MQ=255)
                                             atatCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                                       <  2:488704/90‑1 (MQ=255)
                                               atCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTc                                     <  2:316897/90‑1 (MQ=255)
                                                tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTc                                     <  1:791177/89‑1 (MQ=255)
                                                tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTc                                     >  2:791177/1‑89 (MQ=255)
                                                tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCa                                    >  1:639233/1‑90 (MQ=255)
                                                  aCCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGc                                  >  2:706811/1‑90 (MQ=255)
                                                  aCCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGc                                  >  2:559369/1‑90 (MQ=255)
                                                     gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCgg                               >  2:785769/1‑90 (MQ=255)
                                                     gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCgg                               >  1:136376/1‑90 (MQ=255)
                                                         aCTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCg                           <  2:237489/90‑1 (MQ=255)
                                                          cTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGt                          >  1:822723/1‑90 (MQ=255)
                                                            gCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCACCGCCAAAATTGGTTTCAGCCGGATCGTGc                        >  1:397861/1‑90 (MQ=255)
                                                                      gcGTTTTGTAGACCTGGTCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAt              >  2:721038/1‑90 (MQ=255)
                                                                      gcGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAt              <  2:340729/90‑1 (MQ=255)
                                                                      gcGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAt              >  2:838778/1‑90 (MQ=255)
                                                                        gTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATgg            >  1:408602/1‑90 (MQ=255)
                                                                        gTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATgg            >  1:216566/1‑90 (MQ=255)
                                                                                  ccTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACg  >  1:6336/1‑90 (MQ=255)
                                                                                    tGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGt                    <  2:365252/70‑1 (MQ=255)
                                                                                       |                                                                                    
TCCAGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACG  >  NZ_CP009273/502327‑502498

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