Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,499,7630AC65.4% 71.6 / 58.0 52V192G (GTT→GGT) ydcOBenE family transporter YdcO
Reads supporting (aligned to +/- strand):  ref base A (10/8);  new base C (18/16);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGCTTTCATTGCTGCGATACCCGGTGCGTTTTGCGATGCCATCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAA  >  NZ_CP009273/1499628‑1499890
                                                                                                                                       |                                                                                                                               
gctgctTTCATTGCTGCGATACCCGGTGCGTTTTGCGATGCCATCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCAt                                                                                                                              <  1:476085/139‑1 (MQ=255)
         aTTGCTGCGATACCCGGTGCGTTTTGCGATGCCATCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGa                                                                                                                     <  2:215566/139‑1 (MQ=255)
                        ggTGCGTTTTGCGATGCCATCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTgcg                                                                                                      >  1:513186/1‑139 (MQ=255)
                        ggTGCGTTTTGCGATGCCATCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTgcg                                                                                                      <  2:726796/139‑1 (MQ=255)
                           gCGTTTTGCGATGCCATCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGa                                                                                                   >  2:580594/1‑139 (MQ=255)
                               tttGCGATGCCATCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGac                                                                                               >  2:1210663/1‑139 (MQ=255)
                                        ccATCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcacga                                                                                       <  1:679023/138‑1 (MQ=255)
                                        ccATCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcacga                                                                                       <  2:996258/138‑1 (MQ=255)
                                                     aaaaaGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatcatc                                                                         >  2:1273219/1‑139 (MQ=255)
                                                          gggggAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGTGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGc                                                                    <  1:580594/139‑1 (MQ=255)
                                                          gggggAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCgg                                                                     >  1:196232/1‑138 (MQ=255)
                                                          gggggAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGc                                                                    <  1:1253979/139‑1 (MQ=255)
                                                            gggAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcacg                                                                                        >  2:451278/1‑117 (MQ=255)
                                                            gggAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGacgatcacg                                                                                        <  1:451278/117‑1 (MQ=255)
                                                               aGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTa                                                               >  2:20743/1‑139 (MQ=255)
                                                                 tGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTAcc                                                             >  2:14914/1‑139 (MQ=255)
                                                                      cGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCa                                                          >  2:87420/1‑137 (MQ=255)
                                                                        cTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAgc                                                      >  2:972572/1‑139 (MQ=255)
                                                                        cTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAgc                                                      >  1:1017350/1‑139 (MQ=255)
                                                                           aGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGcgcg                                                   >  1:1127595/1‑139 (MQ=255)
                                                                           aGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGcgcg                                                   >  2:325386/1‑139 (MQ=255)
                                                                           aGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGcgcg                                                   >  2:1033899/1‑139 (MQ=255)
                                                                                 tgAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCa                                              <  2:1006468/138‑1 (MQ=255)
                                                                                   aGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACg                                           >  1:1019664/1‑139 (MQ=255)
                                                                                       aaCGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATcgc                                                                       <  1:891137/107‑1 (MQ=255)
                                                                                       aaCGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATcgc                                                                       >  2:891137/1‑107 (MQ=255)
                                                                                       aaCGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATACCAATAATCATCGCGGAAATTACCGCATAGCGCGGCGCAACGGCCt                                       >  2:443228/1‑139 (MQ=255)
                                                                                             aaaTCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTg                                 <  1:657808/139‑1 (MQ=255)
                                                                                              aaTCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgc                                >  1:145953/1‑139 (MQ=255)
                                                                                                       tAATATAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATAc                       <  2:513186/139‑1 (MQ=255)
                                                                                                       tAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTg                                 <  1:658149/129‑1 (MQ=255)
                                                                                                       tAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTg                                 >  2:658149/1‑129 (MQ=255)
                                                                                                         atatAAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATAc                       >  1:436400/1‑137 (MQ=255)
                                                                                                              aGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGa                                                                                                >  2:211762/1‑59 (MQ=255)
                                                                                                              aGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCaa                 >  1:1277297/1‑138 (MQ=255)
                                                                                                              aGTGAGGAGAACGGGTTTAAAGACGCCATCAGTTGTGACAGCGTCACCTTGCGCGATGa                                                                                                <  1:211762/59‑1 (MQ=255)
                                                                                                               gTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACa               <  1:14914/139‑1 (MQ=255)
                                                                                                               gTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACa               <  2:1260241/139‑1 (MQ=255)
                                                                                                               gTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGATGCCAGCCATACCAGCAACa               >  1:1260241/1‑139 (MQ=255)
                                                                                                               gTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACa               <  1:318239/139‑1 (MQ=255)
                                                                                                               gTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACa               <  2:145953/139‑1 (MQ=255)
                                                                                                               gTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACa               <  1:20743/139‑1 (MQ=255)
                                                                                                               gTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACa               <  1:1273219/139‑1 (MQ=255)
                                                                                                                  gggAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATa             <  1:87420/138‑1 (MQ=255)
                                                                                                                  gggAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATAc            >  1:59748/1‑139 (MQ=255)
                                                                                                                   ggAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGc                                                                    <  2:398591/82‑1 (MQ=255)
                                                                                                                   ggAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGc                                                                    >  1:398591/1‑82 (MQ=255)
                                                                                                                       aaCGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCca       >  2:1078483/1‑139 (MQ=255)
                                                                                                                          gggTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCcacac    <  1:1078483/139‑1 (MQ=255)
                                                                                                                            gTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCt                                       <  2:303631/102‑1 (MQ=255)
                                                                                                                            gTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCt                                       >  1:303631/1‑102 (MQ=255)
                                                                                                                            gTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACaa  >  2:1080988/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                               
GCTGCTTTCATTGCTGCGATACCCGGTGCGTTTTGCGATGCCATCGTCACCAGAAAAAGGGGGAGTGCAACGCTCAGGCTGTGAGCAAACGAAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAA  >  NZ_CP009273/1499628‑1499890

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑