Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,484,1050AG36.1% 59.2 / 16.2 36G55G (GGA→GGGydcUpredicted spermidine/putrescine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (2/21);  new base G (13/0);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 4.54e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA  >  REL606/1484033‑1484179
                                                                        |                                                                          
aCTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAttt                                                                         <  2:963305/76‑1 (MQ=255)
  tgctgcTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTta                                                                       <  2:1105629/76‑1 (MQ=255)
    ctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGGCACTGTTATGGCAGGGGttttttt                                                                     >  2:249115/1‑73 (MQ=255)
       ctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACtt                                                                  <  2:1338270/76‑1 (MQ=255)
        ttGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACttt                                                                 <  2:1109205/76‑1 (MQ=255)
         tGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCtttt                                                                >  2:682457/1‑76 (MQ=255)
         tGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTACtttt                                                                >  1:1095663/1‑76 (MQ=255)
            cccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACttttttc                                                             >  2:456119/1‑73 (MQ=255)
            cccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGAc                                                             <  2:939209/76‑1 (MQ=255)
                cTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGAtt                                                         >  1:1817608/1‑76 (MQ=255)
                cTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGAtt                                                         <  1:182011/76‑1 (MQ=255)
                 tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTGCGAttt                                                        >  2:1923161/1‑76 (MQ=255)
                 tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTACGAttt                                                        >  2:213859/1‑76 (MQ=255)
                 tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTATTTTTTCGAttt                                                        >  1:316540/1‑76 (MQ=255)
                   aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTATTTTTGAGGAatttc                                                      >  1:557353/1‑74 (MQ=255)
                    tgtgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTAc                                                     >  1:1545294/1‑76 (MQ=255)
                       ggTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCAccct                                                  >  2:1143631/1‑74 (MQ=255)
                         tttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGt                                                <  1:1817511/76‑1 (MQ=255)
                         tttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGt                                                <  1:1512254/76‑1 (MQ=255)
                          ttGGCATTTTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTATTTTTATGATTTTACCATGTc                                               >  1:52959/1‑76 (MQ=255)
                           tGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGATTTTACCATTTCg                                              >  1:576566/1‑76 (MQ=255)
                           tGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCACCATGTCg                                              >  1:553998/1‑76 (MQ=255)
                           tGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCg                                              <  2:755455/76‑1 (MQ=255)
                           tGGCACTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTACGATTTTACCATGTCg                                              >  1:1306315/1‑76 (MQ=255)
                             gCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGt                                            <  1:244620/76‑1 (MQ=255)
                              cATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTa                                           <  1:1121388/76‑1 (MQ=255)
                                            gCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGAc                             <  1:294473/76‑1 (MQ=255)
                                               cGCTGCTGACACTGTTGTGGCAGGGGTTTTTTACTTTTGACGATTTCACCATGTCGGGAACGCCGGAACTGACGCt                          >  1:1549401/1‑76 (MQ=255)
                                                gctgctGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTg                         <  2:1727011/76‑1 (MQ=255)
                                                gctgctGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTg                         <  2:538162/76‑1 (MQ=255)
                                                 ctgctgACACTGTTATGGCAGGGGTTTTTTACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTgg                        >  1:1148605/1‑76 (MQ=255)
                                                   gctgACACTGTTATGGCAGGGGTTTTTTACTTTTGACGATTTCACCATGTCGGTAACGCCGGAAATGACGCTGGCg                      >  1:874968/1‑76 (MQ=255)
                                                    ctgACACTGTTATGGCAGGGGTTTTATACTTTTTACGATTTCACCATGTCGGTAACGCCGGAACTGACGCCGGGGa                     >  1:1820269/1‑76 (MQ=255)
                                                     tgACACTGTTATGGCAGGGGTTTTTTACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGaa                    >  2:297685/1‑76 (MQ=255)
                                                      gACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAat                   <  1:554925/76‑1 (MQ=255)
                                                        cacTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAatat                 <  2:1860757/76‑1 (MQ=255)
                                                        cacTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAatat                 <  2:948764/76‑1 (MQ=255)
                                                            gTTATGGCAGGGGTTTTTTACTTTTGACGATTTCACCATGTCGGGAAACCCCGAACTGACGCTGGCGGATATccct             >  2:444839/1‑74 (MQ=255)
                                                             ttATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTg            <  1:1890150/76‑1 (MQ=255)
                                                              tATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgc           <  1:1808615/76‑1 (MQ=255)
                                                               aTGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgcg          <  1:1267334/76‑1 (MQ=255)
                                                               aTGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgcg          >  1:1935864/1‑76 (MQ=255)
                                                                   cAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGt      >  2:1939911/1‑76 (MQ=255)
                                                                       gATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTaa  <  2:300841/76‑1 (MQ=255)
                                                                        |                                                                          
ACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA  >  REL606/1484033‑1484179

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: