Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,812,8001.C68.0% 48.3 / 27.4 25intergenic (‑43/+62)celF/chbR6‑phospho‑beta‑glucosidase/transcriptional regulator ChbR
Reads supporting (aligned to +/- strand):  ref base . (4/4);  new base C (7/10);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCA  >  NZ_CP009273/1812718‑1812877
                                                                                   |                                                                             
tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                        >  2:691038/1‑90 (MQ=255)
tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGccccccc                                                                        >  1:373612/1‑90 (MQ=255)
 cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGcccccccc                                                                       <  1:691038/90‑2 (MQ=255)
  ccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCgg                                                                      <  1:568994/90‑1 (MQ=255)
                 acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATGAATAATCCa                                                        <  1:211785/90‑1 (MQ=255)
                 acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa                                                        <  2:635024/90‑1 (MQ=255)
                 acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa                                                        <  1:682272/90‑1 (MQ=255)
                   gacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCa                                                     >  2:69602/1‑90 (MQ=255)
                     cTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa                                                    >  2:300929/1‑90 (MQ=255)
                     cTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAAc                                                   <  1:28868/90‑1 (MQ=255)
                         tAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAgg                                                <  2:299418/90‑1 (MQ=255)
                           aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                              <  1:832222/90‑1 (MQ=255)
                               cTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAg                                                 >  1:358018/1‑83 (MQ=255)
                               cTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAg                                                 <  2:358018/83‑1 (MQ=255)
                               cTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGat                                          <  2:194562/90‑1 (MQ=255)
                                      aTAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAt                                  <  1:613965/90‑1 (MQ=255)
                                       tAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAt                                  >  1:701207/1‑90 (MQ=255)
                                       tAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAt                                  <  1:267562/90‑1 (MQ=255)
                                          tttCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGt                               <  1:838784/90‑1 (MQ=255)
                                                 cTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtat                        >  2:77349/1‑90 (MQ=255)
                                                   tCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAAc                     >  2:224184/1‑90 (MQ=255)
                                                           gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtata                       >  1:208022/1‑81 (MQ=255)
                                                           gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtata                       <  2:208022/81‑1 (MQ=255)
                                                                cTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                              <  2:222235/52‑1 (MQ=255)
                                                                cTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGAc                                              >  1:222235/1‑52 (MQ=255)
                                                                     atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt    >  1:329703/1‑90 (MQ=255)
                                                                     atCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTc   >  2:679051/1‑90 (MQ=255)
                                                                       cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCa  >  1:732713/1‑90 (MQ=255)
                                                                                   |                                                                             
TCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCA  >  NZ_CP009273/1812718‑1812877

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