Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913310,9340AC52.0% 15.2 / 39.8 25F135V (TTC→GTC) ecpRDNA‑binding transcriptional dual regulator MatA
Reads supporting (aligned to +/- strand):  ref base A (11/1);  new base C (0/13);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 2.69e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-04
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACC  >  NC_000913/310873‑310995
                                                             |                                                             
aTGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCt                                                            >  1:1553510/1‑65 (MQ=255)
 tGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTg                                                           >  1:3852178/1‑65 (MQ=255)
    ggATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCa                                                        >  1:1543874/1‑65 (MQ=255)
      aTAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCAtt                                                      >  1:135434/1‑65 (MQ=255)
      aTAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCAtt                                                      >  1:912945/1‑65 (MQ=255)
      aTAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCAtt                                                      >  1:3720501/1‑65 (MQ=255)
            tccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATcac                                                >  1:2849482/1‑65 (MQ=255)
            tccattcccCGGTCCGTGTTTTTTTCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATcac                                                <  1:3622553/65‑1 (MQ=255)
            tccattcccCGGTCCGTGATTTTTTCTTTCTTAATATTTTTTCCCAGGACCCTGCCATTAATcac                                                <  1:3649642/65‑1 (MQ=255)
                 ttccCGGTCCGTGTTTTTTTCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGa                                           <  1:1905177/65‑1 (MQ=255)
                 ttccCGGTCCGTGTTTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                           <  1:2002497/65‑1 (MQ=255)
                 ttccCGGTCCGTGTTTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGa                                           <  1:1016968/65‑1 (MQ=255)
                 ttccCGGTCCGTGATTTTATTTTTCTTAATATCTTTTCCCAGGACCCTGCCATTAATCACATAGa                                           <  1:554402/65‑1 (MQ=255)
                 ttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                           <  1:3879905/65‑1 (MQ=255)
                 ttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                           <  1:195127/65‑1 (MQ=255)
                 ttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                           <  1:75573/65‑1 (MQ=255)
                 tccccGGCCCGTGATTTTATCTTTCTTAATATCTTTACCAAGAACCCTGCCATTAATCACATAGa                                           <  1:2452459/63‑1 (MQ=255)
                    ccGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCa                                        >  1:972720/1‑65 (MQ=255)
                     cGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAg                                       >  1:2443919/1‑65 (MQ=255)
                      ggTCCGTGTTTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGt                                      <  1:3323234/65‑1 (MQ=255)
                      ggTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGt                                      <  1:1230877/65‑1 (MQ=255)
                      ggTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGt                                      <  1:2353591/65‑1 (MQ=255)
                      cgTCCTTGTTTTTATTTTTCTTAATATTTTTCCGAAGAACCCTGCCATTAATCACATAGACCAGt                                      <  1:965385/64‑1 (MQ=255)
                         ccGTGTTTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCt                                   <  1:908702/65‑1 (MQ=255)
                              aTTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCt                              <  1:6560/65‑1 (MQ=255)
                              aTTTTATCTTTCTAAATATCTTTTCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCt                              <  1:3151505/65‑1 (MQ=255)
                                  tATCTTTCTTAATTTCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAata                          <  1:843768/65‑1 (MQ=255)
                                    tttttcttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATc                        <  1:460018/63‑1 (MQ=255)
                                    tctttcttAATTTCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATc                        <  1:2942182/65‑1 (MQ=255)
                                    tctttcttAATATTTTTTCCCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATc                        <  1:25094/65‑1 (MQ=255)
                                    tctttcttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATc                        <  1:2726533/65‑1 (MQ=255)
                                         tttAATATCTTTTCTCAGGCACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                   <  1:1592184/64‑1 (MQ=255)
                                         cttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                   <  1:1681101/65‑1 (MQ=255)
                                         cttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                   <  1:1632195/65‑1 (MQ=255)
                                         cttAAAATCTTTTCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                   <  1:2830579/65‑1 (MQ=255)
                                          ttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACt                  >  1:834494/1‑65 (MQ=255)
                                              tttCTTTCCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:693839/63‑1 (MQ=255)
                                              tttCTTTCCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:1575281/63‑1 (MQ=255)
                                              tatTTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:349940/65‑1 (MQ=255)
                                              tatTTTTTCACGGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:735459/65‑1 (MQ=255)
                                              tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:3368417/65‑1 (MQ=255)
                                              tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:344257/65‑1 (MQ=255)
                                              tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:2716304/65‑1 (MQ=255)
                                              tatCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:3694145/65‑1 (MQ=255)
                                              tatCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              >  1:1647726/1‑65 (MQ=255)
                                              tatCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:1703076/65‑1 (MQ=255)
                                              tatCTTTTCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:1951409/65‑1 (MQ=255)
                                              tatCTTTTCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:1288830/65‑1 (MQ=255)
                                              tatCTTTCCTAAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:204610/65‑1 (MQ=255)
                                              catCTTTCCTCGGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:459169/64‑1 (MQ=255)
                                              attCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt              <  1:3722186/63‑1 (MQ=255)
                                                tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTcc            >  1:226543/1‑65 (MQ=255)
                                                  ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg          <  1:608309/65‑1 (MQ=255)
                                                  ttttCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg          <  1:3045630/65‑1 (MQ=255)
                                                  ttttCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg          <  1:3847329/65‑1 (MQ=255)
                                                   tttCTCAGGGCCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc         <  1:2056201/65‑1 (MQ=255)
                                                   tttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc         <  1:3257124/65‑1 (MQ=255)
                                                   tttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc         <  1:3374785/65‑1 (MQ=255)
                                                   tttCCCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc         <  1:138053/65‑1 (MQ=255)
                                                   tttCACAGGCCCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc         <  1:349204/65‑1 (MQ=255)
                                                   ttcctcGGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc         <  1:3182853/62‑1 (MQ=255)
                                                    ttCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCg        <  1:219130/65‑1 (MQ=255)
                                                     tctcAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGt       <  1:3366778/65‑1 (MQ=255)
                                                        cAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGtata    <  1:2167170/65‑1 (MQ=255)
                                                        cAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGtata    <  1:2160181/65‑1 (MQ=255)
                                                          ggAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATAcc  >  1:3210277/1‑65 (MQ=255)
                                                             |                                                             
ATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATACC  >  NC_000913/310873‑310995

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