Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,419,8010TG56.5% 11.5 / 26.4 23L41V (TTG→GTG) cycAD‑serine/D‑alanine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base T (6/4);  new base G (13/0);  total (19/4)
Fisher's exact test for biased strand distribution p-value = 2.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-04
Rejected as consensus: Frequency below/above cutoff threshold.

ATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTC  >  NZ_CP009273/4419712‑4419887
                                                                                         |                                                                                      
atCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt                                                                                        <  2:286376/90‑1 (MQ=255)
atCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt                                                                                        <  1:262927/90‑1 (MQ=255)
atCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt                                                                                        <  1:653064/90‑1 (MQ=255)
 tCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGGCCCCTTGGTCAGGGgtt                                                                                       >  1:37255/1‑90 (MQ=255)
      cTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGGACGGGGTGGtttt                                                                                  >  2:338852/1‑89 (MQ=255)
        ccGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTg                                                                                >  2:170023/1‑90 (MQ=255)
         cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTgg                                                                               >  2:233450/1‑90 (MQ=255)
         cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGGGCCCTTGGGACGGGGGGGTTTTTgg                                                                               >  1:341955/1‑90 (MQ=255)
         cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGGGCCATTGGGACGGGGGTGTTTATgg                                                                               >  2:653064/1‑90 (MQ=255)
                        tACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGGCCTTTGGACGGGGGTGTTTTTGGGGGCTGGCAAAAAGa                                                                >  1:792432/1‑90 (MQ=255)
                             cgcAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGc                                                           <  1:741116/90‑1 (MQ=255)
                                 aTCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTGGGGGCTTGCAAAAAGATTAGCCTTg                                                       >  1:76720/1‑90 (MQ=255)
                                 aTCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGGACGGGGGTGTTTATGGGGGCTGGCAAAACGATTAGCCTTg                                                       >  2:155745/1‑90 (MQ=255)
                                      acaAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCggg                                                  >  2:67041/1‑90 (MQ=255)
                                                     cAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGGCTGGCAAAAAGATTAGCCTTTCCGGGGCGGCGGTCAttttt                                   >  1:249681/1‑89 (MQ=255)
                                                       gCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTTGGGGCTGGCCAAAAGATTAGCCTTGCCGGGCCGGCGGTCCTTTTTGt                                 >  1:546977/1‑90 (MQ=255)
                                                          tATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTta                              >  1:280117/1‑90 (MQ=255)
                                                            ttGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTtata                            >  2:267108/1‑90 (MQ=255)
                                                            ttGCCATTGGCGGTGCCATTGGGACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTtata                            >  2:164469/1‑90 (MQ=255)
                                                                 aTTGGCGGTGCCATTGGTACGGGGGTGGTTTTGGGGTCTGGCAAAACGATTAGCCTTGCCGGGGCGGCGGTCCTTTTTGTTTATATGGTc                       >  1:643531/1‑90 (MQ=255)
                                                                    ggCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCAtt                    <  2:365338/90‑1 (MQ=255)
                                                                    ggCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCAtt                    <  1:170023/90‑1 (MQ=255)
                                                                                  tACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCtt      >  2:261698/1‑90 (MQ=255)
                                                                                  tACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCtt      >  1:163785/1‑90 (MQ=255)
                                                                                  tACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGACCATTGGTTTTATGCCCtt      <  2:38955/90‑1 (MQ=255)
                                                                                     ggggTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCttttt   >  2:503015/1‑90 (MQ=255)
                                                                                      gggTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTc  >  1:177182/1‑90 (MQ=255)
                                                                                         |                                                                                      
ATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTTTC  >  NZ_CP009273/4419712‑4419887

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