Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,484,1050AG37.5% 66.6 / 26.3 48G55G (GGA→GGG) ydcUpredicted spermidine/putrescine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (5/25);  new base G (18/0);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 4.60e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA  >  REL606/1484030‑1484179
                                                                           |                                                                          
tttACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGa                                                                            <  2:1850014/76‑1 (MQ=255)
 ttACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAt                                                                           <  1:608032/76‑1 (MQ=255)
    ctgctgCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGtttt                                                                        >  2:1799802/1‑76 (MQ=255)
     tgctgcTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGttttt                                                                       >  1:406249/1‑75 (MQ=255)
     tgctgcTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTta                                                                       <  1:178428/76‑1 (MQ=255)
        tgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTAc                                                                    >  2:726540/1‑76 (MQ=255)
        tgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATAc                                                                    <  1:1119925/76‑1 (MQ=255)
         gctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACt                                                                   >  1:986671/1‑76 (MQ=255)
          ctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACtt                                                                  >  2:1479783/1‑76 (MQ=255)
          ctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACtt                                                                  <  1:1341810/76‑1 (MQ=255)
            tGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACCCTGTTATGGCAGGGGGTTTTTTCtttt                                                                >  2:242009/1‑76 (MQ=255)
              gCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGa                                                              <  1:414882/76‑1 (MQ=255)
                ccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACg                                                            <  2:183757/76‑1 (MQ=255)
                 cACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACCCTGTTATGGCAGGGGGTTTTTCCTTTTCACGa                                                           >  1:1789727/1‑76 (MQ=255)
                    tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTTACGAttt                                                        >  1:1657910/1‑76 (MQ=255)
                     aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTACGAtttt                                                       >  1:366123/1‑75 (MQ=255)
                     aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGAtttt                                                       >  2:1298106/1‑75 (MQ=255)
                        gtgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTACGGTTTTAcc                                                    >  2:1771098/1‑76 (MQ=255)
                        gtgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCAcc                                                    <  1:1150777/76‑1 (MQ=255)
                         tgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTGTCGATTTTACCa                                                   >  2:1694007/1‑76 (MQ=255)
                         tgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTAATTTTTTTGTTTTCACCa                                                   >  1:3395/1‑76 (MQ=255)
                         tgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCa                                                   <  2:1716805/76‑1 (MQ=255)
                          ggTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTTACCAt                                                  >  2:483759/1‑76 (MQ=255)
                              tGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCg                                              <  2:1587900/76‑1 (MQ=255)
                                gCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGt                                            <  1:999636/76‑1 (MQ=255)
                                 cATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTa                                           <  2:950906/76‑1 (MQ=255)
                                  aTTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTACGATTTCACCATGTTGGGaa                                          >  2:957090/1‑76 (MQ=255)
                                  aTTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTaa                                          <  2:468365/76‑1 (MQ=255)
                                   ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAAc                                         <  2:1832694/76‑1 (MQ=255)
                                    tGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACg                                        <  1:1041998/76‑1 (MQ=255)
                                        tATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCATGTCGGGAACGCCgg                                    >  2:1440837/1‑76 (MQ=255)
                                        tATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCgg                                    <  2:1286970/76‑1 (MQ=255)
                                           ttCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAAc                                 <  1:646163/76‑1 (MQ=255)
                                                cTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACg                            <  2:1446062/76‑1 (MQ=255)
                                                 tCGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGc                           >  2:1086137/1‑76 (MQ=255)
                                                  cGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCt                          >  1:274991/1‑76 (MQ=255)
                                                  cGCTGCTGACACTGTTATGGCAGGGATTTTTTACTTTTGACGATTTCACCATGTCGGGAACGCCGGAACTGACGCt                          >  1:1018386/1‑76 (MQ=255)
                                                      gctgACACTGTTATGGCAGGGGTTTTATACTTTTTACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCg                      >  2:115948/1‑76 (MQ=255)
                                                      gctgACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCg                      >  2:844248/1‑76 (MQ=255)
                                                      gctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCg                      <  1:1166573/76‑1 (MQ=255)
                                                          acacTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAata                  <  2:1319055/76‑1 (MQ=255)
                                                           cacTGTTATGGCAGGGGTTTTATACTTTTGACGATTTc                                                       >  2:1602114/1‑38 (MQ=255)
                                                           cacTGTTATGGCAGGGATTTTATACTTTTGACGATTTc                                                       <  1:1602114/38‑1 (MQ=255)
                                                            acTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGGGAATATc                >  2:466897/1‑76 (MQ=255)
                                                             cTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATcc               <  1:1179140/76‑1 (MQ=255)
                                                              tGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCg              <  1:146841/76‑1 (MQ=255)
                                                               gTTATGGCAGGGGTTTTTTACTTTTTACGATTTCACCATGTCGGTAACGCCGGGACTGACGCTGGCGAATATccct             >  2:1724837/1‑74 (MQ=255)
                                                                ttATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGGAACGCCGGAAATGACGCCGGCGAATATCCGTg            >  2:1083338/1‑76 (MQ=255)
                                                                ttATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTg            <  2:1639629/76‑1 (MQ=255)
                                                                 tATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGGCGCTGGCGAATATCCGTgc           <  1:726540/76‑1 (MQ=255)
                                                                 tATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgc           >  2:1105816/1‑76 (MQ=255)
                                                                      cAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGt      >  1:679941/1‑76 (MQ=255)
                                                                        gggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGttt    <  1:756581/76‑1 (MQ=255)
                                                                         ggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTa   <  2:441742/76‑1 (MQ=255)
                                                                          gATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTaa  >  1:511819/1‑76 (MQ=255)
                                                                           |                                                                          
TTTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA  >  REL606/1484030‑1484179

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑