Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0129672,016,0910TA29.7% 66.5 / 16.6 37intergenic (‑123/+56)yeeE/plaPthiosulfate utilization transporter TsuA/YeeE/putrescine/proton symporter PlaP
Reads supporting (aligned to +/- strand):  ref base T (12/14);  new base A (0/11);  total (12/25)
Fisher's exact test for biased strand distribution p-value = 6.62e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTT  >  NC_012967/2016020‑2016160
                                                                       |                                                                        
aaTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAAt                                                                          <  1:1708575/72‑1 (MQ=255)
aaTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAAt                                                                          <  1:306267/72‑1 (MQ=255)
aaTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAAt                                                                          >  1:4498267/1‑72 (MQ=255)
 aTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATa                                                                         >  2:2077843/1‑72 (MQ=255)
   gCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaa                                                                      <  2:3735028/72‑1 (MQ=255)
      gACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaaaa                                                                   <  1:1595046/72‑1 (MQ=255)
      gACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaaaa                                                                   >  2:4961128/1‑72 (MQ=255)
      gACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaaaa                                                                   <  1:5067151/72‑1 (MQ=255)
      gACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaaaa                                                                   >  1:4494557/1‑72 (MQ=255)
       aCCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaaaaa                                                                  >  2:7693901/1‑72 (MQ=255)
         cAATANCCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaacaaaac                                                                >  1:862618/1‑70 (MQ=255)
         cAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAAcc                                                               <  1:4503112/72‑1 (MQ=255)
          aaTAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaaaaaaca                                                               >  2:5390804/1‑69 (MQ=255)
           aTAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAAcacc                                                             >  2:7741707/1‑69 (MQ=255)
            tAACCAAATTCTCTAAATTAGAACAAAGGGTTATTAATGAGGATTTCGATAAAAAAGAAAAAAAAAACCGGa                                                            <  2:5844577/72‑1 (MQ=255)
             aaCAAAATTCTCTAAATAAGAACAAATGTTTATTAAGGAGAATTTCGAAAAAAAAAAAAAAAAAAACCGGAg                                                           <  1:3202923/72‑1 (MQ=255)
              aCCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGc                                                          <  1:3963554/72‑1 (MQ=255)
                aaaaTTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCa                                                        <  2:6282353/71‑1 (MQ=255)
                 aaaTTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGACTACAAAAAACCGGAGcac                                                        >  2:6549912/1‑70 (MQ=255)
                  aaTTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCCGGGGCCaa                                                       >  2:6300591/1‑72 (MQ=255)
                  aaTTCTATAAATTAGAAAAAATGTTTATTAATGAGGATTTCGAAAAAAAAGAAAAAAAAAACCGGAGCCaaa                                                      <  2:4059026/72‑1 (MQ=255)
                   aTTCTCAAAATAAGAAAAAATGTTTATAAAGGAGAATTTCGAAAAAAAAAAAAAAAAAAACCGGAGCCAAat                                                     <  1:8136016/72‑1 (MQ=255)
                    ttCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCCGGGGCaaaaa                                                    >  2:8463824/1‑70 (MQ=255)
                       tctAAATTAGAAAAAATGTTTATTAATGAGGATTTCAAAAAAAAAAAAAAAAAAAACCGGAGCCAAATattc                                                 <  1:4028703/72‑1 (MQ=255)
                          aaaTTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATattcatt                                              >  2:6992121/1‑72 (MQ=255)
                          aaaTTAGAAAAAATGGTTATTAATGAGGATTTCGAAAAAAAAAAAAAAAAAAACCGGAGCCAAATattcatt                                              <  2:2135817/72‑1 (MQ=255)
                           aaTTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATattcattc                                             <  1:4906536/72‑1 (MQ=255)
                           aaTTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATattcattc                                             <  2:7314835/72‑1 (MQ=255)
                                 aaaaaaTGGTTATTAATGAGGATTTCGATAAAAAAAAAAAAAAAAACCGGAGCCAAATATTCATTCGACTcc                                      <  2:6290352/69‑1 (MQ=255)
                                  aaaaaTGTTAATAAAGGAGGATTTCAAAAAAAAAAAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCg                                     <  1:6429526/70‑1 (MQ=255)
                                    aaaTGTTTATTAATGAGGATTTCAATAAAAAAGAAAAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGt                                   <  2:6176762/72‑1 (MQ=255)
                                     aagggTTATTAAGGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGtt                                  <  2:7316353/69‑1 (MQ=255)
                                      aTGGTTATTAAGAAGGATTTCAATAAAAAAGAAGAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTg                                 <  2:6979821/72‑1 (MQ=255)
                                        ggTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGat                               <  2:9213087/72‑1 (MQ=255)
                                         gTTATTAATAAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGatt                              <  2:6742938/72‑1 (MQ=255)
                                         gTTATTAAAGAAGATTTTAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGat                               <  1:4040972/72‑1 (MQ=255)
                                            ataaaTGAGGATTTCAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGattatt                           <  1:4173168/69‑1 (MQ=255)
                                             ttaATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTGCGGTTGAttatta                          <  1:6933275/72‑1 (MQ=255)
                                               aaTGAGGATTTCGATAAAAAAAAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGa                        <  2:2401633/72‑1 (MQ=255)
                                                    ggATTTCGATAAAAAAGAATAAAAAAACCCGGGGCCAAAAATCTTTTTACTCCCGGTTGGTTTTTTGaaaaa                    >  1:8928496/1‑71 (MQ=255)
                                                    ggATTTAGAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAg                    <  1:1643288/72‑1 (MQ=255)
                                                           aaaaaaaaaaaaaaaaaaaCCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTAt            <  1:8076177/69‑1 (MQ=255)
                                                             aaaaaaaaaaaaaaaaaCCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGc          <  2:8274469/71‑1 (MQ=255)
                                                                aaaaaaaaaaaaaaCCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTAc       <  1:395662/72‑1 (MQ=255)
                                                                aaaaGAAAAAAAAACCCGGAGCNAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTAc       <  1:2562362/72‑1 (MQ=255)
                                                                 aaaGAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACg      <  2:8578120/72‑1 (MQ=255)
                                                                   aGAAAAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTc    <  1:4589930/72‑1 (MQ=255)
                                                                     aaTAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCtt  <  1:2423118/72‑1 (MQ=255)
                                                                       |                                                                        
AATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAACCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTT  >  NC_012967/2016020‑2016160

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