Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,341,2370TG58.8% 11.1 / 14.9 17V255G (GTG→GGG) rapZRNase adapter RapZ
Reads supporting (aligned to +/- strand):  ref base T (0/7);  new base G (10/0);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 5.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCA  >  NZ_CP009273/3341154‑3341325
                                                                                   |                                                                                        
gCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTatat                                                                                    <  2:172570/90‑1 (MQ=255)
   aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGGGTATATTGc                                                                                 >  1:85745/1‑90 (MQ=255)
        tACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGGGTATATTGCAGAGc                                                                            >  1:65507/1‑90 (MQ=255)
             aTGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGGGTATATTTCAGAGCAACTg                                                                       >  1:57394/1‑90 (MQ=255)
             aTGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCCCTCGGGGGATATTGCAGAGCAACTg                                                                       >  2:17275/1‑90 (MQ=255)
               gCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGc                                                                     <  2:45459/90‑1 (MQ=255)
                                gTAGATACTTGCCGGTCGCCCTGGGGTGTGCCGGCGGGCAGCACCGTTCGGTGTATATTGCAGAGCACCTGGCCCAATTCCTCCGCTcgc                                                    <  2:224291/90‑1 (MQ=255)
                                         tGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTaaaa                                           <  1:40456/90‑1 (MQ=255)
                                         tGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTaaaa                                           <  1:36760/90‑1 (MQ=255)
                                              gTCGCCATTGGTTGTACCGGCGGGAAGAACCGGTCGGGGTATATTTCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGGAAAAAAGTc                                      >  1:194561/1‑90 (MQ=255)
                                                cGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCa                                    <  2:237913/90‑1 (MQ=255)
                                                            tACCGGCGGGAAGGACCGTTCGGGGGATATTGCAAAGCAACTGGCAGACTACTTCCCCCCGCGCGGGGAAAAAGGCCCGGCACGCCCTcc                        >  2:173517/1‑89 (MQ=255)
                                                            tACCGGCGGGAAGCACCGTTCGGGGTATATTGCAGAGCAACTGGCAGACCACTTCCGCTCGCGCGGGAAAAACGTCCAGTCACGCCCTCg                        >  1:77189/1‑90 (MQ=255)
                                                            tACCGGCGGGAAGCACCCGTCGGGGTATATTTCAGAGCAACTGGCAGACTACTTCCCCCCCCGCGGGAAAAAAGTCCAGTCACGCCCTCg                        >  2:3777/1‑90 (MQ=255)
                                                                ggcggGAAGCACCGTTCGGGGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGGAAAAAAGTCCAGTCACGCCATCGTACg                    >  2:219302/1‑90 (MQ=255)
                                                                 gcggGAAGCACCGTTCGGGGTATATTGCAGAGCAACTGGCAGACTACTTCCGCCCGCGCGGGAAAAAAGTCCCGTCACGCCACCGTACGc                   >  2:223853/1‑90 (MQ=255)
                                                                       aGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGaaa             <  2:52930/90‑1 (MQ=255)
                                                                                  gtgtATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCa  >  1:282330/1‑90 (MQ=255)
                                                                                   |                                                                                        
GCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCA  >  NZ_CP009273/3341154‑3341325

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