Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,738,4480AC53.8% 7.4 / 13.7 13V84G (GTG→GGG) trmDtRNA (guanosine(37)‑N1)‑methyltransferase TrmD
Reads supporting (aligned to +/- strand):  ref base A (6/0);  new base C (0/7);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 5.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCC  >  NZ_CP009273/2738370‑2738533
                                                                              |                                                                                     
aGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCt                                                                            >  2:488266/1‑90 (MQ=255)
                 cGTTGCCAGTTCGCTGACGCCCGCTTGATCAGGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCttt                                                           >  2:41708/1‑90 (MQ=255)
                     gCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTtgct                                                       <  2:183900/90‑1 (MQ=255)
                                     ccGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTccc                                       >  1:518762/1‑90 (MQ=255)
                                         ttGATCAAGCTTGCGTCCCTGTGGTGACAAATAAACCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCaa                                   <  2:202521/90‑1 (MQ=255)
                                                  cTTGGGTCCTGGTGGGAACAGAAAAACCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc                          <  2:300778/90‑1 (MQ=255)
                                                  cTTGCGTCCTTGTGGTGACAGAAAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc                          <  2:554205/90‑1 (MQ=255)
                                                     gCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcat                       >  2:504407/1‑90 (MQ=255)
                                                        ccccTGTGGTGACAGAAAAATCCCCTTTGCGCTTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACcatcat                    <  2:470662/89‑1 (MQ=255)
                                                          ccTGTGGTGACAGATAAATCCCCTTGGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTa                  <  2:163811/90‑1 (MQ=255)
                                                           cTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTaa                 >  1:114351/1‑90 (MQ=255)
                                                                      aaTAAATCCCCTTTCCCCTTCCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTc      <  2:459737/89‑1 (MQ=255)
                                                                          aaTCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgcc  >  2:141523/1‑90 (MQ=255)
                                                                              |                                                                                     
AGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCC  >  NZ_CP009273/2738370‑2738533

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