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

GACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCG  >  NZ_CP009273/2738359‑2738537
                                                                                         |                                                                                         
gACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCa                                                                                           >  1:414546/1‑90 (MQ=255)
       cacCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTgc                                                                                    <  2:236984/90‑1 (MQ=255)
                    ttttGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACccgcc                                                                       <  1:71540/90‑1 (MQ=255)
                            cGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCttt                                                               >  2:508867/1‑90 (MQ=255)
                                gCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCCCCCGCCGCGGCTTTtgct                                                           <  2:414546/90‑1 (MQ=255)
                                        gcgaaCGCCCCTTTAATCAAGTTGGCCCCCGTGTGGTCACAGAAAAACCCCTTTTCCGCCTCCCCCCGCCGCGGCTTTTGCTGCATGAAt                                                   <  2:507043/86‑1 (MQ=255)
                                                   cTTGATCAAGCTTGCGCCCTTGTGGTGACAGATAAACCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCa                                        <  2:504308/90‑1 (MQ=255)
                                                     tGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAg                                      <  1:595503/90‑1 (MQ=255)
                                                             cTTGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc                              <  1:508197/90‑1 (MQ=255)
                                                             cTTGCGCCCCTGTGGTAACAGAAAAACCCCCTTTCCCCCTCCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCAc                              <  2:588629/90‑1 (MQ=255)
                                                               tGCGTCCCTGTGGTGACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACca                            <  1:374603/90‑1 (MQ=255)
                                                                             gACAGATAAATCCCCTTTGCGCCTTCCCCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcccc              <  1:137298/90‑1 (MQ=255)
                                                                             gACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATcccc              >  2:531409/1‑90 (MQ=255)
                                                                               cAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCgg            >  1:327695/1‑90 (MQ=255)
                                                                                         aCCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTccgccgccg  >  2:350548/1‑90 (MQ=255)
                                                                                         |                                                                                         
GACCGCACACCAGAATCAATTTTTGATTCGTTGCCAGTTCGCTGACGCCCGCTTGATCAAGCTTGCGTCCCTGTGGTGACAGATAAATCACCTTTGCGCCTTCACCCGCCGCGGCTTTTGCTGCATGAATGGCGTCCCGCAAGGGTTGCACCATCATTAACATCCCCGGTCCGCCGCCG  >  NZ_CP009273/2738359‑2738537

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