Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,751,0280CT40.5% 37.2 / 19.0 37A54T (GCA→ACA) yiaYL‑threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (15/7);  new base T (0/15);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 2.08e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGC  >  NZ_CP009273/3750942‑3751109
                                                                                      |                                                                                 
tttCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCttt                                                                                >  2:70682/1‑90 (MQ=255)
tttCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCttt                                                                                >  1:805/1‑90 (MQ=255)
      aCCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGcac                                                                          <  2:306273/90‑1 (MQ=255)
         tGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGcac                                                                          <  2:943561/87‑1 (MQ=255)
         tGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGcac                                                                          >  1:943561/1‑87 (MQ=255)
          gcggcgACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATcg                                                                      >  1:173268/1‑90 (MQ=255)
          gcggcgACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATcg                                                                      >  1:90320/1‑90 (MQ=255)
          gcggcgACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATcg                                                                      >  1:260148/1‑90 (MQ=255)
                          gggggggtTTAGGTGGGGTCCCATAAAAAAAAACCCAAAAAATTTTGCTTTCTTCCAGTTTTTTTTCCACCCCGCCCGCCATACCTAAtt                                                      <  2:731311/86‑1 (MQ=255)
                           ggtggGTTTAGGTTGGGTCCCATCAAAAAAAACGCAAAAAATTTTGCTTTCTTCCAGTGTTTTTTGCACACCGCCCGCCATACCTAAttt                                                     <  2:805/90‑1 (MQ=255)
                            gtggGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTc                                                    >  2:806489/1‑90 (MQ=255)
                                ggTTGGGTGGGGTCCCATAAAAAAAAACGCAAAAAATTTTGCGTTCTCCCAGTGTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTa                                                <  2:164180/90‑1 (MQ=255)
                                ggTTAGGTTGGGTGCCATAAAAAAAAACGAAAAAAATTTTGCGTTCTTCCAGTGTTTTTTGCCCACCGCCCGCCATACCTAATTTCGTTa                                                <  2:90469/90‑1 (MQ=255)
                                ggTTAGGTTGGGTCCCATAAAAAATACCAAAAAAAATTTTGCTTTCTTCCAGTTTTTTTTCCACATCGCCCGCCATACCTAATTTCGTTa                                                <  2:180624/90‑1 (MQ=255)
                                ggTTAGGTGGGGTGCCATAAAAAAAAACCCAAAAAATTTTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTa                                                <  1:916650/90‑1 (MQ=255)
                                         gggTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGt                                       >  1:746305/1‑90 (MQ=255)
                                              ccATCAAAAATAACGCAAAAAATATTGCTTTTTTCCAGTGTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTg                                  <  1:895415/90‑1 (MQ=255)
                                              ccATCAAAAAAAACGAAAAAAATATTGCTTTCTTCCAGTCTTTTTTGCACATCCCCCGCCATACCTAATTTCGTTAACATATTGTCAGTg                                  <  2:260148/90‑1 (MQ=255)
                                              ccATAAAAAAAAACGCAAAAAATTTTGCTTTCTTCCAGTGTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTg                                  <  1:489645/90‑1 (MQ=255)
                                               catcatAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCa                                                               <  1:795660/60‑1 (MQ=255)
                                               catcatAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCa                                                               >  2:795660/1‑60 (MQ=255)
                                               catcatAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGa                                 >  2:741542/1‑90 (MQ=255)
                                                 tcaaAAATAACGCAAAAAATATTGCTTTCTTCCAGTGTTTTTTCCACACCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACa                               <  1:806489/90‑1 (MQ=255)
                                                    aaaataaCGCTAAAAATTTTGCTTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAAtt                            <  1:807289/89‑1 (MQ=255)
                                                     aaataaCGCAAAAAATATTGCTTTCTTCCAGTGTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTa                           <  1:200882/90‑1 (MQ=255)
                                                      aatacCGCAAAAAATTTTGCTTTCTCCCGTTTTTTTTTCCACATCCCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTaa                          <  2:13001/90‑1 (MQ=255)
                                                          aCGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATcgcccgcc                                                                >  2:7926/1‑48 (MQ=255)
                                                          aCGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATcgcccgcc                                                                <  1:7926/48‑1 (MQ=255)
                                                          aCGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTa                      >  1:232121/1‑90 (MQ=255)
                                                          aCGCAAAAAATATTGCTTTTTTCCAGTGTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTa                      <  2:232121/90‑1 (MQ=255)
                                                           cacaaaaaaTTTTGCTTTCTTCCAGTTTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTAc                     <  2:90320/86‑1 (MQ=255)
                                                                     atTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTa                >  2:118337/1‑85 (MQ=255)
                                                                     atTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTa                <  1:118337/85‑1 (MQ=255)
                                                                       tGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCAt         >  1:692496/1‑90 (MQ=255)
                                                                             cttcCAGTGTTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTa                <  2:986980/77‑1 (MQ=255)
                                                                             cttcCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTa                >  1:986980/1‑77 (MQ=255)
                                                                              ttcCAGTGTTTTTTCCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGc  <  1:261711/90‑1 (MQ=255)
                                                                                      |                                                                                 
TTTCAAACCTGCGGCGACGTTTTCCGTGGTGGGGTTAGGTTGGGTGCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGC  >  NZ_CP009273/3750942‑3751109

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