Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273645,9160TG52.2% 26.4 / 34.7 23T11P (ACT→CCT) citEcitrate (pro‑3S)‑lyase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (2/9);  new base G (12/0);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 6.73e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATT  >  NZ_CP009273/645780‑646050
                                                                                                                                        |                                                                                                                                       
cgggcggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGttt                                                                                                                                       <  2:169606/139‑1 (MQ=255)
                 aCGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGAGGTTTTTGGTTTTGCCGCGag                                                                                                                      >  2:70486/1‑138 (MQ=255)
                   gCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTCGTTTTTGCAGCGAAGc                                                                                                                    >  1:365663/1‑139 (MQ=255)
                       tGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGAGGTTTTCGTTGTTGCAGCGAAGCGGa                                                                                                                 >  1:353229/1‑138 (MQ=255)
                                 tCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTCGTTGTTGCAGCGTAACGGAAACCATTGGc                                                                                                       >  2:216772/1‑138 (MQ=255)
                                 tCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTTGTTGTTGTAGCGAAGAGGGAAACAtttgg                                                                                                       >  2:269087/1‑135 (MQ=255)
                                 tCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGCAGCGGAGCGGAAATCATTGGc                                                                                                       >  2:353356/1‑138 (MQ=255)
                                   ttCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGGGGTTTTCGTTGTTGCCGCGGAAAGGAAAACATTGGCAAt                                                                                                    >  2:131244/1‑139 (MQ=255)
                                   ttCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGGTGTTGCAGCGAAGCGGGAAACATTGGCCAt                                                                                                    >  2:186644/1‑139 (MQ=255)
                                     cGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAAt                                                                                                    <  2:359743/137‑1 (MQ=255)
                                     cGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAAt                                                                                                    <  1:353356/137‑1 (MQ=255)
                                         gTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGTTTTACGTTGTTTCCGCGAAGCGGAAATCATTGGCAATCCTccc                                                                                              >  1:318009/1‑139 (MQ=255)
                                              aCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGc                                                                                         <  1:85794/139‑1 (MQ=255)
                                                                    gATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGgcg                                                                   <  2:269123/138‑1 (MQ=255)
                                                                                                   gCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTa                                    <  1:269087/138‑1 (MQ=255)
                                                                                                       cAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATg                               <  2:370209/139‑1 (MQ=255)
                                                                                                        aGGCACAAACAACATGCTGCGGCGGGTGCGCGGTTTTTGTTGTTGCCGCGAAGCGGAAATCCTTGGCAATCCTCCCACGGGCGAGGCGGGGTACCGCTGGGGCGTGGCCGCAGGGCTTCCCGTCGGGGACGGaaaaag                               >  2:32569/1‑136 (MQ=255)
                                                                                                        aGGCACAAACAACATGCTGCGGCGGGGGCGAGGTTTTCGTTGTTGCAGCGAAGCGGGAAACATTTGCAATCCTCCCCTTGCCGAGGCCGGGAACCCCTGGCGCGGGCCCGCCGGGGTTCCCCTCGGGCACCTAAAATg                               >  2:387880/1‑138 (MQ=255)
                                                                                                            acaaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTc                          <  2:318009/139‑1 (MQ=255)
                                                                                                               aacaacATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAg                       <  1:70486/139‑1 (MQ=255)
                                                                                                                          gcggcgGGTGCGAGGTTTACGTTGTTGCAGCGAAGCGGAAATAATTGGCAATCCCCCCCTGGCAGAGCCGGGGGACCGCTGGCGCGTGCCAGACGGGGTTCCCGTCGTGCCCGGAAAATGGAGGCCAGGGCGCCctttt            >  2:341951/1‑137 (MQ=255)
                                                                                                                              cggGGGCGAGTTTTTCGTTGTTGCAGCGAAGCGGAAATAATTGGCAATCCTCCCCTGGCAGAGGCGGGGTACCGCTGGCGCGGGCCAGCAGGGCTTCCAGTCGTGCCCGGAAAATGCAGTCCAGGGCGCCTTTGTCATc        >  2:335097/1‑139 (MQ=255)
                                                                                                                              cggGGGCGAGGTTTTCGTTGTTGCAACGGAGAGGAAAACATTGGGAATCCTCCCAAGGCAGAGCCGGGGGACCCCCGGCGCGGGCCAGCAGGGCTTCCAGTCGGGCGCGTAAAATGCAGTCCAGTGCGCCTTTGGCATc        >  2:14631/1‑139 (MQ=255)
                                                                                                                                    cGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACAtt  <  1:387880/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
CGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATT  >  NZ_CP009273/645780‑646050

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