Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273295,2160GC37.8% 48.2 / 66.5 45A477G (GCC→GGC) paoCaldehyde oxidoreductase molybdenum‑binding subunit PaoC
Reads supporting (aligned to +/- strand):  ref base G (8/20);  new base C (6/11);  total (14/31)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCAC  >  NZ_CP009273/295127‑295303
                                                                                         |                                                                                       
caGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGc                                                                                        <  1:385503/91‑1 (MQ=255)
 aGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGcc                                                                                       <  1:717958/91‑1 (MQ=255)
 aGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGcc                                                                                       <  1:477765/91‑1 (MQ=255)
   aaCCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAg                                                                                     <  1:26956/91‑1 (MQ=255)
    aCCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGa                                                                                    >  2:400500/1‑91 (MQ=255)
    aCCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGa                                                                                    >  1:304000/1‑91 (MQ=255)
    aCCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGa                                                                                    >  1:598308/1‑91 (MQ=255)
    aCCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGa                                                                                    <  1:453141/91‑1 (MQ=255)
    aCCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGa                                                                                    >  2:20967/1‑91 (MQ=255)
    aCCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGa                                                                                    >  1:495822/1‑91 (MQ=255)
      ccGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAAt                                                                                  <  1:149837/91‑1 (MQ=255)
      ccGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAAt                                                                                  <  2:494274/91‑1 (MQ=255)
      ccGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAAt                                                                                  <  2:337798/91‑1 (MQ=255)
      ccGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGAAt                                                                                  <  1:320237/91‑1 (MQ=255)
                  ggAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGcc                                                                      <  2:136294/91‑1 (MQ=255)
                  ggAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGAATGGTGTAGCTGcc                                                                      <  1:20967/91‑1 (MQ=255)
                  ggAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGAATGGTGTAGCTGcc                                                                      <  1:586260/91‑1 (MQ=255)
                  ggAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGAATGGTGTAGCTGcc                                                                      <  2:332613/91‑1 (MQ=255)
                  ggAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGAATGGTGTAGCTGcc                                                                      <  1:195766/91‑1 (MQ=255)
                   gAAACTGGAATCGCCGATGTGAACCGCATCCTGCTCCAGCGGTACGCCAATCATTTCCGCTGCCTTCTGGGCCATAATGGTGTAGCTGCCg                                                                     <  2:59134/91‑1 (MQ=255)
                       cTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGaa                                                                                   <  1:576486/73‑1 (MQ=255)
                       cTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGaa                                                                                   >  2:576486/1‑73 (MQ=255)
                        tGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGAATGGTGTAGCTGCCGGTGcc                                                                <  2:656738/91‑1 (MQ=255)
                         ggAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCa                                                               <  2:113364/91‑1 (MQ=255)
                         ggAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGAATGGTGTAGCTGCCGGTGCCa                                                               <  2:492175/91‑1 (MQ=255)
                             tCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGt                                                           <  1:714432/91‑1 (MQ=255)
                             tCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGt                                                           <  1:350291/91‑1 (MQ=255)
                              cGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGAATGGTGTAGCTGCCGGTGCCAATGTc                                                          <  1:660494/91‑1 (MQ=255)
                              cGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGAATGGTGTAGCTGCCGGTGCCAATGTc                                                          <  1:599029/91‑1 (MQ=255)
                                      tGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGt                                                  <  2:613836/91‑1 (MQ=255)
                                       gAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGCATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTc                                                 >  2:220444/1‑91 (MQ=255)
                                            gCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGttt                                            >  1:345835/1‑91 (MQ=255)
                                             cAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTc                                           <  1:400500/91‑1 (MQ=255)
                                              aaCCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCt                                          <  1:369699/91‑1 (MQ=255)
                                               aCCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTa                                         >  1:148253/1‑91 (MQ=255)
                                                      ccAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAAc                                  <  2:357355/91‑1 (MQ=255)
                                                      ccAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAAc                                  <  2:721878/91‑1 (MQ=255)
                                                       cAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGCCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACg                                 >  2:280480/1‑91 (MQ=255)
                                                             tACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCg                           >  2:254073/1‑91 (MQ=255)
                                                              aCGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGt                          >  2:740596/1‑91 (MQ=255)
                                                                   aaGCATTTCCGCTGCCGTCTGGCCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGt                     >  2:807027/1‑91 (MQ=255)
                                                                     gCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTc                   <  1:323678/91‑1 (MQ=255)
                                                                        tttCCGCTGCCGTCTGGCCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGCGCCGTTTTGTTCGAg                >  2:125822/1‑91 (MQ=255)
                                                                             gCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGaa           <  2:144465/91‑1 (MQ=255)
                                                                                      tGGCCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCAc  <  1:774168/91‑1 (MQ=255)
                                                                                         |                                                                                       
CAGAACCCGCAGAAACCGGGAAACTGGAATCGCCGAGGTGAACCGCAACCTGCTCCAGCGGTACGCCAAGCATTTCCGCTGCCGTCTGGGCCAGAATGGTGTAGCTGCCGGTGCCAATGTCGGTCATGTCCGTTTCTACGGTAACGGTGCCGTTTTGTTCGAGGTGAACCCGAGCAC  >  NZ_CP009273/295127‑295303

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