Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,419,8010TG52.4% 18.8 / 28.3 21L41V (TTG→GTG) cycAD‑serine/D‑alanine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base T (1/9);  new base G (11/0);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 3.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-05
Rejected as consensus: Frequency below/above cutoff threshold.

TCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTT  >  NZ_CP009273/4419713‑4419875
                                                                                        |                                                                          
tCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgtt                                                                           <  1:161877/90‑1 (MQ=255)
tCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgtt                                                                           <  1:320784/90‑1 (MQ=255)
        cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGGCATTGGTACGGGGGTGTTTAttg                                                                   >  2:258300/1‑88 (MQ=255)
            tGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTc                                                               <  1:143289/90‑1 (MQ=255)
                        aCGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGAt                                                   <  2:57160/90‑1 (MQ=255)
                            cgcAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGGCTGGCAAAACGATTAGc                                               >  1:73579/1‑90 (MQ=255)
                                     acaAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCggg                                      >  1:87275/1‑90 (MQ=255)
                                              cATATTCAGCTTATTGCCATTGGCGGTGCCCTTTGGGCGGGGGTGTTTTTGGGGGCTGGCAAAAAGATTTTCCTTTCCGGGCCGTTGAtg                             >  1:255025/1‑89 (MQ=255)
                                              cATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTGGGGGCTGGCAAAACGATTAGCCTTTCCGGGCCGTCGATc                             >  2:148159/1‑90 (MQ=255)
                                                      gCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGt                     <  2:38581/90‑1 (MQ=255)
                                                       cTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGtt                    <  1:368277/90‑1 (MQ=255)
                                                           ttGCCATTGGCGGTGGCATTTGGACGGGGGTGTTTTTGGGGTCTGGCAAAAAGATTAGCCTTGCCGGGGCGTCGATCATTTTTGTTttta                >  2:295216/1‑90 (MQ=255)
                                                             gCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATg              <  1:118335/90‑1 (MQ=255)
                                                                aTTGGCGGTGCCATTGGGACGGGGGTGTTTTTGGGGGCTGGCAAAACGATTAGCCTTGCCGGGCCGGCGATCCTTTTTGTTTATATGATc           >  2:198526/1‑90 (MQ=255)
                                                                aTTGGCGGTGCCATTGGGACGGGGGTGTTTATGGGGGCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTTGTTTATATGATc           >  1:278491/1‑90 (MQ=255)
                                                                   ggCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCAtt        <  2:255025/90‑1 (MQ=255)
                                                                   ggCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCAtt        <  1:258300/90‑1 (MQ=255)
                                                                     cGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTgg      <  2:208653/90‑1 (MQ=255)
                                                                       gTGCCATTGGTACGGGGGTGTTTTTGGGGGCTGGCAAAACGATTAGCCTTGCCGGGGCGTCGGTCATTTTCGTTTTTATGATCATTGGtt    >  2:129806/1‑90 (MQ=255)
                                                                         gCCATTGGTACGGGGGTGTTTTTGGGGGCTGGCAAAACGATTAGCCTTGCCGGGCCGGCGGTCATTTTTCGTTATATGATCATTGGgttt  >  2:391803/1‑90 (MQ=255)
                                                                         gCCATTGGTACGGGGGTGTTTATGGGGGCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTTGTTTATATGATCATTGGtttt  >  2:116064/1‑90 (MQ=255)
                                                                                    ggggTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAtt                          <  2:149799/55‑1 (MQ=255)
                                                                                    ggggTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAtt                          >  1:149799/1‑55 (MQ=255)
                                                                                        |                                                                          
TCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTT  >  NZ_CP009273/4419713‑4419875

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