Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,367,4000AC56.5% 17.0 / 24.5 23G361G (GGT→GGGmenEo‑succinylbenzoate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (7/3);  new base C (0/13);  total (7/16)
Fisher's exact test for biased strand distribution p-value = 4.89e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-04
Rejected as consensus: Frequency below/above cutoff threshold.

AACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAG  >  NZ_CP009273/2367314‑2367485
                                                                                      |                                                                                     
aaCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAccc                                                                                    <  1:152499/90‑1 (MQ=255)
aaCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAccc                                                                                    <  2:602653/90‑1 (MQ=255)
      ttGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGGGAAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCg                                                                              >  2:328353/1‑90 (MQ=255)
        gTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGcc                                                                            <  2:574934/90‑1 (MQ=255)
          cGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc                                                                          <  2:214487/90‑1 (MQ=255)
           ggCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACt                                                                         <  1:418767/90‑1 (MQ=255)
            gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTg                                                                        <  2:153780/90‑1 (MQ=255)
                      cGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt                                                              <  1:191293/90‑1 (MQ=255)
                                cTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACg                                                    >  2:240779/1‑90 (MQ=255)
                                   cAAACCCGCAGGATGTGCAGCAATTACGCGCTCGCCTTCTTCCGGCTGAATCCCCCCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc                                                 <  2:567096/90‑1 (MQ=255)
                                   cAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc                                                 <  2:106726/90‑1 (MQ=255)
                                         cGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAAt                                           >  1:324393/1‑90 (MQ=255)
                                          gCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATg                                          >  1:286991/1‑90 (MQ=255)
                                             ggATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc                                       <  2:244189/90‑1 (MQ=255)
                                                  tgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt                                  <  1:247995/90‑1 (MQ=255)
                                                    aagcaATTGCGCGGTCGGTTTCTTCCGGGTGAATCCCCTCGCCCCCAATGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGc                                <  1:374385/89‑1 (MQ=255)
                                                     agcaATTACGCGCTCGCCTTCTTCCGGCTGAATCCCTTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGcc                               <  1:55377/90‑1 (MQ=255)
                                                           tACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATg                         >  1:651580/1‑90 (MQ=255)
                                                               cgcTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtc                     <  2:579650/90‑1 (MQ=255)
                                                                          ttcCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATc          <  1:62134/90‑1 (MQ=255)
                                                                               gCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATcgcgcg     >  1:490726/1‑90 (MQ=255)
                                                                                  gAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag  <  1:4223/90‑1 (MQ=255)
                                                                                  gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag  >  1:598948/1‑90 (MQ=255)
                                                                                      |                                                                                     
AACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAG  >  NZ_CP009273/2367314‑2367485

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