Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273294,6710AC42.3% 35.6 / 10.6 26Y659D (TAT→GAT) paoCaldehyde oxidoreductase molybdenum‑binding subunit PaoC
Reads supporting (aligned to +/- strand):  ref base A (12/3);  new base C (0/11);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 4.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGG‑TATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCC  >  NZ_CP009273/294533‑294806
                                                                                                                                           |                                                                                                                                       
cgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGG‑TATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATa                                                                                                                                         >  2:594420/1‑139 (MQ=255)
cgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGG‑TATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATa                                                                                                                                         >  2:476726/1‑139 (MQ=255)
cgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGG‑TATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATa                                                                                                                                         >  2:517492/1‑139 (MQ=255)
                    aCGCCGCACAGGCCCAGCTCCCCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGCTCGG‑TATCCTCCAGGAAAATCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTaa                                                                                                                     <  1:101892/139‑1 (MQ=255)
                          cacaGTCCCAGCTCACCGCCCCCTTTGCCCTCCTCCGGGGAGGATATGGGTGGGT‑TATCACCCAGGAAAATCCCCTCCTTTTTTGGAATATCCGAAGGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGt                                                                                                               <  2:386578/139‑1 (MQ=255)
                                         ccGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGG‑TATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt                                                                                                >  1:176463/1‑139 (MQ=255)
                                                       cTTCATCGGGGAGGATATGGGGTCGG‑TATCACCCAGAAAAATCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc                                                                                  <  1:150599/139‑1 (MQ=255)
                                                          cATCGGGGAGGATATGGGGTCGG‑TATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc                                                                                >  1:20334/1‑138 (MQ=255)
                                                                    tatagGGGGTGGG‑TATCACCAAGAAAAATCCCCCCCTGTTTGGGGATACCCGAAGGAACCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg                                                                     <  1:309741/138‑1 (MQ=255)
                                                                        tgggttggg‑tTTCATCAAGGAAAATCCCCTCCTGTTTTGGGTTATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcc                                                                 <  2:479418/139‑1 (MQ=255)
                                                                         tgggTCGGTTCTCACCCAGGAAAACCCCCCCCTTGTTGGGGATCGCGGCAGGCCCCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcc                                                                 <  2:405753/138‑1 (MQ=255)
                                                                          ggctagg‑TATCACCAAGGAAAATCCCCTCCTGTTTTGGAATATCCGCATGACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccca                                                               <  1:199527/134‑1 (MQ=255)
                                                                            gTCGG‑TATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatg                                                             >  2:257515/1‑139 (MQ=255)
                                                                                   atcaccCAGAAAAACCCCCCCCTTTTTTGGTATACCCGCATGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                        <  2:346513/138‑1 (MQ=255)
                                                                                   accacccAGAAAAACCCCCCCCTTTTTGGGGATACCCAAATGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                        <  2:149282/133‑1 (MQ=255)
                                                                                      atcCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt                                                     >  2:434860/1‑138 (MQ=255)
                                                                                            gAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGGCATTgcg                                              >  1:614964/1‑139 (MQ=255)
                                                                                                 tCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCCAt                                         >  2:259611/1‑139 (MQ=255)
                                                                                                      tcctGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                    >  2:460270/1‑139 (MQ=255)
                                                                                                               gggATACCCGCATGCCCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg                           <  1:295006/139‑1 (MQ=255)
                                                                                                                gaatacCCGCATGACCCGCCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                          <  2:348966/137‑1 (MQ=255)
                                                                                                                     tCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                     >  1:348966/1‑139 (MQ=255)
                                                                                                                                 ggCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg                           >  2:138278/1‑121 (MQ=255)
                                                                                                                                 ggCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg                           <  1:138278/121‑1 (MQ=255)
                                                                                                                                  gccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGAt        <  1:476726/136‑1 (MQ=255)
                                                                                                                                        cATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg                                                                     >  2:266633/1‑72 (MQ=255)
                                                                                                                                        cATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg                                                                     <  1:266633/72‑1 (MQ=255)
                                                                                                                                        cATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTcc  >  1:529432/1‑139 (MQ=255)
                                                                                                                                        cATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTcc  >  1:149013/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                       
CGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGG‑TATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCC  >  NZ_CP009273/294533‑294806

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑