Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,454,3620TG56.2% 10.9 / 16.0 16V351G (GTG→GGG) paaC3‑hydroxyacyl‑CoA dehydrogenase PaaC
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base G (9/0);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 8.74e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-03
Rejected as consensus: Frequency below/above cutoff threshold.

TTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAAC  >  NZ_CP009273/1454278‑1454443
                                                                                    |                                                                                 
ttgaGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGGAGTGATCGATAAAATGGCGGGGAAAggggtgg                                                                              >  1:311402/1‑89 (MQ=255)
        caAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGGGACCATtgct                                                                      >  2:335925/1‑90 (MQ=255)
            gCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGGGGGGACCATtgctgctg                                                                  >  2:549197/1‑90 (MQ=255)
              gagaCGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGGGGTGACCATTGCTGCTGCa                                                                >  1:131328/1‑90 (MQ=255)
                    gCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTg                                                          <  2:186476/90‑1 (MQ=255)
                    gCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGGGGTGACCATTGCTGCTGCAGCGGTg                                                          >  2:307690/1‑90 (MQ=255)
                              tGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGGGGGGACCATTGCTGCTGCAGCGGGGAACCCGGACt                                                >  2:548883/1‑90 (MQ=255)
                                             cACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAgg                                 <  1:270778/90‑1 (MQ=255)
                                                        gtAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTAc                      <  1:380832/90‑1 (MQ=255)
                                                        gtAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTAc                      <  1:429219/90‑1 (MQ=255)
                                                               tCGATAAAATGGCGGGGAAGGGGGGGACCATTGCTGCTGCAGCGGGGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAAc               >  2:461650/1‑90 (MQ=255)
                                                                  aTAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAAcagc            <  1:548883/90‑1 (MQ=255)
                                                                    aaaaTGGCGGGCAAGGGGGGGACCATTGCTGCTGCAGCGGTGAAACCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCCAcagccg          >  1:285812/1‑88 (MQ=255)
                                                                            gggCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAAc  <  2:73821/90‑1 (MQ=255)
                                                                            gggCAAGGGGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGGCCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAAc  >  1:347704/1‑90 (MQ=255)
                                                                            gggCAAGGGGGGGGCCATTGCTGCTGCAGCGGGGAACCCGGACTCAGGGACCCGCCAGGCCATTTATTACCTGCAACAGCAGGGCAAAAc  >  1:441669/1‑90 (MQ=255)
                                                                            gggCAAGGGGGGGACCATTGCTGCTGCAGCGGGGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAAc  >  1:110673/1‑90 (MQ=255)
                                                                                    |                                                                                 
TTGAGACACAAGGCGAGACGGCACAGGCGCTGGCAATACGACTGGCACGCCCGGTGGTAGTGATCGATAAAATGGCGGGCAAGGTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAAC  >  NZ_CP009273/1454278‑1454443

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