Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,341,2370TG54.5% 34.4 / 42.1 33V255G (GTG→GGG) rapZRNase adapter RapZ
Reads supporting (aligned to +/- strand):  ref base T (4/11);  new base G (18/0);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 7.05e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CTTGAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCA  >  NZ_CP009273/3341149‑3341325
                                                                                        |                                                                                        
cTTGAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGgtg                                                                                         <  1:1112083/90‑1 (MQ=255)
cTTGAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGgtg                                                                                         <  1:1204168/90‑1 (MQ=255)
cTTGAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGgtg                                                                                         <  2:658965/90‑1 (MQ=255)
   gAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGGACCGGTCGGGGtat                                                                                      >  2:988976/1‑90 (MQ=255)
     gCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGGACCCTTCGGGGTatat                                                                                    >  2:1405318/1‑90 (MQ=255)
        aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGGACCCGTCGGGGTATAtttc                                                                                 >  2:968339/1‑88 (MQ=255)
        aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGCACCCGTCCGGGTATAtttc                                                                                 >  2:689460/1‑88 (MQ=255)
        aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGgtg                                                                                         <  2:589757/82‑1 (MQ=255)
        aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGgtg                                                                                         >  1:589757/1‑82 (MQ=255)
        aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGGTTGGGGGGTATTGc                                                                                 >  1:275707/1‑90 (MQ=255)
        aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGGTCCGGGTATAtttc                                                                                 >  1:1597963/1‑88 (MQ=255)
               ccTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCaa                                                                          <  2:55913/90‑1 (MQ=255)
                cTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGCACCGTTCGGGGGTTATTTCAGAGCaaa                                                                         >  1:373642/1‑89 (MQ=255)
                  aTGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGCACCGTTCGGGGTTTATTGGAGAGCAAATg                                                                       >  2:753021/1‑90 (MQ=255)
                  aTGCTGGAAACCAACAACCGTAGCTACCTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCCTTCGGGGGATATTGCAGAGCAAATg                                                                       >  2:1234762/1‑90 (MQ=255)
                    gCTGGACAGCAACCACAGTAGCTACTTGACGGTGGCAATTGGTTGTACCGGCGGGAAGGCACGTTCGCTGTATATTGCAGAGCAACTGGc                                                                     <  1:476728/90‑1 (MQ=255)
                    gCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGCACCCTTCGGGGTATATTTCAGAGCAACTGGc                                                                     >  2:1245113/1‑90 (MQ=255)
                    gCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGc                                                                     <  2:158430/90‑1 (MQ=255)
                            ccAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGCACCCGTCGGGGTATATTTCAGAGCAACTGGCAGactact                                                             >  2:715828/1‑90 (MQ=255)
                            ccAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCCGTCGGGGGATATTTGAGAGCAACTGGCAGgctact                                                             >  1:404411/1‑90 (MQ=255)
                                        cagaCTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTcgcgcg                                                 <  1:891841/87‑1 (MQ=255)
                                         cTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGAACCCCTCGGGGTATATTGCAGAGCAACTGGCAGGCTACTTCCGCCCGCGCgg                                                >  1:718611/1‑90 (MQ=255)
                                            cTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGGGGGTATTTCAGAGCAAATGGCCGACTAATTCCGCTCCCGCGGGaa                                             >  2:1459841/1‑90 (MQ=255)
                                                     cGCCATTGGTTGTACCGGCGGGAAGCACCCTTCGGGGTATATTTCAGAGCAACTGGCAGACTACTTCCCCTCCCGCGGGAAAAAAGTccc                                    >  2:1588772/1‑89 (MQ=255)
                                                          ttGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCAc                               <  1:1405318/90‑1 (MQ=255)
                                                          ttGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCAc                               <  1:968339/90‑1 (MQ=255)
                                                            ggTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGc                             <  2:1040655/90‑1 (MQ=255)
                                                                 tACCGGCGGGGAGCACCGCTCGGTGGATATTTCAGAGGAAATGGCCGACTACTTCCCCCCGCGCGGGAAAAACCACCAGTCACGCCATCg                        >  1:1397364/1‑90 (MQ=255)
                                                                 tACCGGCGGGAAGCACCGTTCGGGGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGGAAAAAAGTCCCGTCACGCCATCg                        >  2:459840/1‑90 (MQ=255)
                                                                   ccGGCGGGAAGCACCGTTCGGGGTATATTGCAGAGCAACTGGCAGACTACTTCCCCTCGCGCGGGAAAAAAGGCCAGTCACGCCATCGTa                      >  2:1292443/1‑90 (MQ=255)
                                                                     ggcggGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACg                    <  1:1245113/90‑1 (MQ=255)
                                                                      gcggGAAGCACCGTTCGGGGTATATTTCAGAGCAACTTGCAGACCACTTCCGCTCCCGCGGGAAAAAAATCCAGTCACGCCATCGTACGc                   >  1:1089854/1‑90 (MQ=255)
                                                                            aGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGaaa             <  2:1077862/90‑1 (MQ=255)
                                                                            aGCACCATTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGAAAAAACGTCCAGTCACGCCATCGTACGCTGGaaa             >  2:615336/1‑90 (MQ=255)
                                                                                 cGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGt        <  1:689460/90‑1 (MQ=255)
                                                                                  gTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTa       >  2:360246/1‑90 (MQ=255)
                                                                                  gTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTa       >  2:660955/1‑90 (MQ=255)
                                                                                   ttCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTaa      <  2:373642/90‑1 (MQ=255)
                                                                                       gtgtATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCa  >  2:684439/1‑90 (MQ=255)
                                                                                        |                                                                                        
CTTGAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCA  >  NZ_CP009273/3341149‑3341325

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