Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273296,4220GC22.2% 170.3 / 13.3 63P75R (CCG→CGG) paoCaldehyde oxidoreductase molybdenum‑binding subunit PaoC
Reads supporting (aligned to +/- strand):  ref base G (13/36);  new base C (14/0);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 5.37e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGG  >  NZ_CP009273/296335‑296511
                                                                                        |                                                                                         
gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTGATGAAAGCCAGGGCGCCCCgc                                                                                         >  2:373127/1‑90 (MQ=255)
gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCCGTACGCCCCgc                                                                                         >  1:63175/1‑90 (MQ=255)
gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAACCCGGACGCCCCgc                                                                                         >  1:34455/1‑90 (MQ=255)
  ggggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACACCCCGTACGCCCCgcgc                                                                                       >  1:410783/1‑90 (MQ=255)
  ggggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGACAGCCCGTACGCCCCgcgc                                                                                       >  2:41554/1‑90 (MQ=255)
     gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCttt                                                                                    >  1:142596/1‑90 (MQ=255)
     gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCGGACGCCCCGCGCctt                                                                                    >  1:435246/1‑90 (MQ=255)
       cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTACGCCCGGCGCttttt                                                                                  >  1:142431/1‑90 (MQ=255)
       cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCttttt                                                                                  >  1:15777/1‑90 (MQ=255)
       cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGGATGGACGCCACTACAGCCCGGGGTtttt                                                                                   >  2:231980/1‑90 (MQ=255)
       cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCACAACCCCCCGGGCttttt                                                                                  >  2:173975/1‑90 (MQ=255)
       cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAGCCAGTACGCCCGGCGCttttt                                                                                  >  1:50159/1‑90 (MQ=255)
        gccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGAAAGCCCGTACGCCCCGCGCcttttt                                                                                 >  2:76133/1‑89 (MQ=255)
              cAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCgcggc                                                                           <  1:439328/90‑1 (MQ=255)
                  cTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTcc                                                                       <  1:173975/90‑1 (MQ=255)
                  cTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTcc                                                                       <  2:15777/90‑1 (MQ=255)
                  cTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTcc                                                                       <  1:41554/90‑1 (MQ=255)
                   tGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTACGCCCCGCGCTTTTTTTGCGGGGGCCg                                                                      >  1:351674/1‑90 (MQ=255)
                        gtgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCa                                                                 <  1:47938/90‑1 (MQ=255)
                        gtgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCa                                                                 <  1:298927/90‑1 (MQ=255)
                          gtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGgc                                                                                         <  2:286277/64‑1 (MQ=255)
                          gtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGgc                                                                                         >  1:286277/1‑64 (MQ=255)
                             tttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGc                                                            <  2:286641/90‑1 (MQ=255)
                               tGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCgg                                                          <  2:217654/90‑1 (MQ=255)
                               gacagCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCgg                                                          <  2:310148/86‑1 (MQ=255)
                                 tCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTg                                                        <  2:142431/90‑1 (MQ=255)
                                  cGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTACGCCCGGCGCCTTTTTCGCGGCGGCCGGATCAAggggggggg                                                       >  1:338562/1‑87 (MQ=255)
                                     cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGc                                                    <  2:240628/90‑1 (MQ=255)
                                     cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGc                                                    <  2:127358/90‑1 (MQ=255)
                                     cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGc                                                    <  2:123618/90‑1 (MQ=255)
                                     cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGc                                                    <  1:149457/90‑1 (MQ=255)
                                      tttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCg                                                   <  1:91115/90‑1 (MQ=255)
                                              gTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTTGCGGGGGCCGGATCAAGGGGGGGGGGGGGGCCttttg                                           >  2:149457/1‑88 (MQ=255)
                                              gTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGACCCCCCGCGCTTTTTGCGCGGCGGCCGTATCAAGGGGGGGGGGGGGGCCttttg                                           >  1:326079/1‑88 (MQ=255)
                                                gCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCa                                         <  1:140894/90‑1 (MQ=255)
                                                gCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCa                                         <  2:149492/90‑1 (MQ=255)
                                                  cccGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAAt                                       <  1:463312/90‑1 (MQ=255)
                                                  cccGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAAt                                       <  2:153717/90‑1 (MQ=255)
                                                   ccGGCGTTACTGGCGGTAAATACAACCAGTACCCCCCGCGCCTTTTTCGCGGCGGCCGCATCAAAGGGGGGGGGGGGGCCTTTTGCAATg                                      >  2:490639/1‑90 (MQ=255)
                                                     ggCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGc                                    <  2:410783/90‑1 (MQ=255)
                                                      gCGTTACTGGCGGTAATGACAGCCCGTACGCCCCGCGCCTTTTTCGCGGCGGCCGTATCAAGGGGGGGGGGGCGGCCTTTTGCAATgggg                                   >  2:85430/1‑88 (MQ=255)
                                                           aCTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAAcc                              <  1:408942/90‑1 (MQ=255)
                                                              ggCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac                           <  2:472125/90‑1 (MQ=255)
                                                              ggCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac                           <  1:76133/90‑1 (MQ=255)
                                                              ggCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac                           <  2:376003/90‑1 (MQ=255)
                                                              ggCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac                           <  2:34455/90‑1 (MQ=255)
                                                              ggCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgac                           <  2:351674/90‑1 (MQ=255)
                                                                cGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAg                                                               <  2:147904/52‑1 (MQ=255)
                                                                cGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAg                                                               >  1:147904/1‑52 (MQ=255)
                                                                   tAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat                      <  1:490639/90‑1 (MQ=255)
                                                                   tAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat                      <  1:313742/90‑1 (MQ=255)
                                                                   tAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGatat                      <  2:332341/90‑1 (MQ=255)
                                                                    aaTGACAGCCAGTACGCCCGCCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGTCGTCCTTTGGCAATGGCGGAACCGACGAtata                     <  2:204492/90‑1 (MQ=255)
                                                                        aCAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCa                 <  1:492055/90‑1 (MQ=255)
                                                                          aGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATa               >  2:283457/1‑90 (MQ=255)
                                                                           gCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAg              >  2:281256/1‑90 (MQ=255)
                                                                              aGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcg           >  1:335083/1‑90 (MQ=255)
                                                                               gTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcgg          >  2:495431/1‑90 (MQ=255)
                                                                               gTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcgg          <  1:73106/90‑1 (MQ=255)
                                                                                 aCGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGgcggcg        <  2:118711/90‑1 (MQ=255)
                                                                                  cGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGt       >  2:180857/1‑90 (MQ=255)
                                                                                  cGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGt       <  2:50159/90‑1 (MQ=255)
                                                                                    gCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGtt      >  2:183656/1‑90 (MQ=255)
                                                                                        ggCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcg                             >  1:15801/1‑63 (MQ=255)
                                                                                        ggCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcg                             <  2:15801/63‑1 (MQ=255)
                                                                                        ggCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTgggg  >  1:50049/1‑90 (MQ=255)
                                                                                        |                                                                                         
GTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGG  >  NZ_CP009273/296335‑296511

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