Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913310,9340AC43.5% 23.6 / 11.4 23F135V (TTC→GTC) ecpRDNA‑binding transcriptional dual regulator MatA
Reads supporting (aligned to +/- strand):  ref base A (9/4);  new base C (0/10);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 1.60e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTAT  >  NC_000913/310870‑310992
                                                                |                                                          
gTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGaa                                                            <  1:1018232/65‑1 (MQ=255)
    tGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTg                                                        >  1:423123/1‑65 (MQ=255)
      cGGATAATTCCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGcc                                                      <  1:1188198/65‑1 (MQ=255)
       ggATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTgtct                                                     >  1:5008/1‑62 (MQ=255)
       ggATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCa                                                     >  1:289529/1‑65 (MQ=255)
       ggATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCa                                                     >  1:3967027/1‑65 (MQ=255)
       ggATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCa                                                     >  1:1156329/1‑65 (MQ=255)
          tAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTa                                                  <  1:4021664/65‑1 (MQ=255)
               tccatttccCGCGCCGTGATTTTATGCTTCTTAATCTCTTTTCTCAGGAGCCTGCCATTAATcac                                             <  1:1743804/65‑1 (MQ=255)
               tccattcccCGGTCCGTGATTTTATCTTTCTAAATATCTTTTCTCAGGACCCTGCCATTAATcac                                             <  1:1869684/65‑1 (MQ=255)
                    ttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                        <  1:1779567/65‑1 (MQ=255)
                    ttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGa                                        <  1:4652820/65‑1 (MQ=255)
                    ttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGa                                        <  1:188188/65‑1 (MQ=255)
                    ttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGAAACCTGCCATTAATCACATAGa                                        <  1:4554964/65‑1 (MQ=255)
                    tccccGTTCCGTGATTTTATCTTTCTTAATTTCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                        <  1:1648001/63‑1 (MQ=255)
                    tccccGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCTATTAATCACATAGa                                        <  1:3679060/63‑1 (MQ=255)
                    tccccGGCCCGTGATTTTATCTTTCTTAATATCTTTGCTCAGGAACCTGCCATTAATCACATAGa                                        <  1:346835/63‑1 (MQ=255)
                       ccGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCa                                     >  1:4283119/1‑65 (MQ=255)
                         ggTCCGTGTTTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGt                                   <  1:1138987/65‑1 (MQ=255)
                         ggTCCGTGATTTTATCTTTCTAAATATCTTTTCACAGGACCCTGCCATTAATCACATAGACCAGt                                   <  1:1096623/65‑1 (MQ=255)
                            ccGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCt                                <  1:3560610/65‑1 (MQ=255)
                                     tATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAata                       <  1:1630834/65‑1 (MQ=255)
                                     tATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAata                       >  1:140396/1‑65 (MQ=255)
                                     tATCTTTCTTAATATCTTTTCTAAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAata                       <  1:366895/65‑1 (MQ=255)
                                       tctttcttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATc                     <  1:2969379/65‑1 (MQ=255)
                                       tctttcttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATc                     <  1:1557392/65‑1 (MQ=255)
                                       tctttcttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATc                     <  1:200584/65‑1 (MQ=255)
                                          ttcttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACg                  <  1:3884134/65‑1 (MQ=255)
                                            cttAATGTCTTTTCTCCGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                <  1:3392745/65‑1 (MQ=255)
                                            cttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                <  1:4499607/65‑1 (MQ=255)
                                            cttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                <  1:3080430/65‑1 (MQ=255)
                                                atttCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAg            <  1:2329401/62‑1 (MQ=255)
                                                atatCTTTTCGCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAg            <  1:50710/65‑1 (MQ=255)
                                                 tttCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:1930835/63‑1 (MQ=255)
                                                 tttCTTTTCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:1730140/63‑1 (MQ=255)
                                                 tttCTTTTCCCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:2482791/63‑1 (MQ=255)
                                                 tctCTTTTCCCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:3527805/63‑1 (MQ=255)
                                                 tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:4080812/65‑1 (MQ=255)
                                                 tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:1750258/65‑1 (MQ=255)
                                                 tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:4006914/65‑1 (MQ=255)
                                                 tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:4070518/65‑1 (MQ=255)
                                                 tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:43817/65‑1 (MQ=255)
                                                 tatCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:3998137/65‑1 (MQ=255)
                                                 tatCTTTTCCCGGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:3329331/65‑1 (MQ=255)
                                                 tatCTTTTCCCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:4594841/65‑1 (MQ=255)
                                                 tatCTTTTCCCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:2273280/65‑1 (MQ=255)
                                                 tatCTTTTCCAAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt           <  1:3175042/65‑1 (MQ=255)
                                                   tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTcc         >  1:3747691/1‑65 (MQ=255)
                                                   tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTcc         >  1:695526/1‑65 (MQ=255)
                                                     ttttCTCGGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg       <  1:1368535/65‑1 (MQ=255)
                                                     ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg       <  1:4231684/65‑1 (MQ=255)
                                                     ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg       <  1:1222681/65‑1 (MQ=255)
                                                     ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg       <  1:391627/65‑1 (MQ=255)
                                                     ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg       <  1:3494945/65‑1 (MQ=255)
                                                     ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg       <  1:3483571/65‑1 (MQ=255)
                                                     ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg       <  1:646643/65‑1 (MQ=255)
                                                     ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg       <  1:2058945/65‑1 (MQ=255)
                                                     ttttCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg       <  1:3073864/65‑1 (MQ=255)
                                                     tttgCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg       <  1:2549942/65‑1 (MQ=255)
                                                      tttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc      <  1:787433/65‑1 (MQ=255)
                                                      tttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc      <  1:739266/65‑1 (MQ=255)
                                                      tttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc      <  1:1411903/65‑1 (MQ=255)
                                                      tttCACAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc      <  1:329413/65‑1 (MQ=255)
                                                       ttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCg     <  1:2743199/65‑1 (MQ=255)
                                                        cctcAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGt    <  1:3326576/64‑1 (MQ=255)
                                                          tcAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGtat  >  1:2199684/1‑65 (MQ=255)
                                                                |                                                          
GTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTAT  >  NC_000913/310870‑310992

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