Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913310,9340AC57.1% ‑4.3 / 31.0 28F135V (TTC→GTC) ecpRDNA‑binding transcriptional dual regulator MatA
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base C (0/16);  total (7/21)
Fisher's exact test for biased strand distribution p-value = 6.69e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATA  >  NC_000913/310868‑310993
                                                                  |                                                           
ccGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAAc                                                            >  1:2097843/1‑68 (MQ=255)
      tGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCa                                                      >  1:2976239/1‑68 (MQ=255)
             aaTTTCCATTTCCCGGTCCGTGTTTTTATCTTTCTAAATATTTTTTCTCAGGACCCTGCCATTAATca                                               <  1:3380048/68‑1 (MQ=255)
              atttccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATcac                                              <  1:883277/68‑1 (MQ=255)
              atttccatttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATcac                                              <  1:2932034/68‑1 (MQ=255)
              atttccattcccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATcac                                              <  1:2066998/68‑1 (MQ=255)
              attgccattcccGGGACGGGGATTTTATCTTTTTTAATATCTTTTCTCAGGACCCTGCCATTAATcac                                              <  1:706517/68‑1 (MQ=255)
                   catttccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGAAACCTGCCATTAATCACATAGa                                         <  1:3621152/68‑1 (MQ=255)
                   cattcccCGGTCCTTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGa                                         <  1:2110041/68‑1 (MQ=255)
                   cattcccCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGa                                         <  1:1534559/68‑1 (MQ=255)
                        ccgggTCCGTGATTTTATCTTTCTAAACATCTTTTCCCAGGAACCTGCCATTAATCACATAGACCAGt                                    <  1:909044/65‑1 (MQ=255)
                        cccGGTCCTTGTTTTTATTTTTCTTAATATCTTTTCGCAGGACCCTGCCATTAATCACATAGACCAGt                                    <  1:860420/68‑1 (MQ=255)
                        cccGGTCCGTGTTTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGt                                    <  1:1374758/68‑1 (MQ=255)
                        cccGGTCCGTGTTTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGt                                    <  1:1976593/68‑1 (MQ=255)
                        cccGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGt                                    <  1:1631085/68‑1 (MQ=255)
                        cccGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGt                                    >  1:456982/1‑68 (MQ=255)
                        cccGGTCCGTGATTTTATCTTTCTAAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGt                                    <  1:4156066/68‑1 (MQ=255)
                        cccGGCCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGAAACCTGCCATTAATCACATAGACCAGt                                    <  1:1837018/68‑1 (MQ=255)
                         ccGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGtt                                   >  1:2739883/1‑68 (MQ=255)
                             tCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCttt                               >  1:1055669/1‑68 (MQ=255)
                              ccGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCtttt                              <  1:3860557/68‑1 (MQ=255)
                                    ttttttGTTTCTTAATCTCTTTTCTCAGCACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAata                        <  1:1859033/68‑1 (MQ=255)
                                    ttttATCTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAata                        <  1:2569769/68‑1 (MQ=255)
                                      ttATTTTTCTTAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATc                      <  1:2780177/68‑1 (MQ=255)
                                      ttATCTTTCTTAATTTCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATc                      <  1:3072122/68‑1 (MQ=255)
                                         tctttcttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACg                   >  1:1897590/1‑68 (MQ=255)
                                         tctttcttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACg                   <  1:2168961/68‑1 (MQ=255)
                                          ctttcttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGNTCTTNNCTAATATCNCNa                  >  1:1100332/1‑68 (MQ=255)
                                           tttcttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                 <  1:409071/68‑1 (MQ=255)
                                           tttcttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                 <  1:3979306/68‑1 (MQ=255)
                                           tttcttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                 <  1:1439249/68‑1 (MQ=255)
                                           tttcttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                 <  1:3074013/68‑1 (MQ=255)
                                           tttcttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                 <  1:1122902/68‑1 (MQ=255)
                                           tttcttAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGAc                 <  1:430630/68‑1 (MQ=255)
                                              cttAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAGTATCACGACTCa              <  1:2589862/68‑1 (MQ=255)
                                                tattatCTTTCCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt            <  1:1504530/65‑1 (MQ=255)
                                                tAATTTCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt            <  1:1062359/68‑1 (MQ=255)
                                                tAATTTCTTTTCGGGGAAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt            <  1:1193859/68‑1 (MQ=255)
                                                tAATATTTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt            <  1:2263215/68‑1 (MQ=255)
                                                tAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt            <  1:3054462/68‑1 (MQ=255)
                                                tAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt            <  1:1882422/68‑1 (MQ=255)
                                                tAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt            <  1:857466/68‑1 (MQ=255)
                                                tAATATCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt            <  1:830699/68‑1 (MQ=255)
                                                tAATATCTTTTCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt            <  1:481802/68‑1 (MQ=255)
                                                aaaTTTCTTTTCCCAGACACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGt            <  1:3696668/67‑1 (MQ=255)
                                                    atCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg        <  1:3043919/68‑1 (MQ=255)
                                                    atCTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg        <  1:3642223/68‑1 (MQ=255)
                                                    atCTTTTCACAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg        <  1:2131369/68‑1 (MQ=255)
                                                    atCTTTCCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAg        <  1:1015083/68‑1 (MQ=255)
                                                     tttttcCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCGCGACTCAGTCCAGc       <  1:1772475/66‑1 (MQ=255)
                                                     tCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc       >  1:1613088/1‑68 (MQ=255)
                                                     tCTTTTCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc       <  1:2428323/68‑1 (MQ=255)
                                                     tCTTTTCTAAGGAACCTACCATTAATCACATAGTCCAGTTCTTTTCTAATATCACGACTCAGTCCAGc       <  1:3682455/68‑1 (MQ=255)
                                                     tCTTTTCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGc       <  1:3458022/68‑1 (MQ=255)
                                                      cTTTTCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCg      <  1:3886014/68‑1 (MQ=255)
                                                      cTTTTCCCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCg      <  1:3251097/68‑1 (MQ=255)
                                                       ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGt     <  1:3854705/68‑1 (MQ=255)
                                                       ttttCTCAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGt     <  1:316306/68‑1 (MQ=255)
                                                       ttttCTCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGt     <  1:2524780/68‑1 (MQ=255)
                                                       ttttCCCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTGATATCACGACTCAGTCCAGCGt     <  1:1363091/68‑1 (MQ=255)
                                                       ttttCCCAGAACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGt     <  1:2004979/68‑1 (MQ=255)
                                                          tctcAGGACCCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGtata  <  1:3185396/68‑1 (MQ=255)
                                                                  |                                                           
CCGTCATGCGGATAATTTCCATTTCCCGGTCCGTGATTTTATCTTTCTTAATATCTTTTCTCAGGAACCTGCCATTAATCACATAGACCAGTTCTTTTCTAATATCACGACTCAGTCCAGCGTATA  >  NC_000913/310868‑310993

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