Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0129672,016,0910TA46.9% 12.2 / 26.9 32intergenic (‑123/+56)yeeE/plaPthiosulfate utilization transporter TsuA/YeeE/putrescine/proton symporter PlaP
Reads supporting (aligned to +/- strand):  ref base T (6/11);  new base A (0/15);  total (6/26)
Fisher's exact test for biased strand distribution p-value = 1.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTT  >  NC_012967/2016020‑2016160
                                                                       |                                                                      
aaTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAAt                                                                        >  1:3178035/1‑72 (MQ=255)
   gCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaa                                                                     >  1:6752106/1‑72 (MQ=255)
    cAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaa                                                                    >  2:247510/1‑72 (MQ=255)
    cAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaa                                                                    >  1:7716824/1‑72 (MQ=255)
     aGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaaa                                                                   >  1:5473463/1‑72 (MQ=255)
          aaTAACCAAATTCTCTAAATTAGAACAAATTGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCg                                                             <  2:9476223/72‑1 (MQ=255)
                 aaaTTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCaa                                                      <  2:7792694/72‑1 (MQ=255)
                    ttCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAata                                                   <  1:2408627/72‑1 (MQ=255)
                       tctAAATTAGAAAAAATGTTTATTAATGAGGATTTCGAAAAAAAAGAATAAAAAAACCGGAGCCAAATattc                                                <  1:5315050/72‑1 (MQ=255)
                       tcaaaaTAAGAAAAAATGGTAATTAATGAGGATTTCGATAAAAAAAAATAAAAAAACCGGAGCCAAATattc                                                <  2:330811/69‑1 (MQ=255)
                       ctaaaaTTAAAAAAAATGGTTTTTAAAGGGGATTTCAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATatt                                                 <  1:6156457/69‑1 (MQ=255)
                        taaaaTTAAAAAAAAGTGTTATTAAAGGGGATTTAAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATattc                                                <  1:3447722/70‑1 (MQ=255)
                          aaaTTAGAAAAAATGGTTTTTAATGAGGTTTTCGAAAAAAAAGAAAAAAAAAACCGGAGCCAAATattcatt                                             <  2:7390903/72‑1 (MQ=255)
                            aTTAGAAAAAGTGGTTATTAATGAGGATTACGATAAAAAAGAAAAAAAAAGCCGGAGCCAAATATTCATTCg                                           <  2:8157229/72‑1 (MQ=255)
                              tAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGAc                                         <  2:7769801/72‑1 (MQ=255)
                              tAGAAAAAATGGTTATTAATGAGGATTTCAATAAAAAAAAAAAAAAAAACCGGAGCCAAATATTCATTCGAc                                         <  1:5664925/72‑1 (MQ=255)
                                 aaaaaaGGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTcc                                      <  2:7832947/69‑1 (MQ=255)
                                  aaaaaTGGTTTATAAAGGGGGTTTCGAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCTACTcc                                      <  1:6779932/70‑1 (MQ=255)
                                   aaaaTGGTTTTTAAGGGGGTTTTGAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCg                                     <  2:3208190/71‑1 (MQ=255)
                                    aaaTGGTTTATAAAGGGGATTTCAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCgg                                    <  1:8488087/72‑1 (MQ=255)
                                    aaaTGGTTATTAATGGGAATTTCGATAAAAAAGAAAAAAAAAGCCGGAGCCAAATATTCATTCGACTCCGGt                                   <  2:7190886/72‑1 (MQ=255)
                                    aaaTGGTTATTAAAGAGGATTTAAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCgg                                    <  1:8281487/72‑1 (MQ=255)
                                     aaTGGTTATTAATAGTGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGtt                                  <  2:9845289/72‑1 (MQ=255)
                                      aTGGTTATTAATGAGGATTTCAATAAAAAAAAAAAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTg                                 <  1:9193961/72‑1 (MQ=255)
                                        gtttataaAGGAGGATTTCGATAAAAAAAAATAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGat                               <  2:280828/70‑1 (MQ=255)
                                        gtttataaAGGAGGATTCGGAAAAAAAAAAAAAAAAAAGCCGGAGCCAAATATTCATTCGACTCAGGTTGat                               <  1:5501469/70‑1 (MQ=255)
                                        ggTTATGAATGAGGATTTCGATAAAAAAGAAAAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGat                               <  1:1961538/72‑1 (MQ=255)
                                          ttattaATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGatta                             <  1:4616695/72‑1 (MQ=255)
                                          ttattaATGAGGATTTCAATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGatta                             <  2:8575834/72‑1 (MQ=255)
                                             taaaTGAGATTTCCAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGAttatta                          <  1:8467722/70‑1 (MQ=255)
                                              taATGAGGATTTCGATAAAAAAGAATAAAAAAACCCGATCCAAAAAATTATTTGACTCCCGTTGGTTATTAg                         >  1:2132961/1‑72 (MQ=255)
                                               aaTGAGGATTTCGATAAAAAAAAAAAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGa                        <  2:7628389/72‑1 (MQ=255)
                                                    gagtttCGAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAg                    <  2:4234711/69‑1 (MQ=255)
                                                     gATTTAGAAAAAAAAGAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAga                   <  2:3384463/72‑1 (MQ=255)
                                                        ttCGATAAAAAAGAATAAAAAAACCGGAGACAAATATTCATTCGACTCCGGTTGATTATTAGAAAAgatgat               >  2:9371540/1‑72 (MQ=255)
                                                           gaaaaaaaaaaaaaaaaaaCCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTAt            <  1:9424381/69‑1 (MQ=255)
                                                           aaaaaaaaaGAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTAt            <  1:2006600/69‑1 (MQ=255)
                                                            aaaaaaaaaaaaaaaaaaCCCGGGGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATg           <  2:3613189/70‑1 (MQ=255)
                                                            aaaaaaaaaaaaaaaaaaCCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATg           <  2:2580260/70‑1 (MQ=255)
                                                             gaaaaaaaaaaaaaaaaCCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGc          <  2:3462807/71‑1 (MQ=255)
                                                                aaaaaaaaaaaaaaACCGGAGCCAAATATTCATTCGANTCCGGTTGATTATTAGAAAAGATGATTATGCTAc       <  1:529526/72‑1 (MQ=255)
                                                                     aaTAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCtt  <  2:8908456/72‑1 (MQ=255)
                                                                       |                                                                      
AATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTT  >  NC_012967/2016020‑2016160

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: