Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,335,1610TG34.2% 41.3 / 11.3 38intergenic (‑30/+29)yjdL/cadApredicted transporter/lysine decarboxylase 1
Reads supporting (aligned to +/- strand):  ref base T (5/20);  new base G (0/13);  total (5/33)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTC  >  REL606/4335111‑4335211
                                                   |                                                  
cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGt                                                    <  1:384140/51‑1 (MQ=255)
  gCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTgg                                                  <  1:2131633/51‑1 (MQ=255)
    tgtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCa                                                <  1:1678768/51‑1 (MQ=255)
    tgtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCa                                                <  1:2395914/51‑1 (MQ=255)
      tgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAg                                              <  1:5879666/51‑1 (MQ=255)
        aGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCaa                                               >  1:6454308/1‑48 (MQ=255)
        aGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                            <  1:2093806/51‑1 (MQ=255)
                   aTGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGt                                 <  1:5365406/51‑1 (MQ=255)
                    tgtaTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTa                                <  1:3304997/51‑1 (MQ=255)
                         tctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGc                           <  1:5115531/51‑1 (MQ=255)
                              ttATGAGAAAAAAAGGGAAGGGGAAAGCCACTTCCCTTGTACGAGCTAatt                      <  1:5524840/51‑1 (MQ=255)
                                  gAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttt                  <  1:6906414/51‑1 (MQ=255)
                                  gAGCAAAAAAGGGAAGGGGAAAGCCACTTCCCTTGTACGAGCTna                        <  1:6601976/45‑1 (MQ=255)
                                      aaaaagGGGAGGGGGCAAGCCACTTCCCTTGTACGAGCt                          <  1:1548118/39‑1 (MQ=255)
                                      aaaaagGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttttt                <  1:6137118/49‑1 (MQ=255)
                                      aaaaagGGAAAGGGGAAAGCCACTTCCCTTGTACGAGCTaa                        <  1:14813/41‑1 (MQ=255)
                                      aaaaaaGGGAGGGGGAAAGCCACTTCCCTTGTACGAGCTAATTAtttttt               <  1:266969/50‑1 (MQ=255)
                                      aaaaaaGGGAAGGGGGAAAGCCACTTCCCTTGTACGAGCTAATTAttt                  <  1:664767/48‑1 (MQ=255)
                                      aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAattatt                   <  1:3039265/46‑1 (MQ=255)
                                      aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt               <  1:4240626/50‑1 (MQ=255)
                                      aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt               <  1:4600036/50‑1 (MQ=255)
                                      aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttttt                <  1:4683655/49‑1 (MQ=255)
                                      aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTg              <  1:3891388/51‑1 (MQ=255)
                                      aaaaaaGGGAAGTGGAAAGCCACTTCCCTTGTACGAGCTAATTAttttt                <  1:4224562/49‑1 (MQ=255)
                                      aaaaaaGGGAAGTGGAAAGCCACTTCCCTTGTACGAGCTAATTAtttt                 <  1:5568120/48‑1 (MQ=255)
                                      aaaaaaGGGAAGTGGAAAGCCACTTCCCTTGTACGAGCTAATTAttt                  <  1:2918992/47‑1 (MQ=255)
                                      aaaaaaGGGAAGGGGCAAGCCACTTCCCTTGTACGAGCTAattat                    <  1:4415142/45‑1 (MQ=255)
                                      aaaaaaGGGAAGGGGAAAGCCACTTCCCTTGTACGAGc                           <  1:2119312/38‑1 (MQ=38)
                                      aaaaaaGGGAAGGGGAAAGCCACTTCCCTTGTACGAGCTAatta                     <  1:1771533/44‑1 (MQ=38)
                                      aaaaaaGGGAAGGGGAAAGCCACTTCCCTTGTACGAGCTAATTAtttttt               <  1:736136/50‑1 (MQ=255)
                                      aaaaaaGGGAAGGGGAAAGCCACTTCCCTTGTACGAGCTAATTAttttt                <  1:2220817/49‑1 (MQ=255)
                                       aaaaaGGGAAGGGGGAAAGCCACTTCCCTTGTACGAGCTaa                        <  1:279505/41‑1 (MQ=255)
                                          aaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt               >  1:4110775/1‑46 (MQ=255)
                                            gggAAGGGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttct        <  1:2431038/51‑1 (MQ=255)
                                               aaGTGGCAAGCCACTTGCCTTGTACGAGCTAATTATTTTTTGCTTtcttct     <  1:6083750/51‑1 (MQ=255)
                                               aaGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct     >  1:3684460/1‑51 (MQ=255)
                                               aaGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttct     >  1:7152164/1‑51 (MQ=255)
                                                aGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTtcttctc    <  1:4202587/51‑2 (MQ=255)
                                                aGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTcttctt    >  1:3117760/1‑51 (MQ=255)
                                                   tGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc  >  1:2146139/1‑51 (MQ=255)
                                                   |                                                  
CGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTC  >  REL606/4335111‑4335211

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: