Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,485,4430AC56.0% 17.6 / 24.9 25V195G (GTG→GGG) ahrNADPH‑dependent aldehyde reductase Ahr
Reads supporting (aligned to +/- strand):  ref base A (6/5);  new base C (0/14);  total (6/19)
Fisher's exact test for biased strand distribution p-value = 2.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-03
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCTGCTCTTTCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAG  >  NZ_CP009273/4485373‑4485527
                                                                      |                                                                                    
aTTCACCACTTTATCGGCACCCATCGCCAGCACTTCCTGCTCTTTCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTg                                                                   >  1:76733/1‑90 (MQ=255)
               ggCACCCATCGCCAGCACTTCCTGCTCTTTCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGc                                                    <  1:85773/90‑1 (MQ=255)
                  aCCCATCGCCAGCACTTCCTGCTCTTTCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAg                                                     >  1:209650/1‑86 (MQ=255)
                  aCCCATCGCCAGCACTTCCTGCTCTTTCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAg                                                     <  2:209650/86‑1 (MQ=255)
                    ccATCGCCAGCACTTCCTGCTCTTTCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGatat                                               <  2:359967/90‑1 (MQ=255)
                            aGCCCTCCCTGCTCTTCCGCCGAATTAAAACAAAAGGCTCCCCCCCCCCACCCCATTGCGGGCAGAATTTTTATAGCGATATCCCCCAGc                                       <  1:369661/90‑1 (MQ=255)
                            aGCACTTCCTGCTCTTTCGCCGGATTAGAACTAAAGGCTGTCCCCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGc                                       <  2:144678/90‑1 (MQ=255)
                                 ttCCTGCTCTTTCGCCGGATTAGAACTAAAGGCTGTCCCCTCGCACCCCATTGCGGGCAGAAGTTTTATAGCGATATGCCCCAGCCCGcc                                  <  2:421245/90‑1 (MQ=255)
                                 ttCCTGCTCTTCCCCCGGATTAGAACTAAAGGCTGTCCCCCCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGcc                                  <  2:479177/90‑1 (MQ=255)
                                     tGCTCTTTCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATa                              <  2:227996/90‑1 (MQ=255)
                                       ctctTTCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATAcc                            >  1:466294/1‑90 (MQ=255)
                                          tttCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAAt                         >  2:206220/1‑90 (MQ=255)
                                           ttCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAAtt                        <  2:191565/90‑1 (MQ=255)
                                                   aTTAGAACTAAAGGCTGTCCCCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAAc                <  2:303748/90‑1 (MQ=255)
                                                     tAGAACTAAAGGCTGTCCCCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAAcgc              <  2:17842/90‑1 (MQ=255)
                                                           tAAAGGCTGTCCCCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTgg        <  2:141218/90‑1 (MQ=255)
                                                            aaaGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGt       >  2:398962/1‑90 (MQ=255)
                                                             aaGGCTGTCCCCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTa      <  2:115306/90‑1 (MQ=255)
                                                             aaGCCTGTCCCCTCGCATCCCATTGCTTCCAGAATTTTTATAGCGATATGCCCCGCCCCGCCAATACCAATTACCCCAACGCGGCTGGTa      <  2:92606/90‑1 (MQ=255)
                                                               ggCTGTCCCCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGc    <  2:76733/90‑1 (MQ=255)
                                                               ggCTGTCCCCCCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGc    <  1:202221/90‑1 (MQ=255)
                                                               ggCTGTCCCCCCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCACCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGc    <  2:470754/90‑1 (MQ=255)
                                                               ggCTGCCCCCTCGCACCCCTTTGCGTGCAGAAGTTTTATGGCGATATCCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGc    <  2:246140/90‑1 (MQ=255)
                                                                 cTGTCCCCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAg  <  2:59144/90‑1 (MQ=255)
                                                                 cTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAg  >  1:13746/1‑90 (MQ=255)
                                                                      |                                                                                    
ATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCTGCTCTTTCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAG  >  NZ_CP009273/4485373‑4485527

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