Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,008,2300GA69.7% 49.6 / 24.3 33A258A (GCG→GCA) pqiAmembrane integrity‑associated transporter subunit PqiA
Reads supporting (aligned to +/- strand):  ref base G (4/6);  new base A (14/9);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGTTATGGAGCGAAGGCTCTTAT  >  NZ_CP009273/1008092‑1008362
                                                                                                                                          |                                                                                                                                    
aGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTcgcg                                                                                                                                      >  1:9034/1‑139 (MQ=255)
aGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTcgcg                                                                                                                                      <  1:430790/139‑1 (MQ=255)
 ggCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCgcgc                                                                                                                                     <  2:132590/139‑1 (MQ=255)
                cGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCca                                                                                                                      <  2:12396/139‑1 (MQ=255)
                cGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCca                                                                                                                      <  2:574826/139‑1 (MQ=255)
                     ggggCTGCGTTCGTTCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATg                                                                                                                 <  1:225866/139‑1 (MQ=255)
                     ggggCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATg                                                                                                                 >  1:748971/1‑139 (MQ=255)
                                             tACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTg                                                                                         <  2:423805/139‑1 (MQ=255)
                                             tACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTg                                                                                         >  1:791723/1‑139 (MQ=255)
                                               cGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGcc                                                                                       >  2:793660/1‑139 (MQ=255)
                                               cGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGcc                                                                                       >  2:399987/1‑139 (MQ=255)
                                                           ctgctgATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGa                                                                            <  1:793660/138‑1 (MQ=255)
                                                                 atgaACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGAtt                                                                      >  1:745057/1‑138 (MQ=255)
                                                                      cccGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGa                                                                            >  1:58695/1‑127 (MQ=255)
                                                                      cccGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGa                                                                            <  2:58695/127‑1 (MQ=255)
                                                                                cgcgTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCactgc                                                                                                                                  >  1:325261/1‑63 (MQ=255)
                                                                                cgcgTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCactgc                                                                                                                                  <  2:325261/63‑1 (MQ=255)
                                                                                            ccAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGAt                                          >  2:185619/1‑139 (MQ=255)
                                                                                               aaGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGa                                                                            >  1:768469/1‑102 (MQ=255)
                                                                                               aaGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGa                                                                            <  2:768469/102‑1 (MQ=255)
                                                                                               aaGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATtct                                       >  2:833207/1‑139 (MQ=255)
                                                                                                 gggTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTg                                                                             >  2:171916/1‑99 (MQ=255)
                                                                                                 gggTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTg                                                                             <  1:171916/99‑1 (MQ=255)
                                                                                                 gggTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGa                                                                            <  1:737222/100‑1 (MQ=255)
                                                                                                 gggTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGa                                                                            >  2:737222/1‑100 (MQ=255)
                                                                                                    tACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTg                                  <  2:9034/139‑1 (MQ=255)
                                                                                                    tACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGc                                    >  2:469062/1‑137 (MQ=255)
                                                                                                    tACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGc                                    <  1:469062/137‑1 (MQ=255)
                                                                                                           ggCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCAt                           >  1:126398/1‑139 (MQ=255)
                                                                                                            gCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCAtt                          <  2:281463/139‑1 (MQ=255)
                                                                                                            gCGTAGAAACAGCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCAtt                          <  1:185619/139‑1 (MQ=255)
                                                                                                                aGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGt                      >  1:723582/1‑139 (MQ=255)
                                                                                                                       gCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATAtt                                                                                            >  2:108232/1‑62 (MQ=255)
                                                                                                                       gCCTGCAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATAtt                                                                                            <  1:108232/62‑1 (MQ=255)
                                                                                                                            cAGTGGACACTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGTTATGGAGCGaa           <  2:748971/138‑1 (MQ=255)
                                                                                                                                  acacTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGTTATGGAGCGAAGGctct     >  1:829515/1‑138 (MQ=255)
                                                                                                                                    acTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGTTATGGAGCGAAGGCTCTTAt  >  2:661540/1‑139 (MQ=255)
                                                                                                                                    acTCGCACTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGTTATGGAGCGAAGGCTCTTAt  >  2:590279/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                    
AGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGTTATGGAGCGAAGGCTCTTAT  >  NZ_CP009273/1008092‑1008362

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