Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9630TG55.6% 9.5 / 42.6 36D492A (GAC→GCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/15);  new base G (20/0);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 2.87e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA  >  NZ_CP009273/822886‑823045
                                                                             |                                                                                    
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCGACGCCGGGAg                                                                         >  2:288438/1‑90 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGGAg                                                                         >  1:433664/1‑90 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGGAg                                                                         >  1:42558/1‑90 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGAGGGGGGGTGAAAGCCGGagg                                                                         >  1:1355971/1‑88 (MQ=255)
    ggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGGAAgcgg                                                                      >  2:1860307/1‑87 (MQ=255)
       gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGGAGTCGGGtc                                                                  >  1:683438/1‑90 (MQ=255)
            tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGGCGGAAGGCGGGtcgtcgg                                                             >  2:1502261/1‑89 (MQ=255)
               ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTGGGGtcgtcgggga                                                          >  2:558172/1‑86 (MQ=255)
                 gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAc                                                        >  1:676141/1‑90 (MQ=255)
                     cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa                                                    <  1:1239524/90‑1 (MQ=255)
                         tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACCCCCGGAGGCGGGTCGTCGAGAAAAAGAAAGTc                                               >  2:1254695/1‑90 (MQ=255)
                          ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                              <  2:1515725/90‑1 (MQ=255)
                          ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                              <  2:1376956/90‑1 (MQ=255)
                          ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                              <  1:1251657/90‑1 (MQ=255)
                                 ccATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGTTGTCGGGAAAAAGAAAGTCCGGTTCAt                                       >  2:659908/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                     <  1:659908/90‑1 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                     <  2:368665/90‑1 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                     <  2:493857/90‑1 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGTCGGTTCGGCGGGAAAAAGAAAGTCCGGTTCATGc                                     >  1:233184/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGACGGGTCGTCCGGAAACAGAAAGTCCGGTTCATGc                                     >  1:430838/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCCGGAGGCGGGTTGTCGGGAAAAAAAATGTTCCGGttgtgc                                      >  2:620456/1‑86 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAACAGAAAGTCCGGTTCATGc                                     >  1:435827/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACCCCGGAAGTCGGGTCGGCGAGAAAAAGAATGTCCGGGTCATGc                                     >  1:368665/1‑90 (MQ=255)
                                      cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACCCCGGGAGGCGGGTCGTCGAGAAACAAAAAGTCCGGGTCATTCATc                                  >  1:1802031/1‑90 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGTTCGTCGGGAAAAAAAAAGTCCGGGTCATGCATCAg                                >  1:1578486/1‑90 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                           <  1:1182507/90‑1 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAAc                           >  1:1534743/1‑90 (MQ=255)
                                                gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTGGGGTGGTCGAGAAAAAGAATGTCCGGTTCATGCAGCAGCGGACAgg                        >  1:1095854/1‑90 (MQ=255)
                                                gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACCCCGGGAGGCGGGTCGTCGAGAAAAAAAAAGGCCGGGTCGTGCATCAGCGaaaaag                        >  2:1817480/1‑86 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATTTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc          <  2:1101228/90‑1 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc          <  2:209586/90‑1 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc          <  2:435827/90‑1 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc          <  1:29053/90‑1 (MQ=255)
                                                               ggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCt         <  2:1534743/90‑1 (MQ=255)
                                                                      gAGGGGGTCACCCCCGGAAGCCGGTCCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCCAACAGGCCAGCGCCCGCCGCTGTTTaaa  <  2:1033632/90‑1 (MQ=255)
                                                                      gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa  <  2:1443741/90‑1 (MQ=255)
                                                                             |                                                                                    
GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA  >  NZ_CP009273/822886‑823045

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