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

CTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAAT  >  REL606/1484034‑1484180
                                                                       |                                                                           
ctgctgCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGtttt                                                                         >  1:1437043/1‑76 (MQ=255)
  gctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTtat                                                                       <  1:1678220/76‑1 (MQ=255)
      ctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACtt                                                                   <  1:803794/76‑1 (MQ=255)
        tGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACtttt                                                                 <  1:1307837/76‑1 (MQ=255)
         ggCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTg                                                                <  1:1315336/76‑1 (MQ=255)
            ccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACg                                                             <  1:829822/76‑1 (MQ=255)
             cACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTATTTTTGTCGa                                                            >  1:1600666/1‑76 (MQ=255)
                tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGAttt                                                         >  1:497664/1‑76 (MQ=255)
                tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTAATTTTTACGAttt                                                         >  2:550695/1‑76 (MQ=255)
                 aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGACAGGGGTTTTTTAATTTTTGAGATTTc                                                        >  2:747870/1‑76 (MQ=255)
                  aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCa                                                       >  1:71256/1‑76 (MQ=255)
                    gtgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCAcc                                                     <  1:95088/76‑1 (MQ=255)
                      ggTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTTCGATTTTACCAt                                                   >  1:144065/1‑76 (MQ=255)
                        tttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTAGGATTTTACCCTGt                                                 >  2:1176840/1‑76 (MQ=255)
                        tttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGATTTTACCATGt                                                 >  1:1248455/1‑76 (MQ=255)
                        tttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCATGt                                                 >  1:1518819/1‑76 (MQ=255)
                             cATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTAATTTTTTCGATTTTCCCATGTCggga                                            >  1:346049/1‑74 (MQ=255)
                             cATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTa                                            <  2:1311297/76‑1 (MQ=255)
                             cATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTa                                            <  2:465764/76‑1 (MQ=255)
                              aTTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTaa                                           <  1:374989/76‑1 (MQ=255)
                               ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGATTTTACCATTTCGGTaaa                                          >  1:1199478/1‑75 (MQ=255)
                                  tCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGcc                                       <  1:470500/76‑1 (MQ=255)
                                   cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCg                                      <  2:1061691/76‑1 (MQ=255)
                                       ttCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAAc                                  <  1:373393/76‑1 (MQ=255)
                                          ggCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGa                               <  1:896161/76‑1 (MQ=255)
                                          ggCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGa                               <  1:227722/76‑1 (MQ=255)
                                           gCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTTTTTTGACGGTTTCACCATGTCGGGAACGCCGGAACTGAc                              >  1:243604/1‑76 (MQ=255)
                                              cGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCAAGTCCGGAACGCCGGGACTGACGcg                           >  2:65329/1‑75 (MQ=255)
                                              cGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTTACGATTTCACCATGTCGGTAACGCCGGAACTGACGCt                           >  2:1489543/1‑76 (MQ=255)
                                              cGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCt                           >  2:117742/1‑76 (MQ=255)
                                                 tgctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGc                        <  1:1150583/76‑1 (MQ=255)
                                                  gctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAAAGCCGGGACTGACGCTGGCg                       >  2:1666158/1‑76 (MQ=255)
                                                    tgACACTGTTATGGCAGGGATTGTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGaa                     <  1:699100/76‑1 (MQ=255)
                                                           gTTATGGCAGGGGTTTTTTACTTTTGACGATTTCACCATGTCGGTAACGCCGGAAATGACGCTCGCGAATATCCgg              >  1:1160236/1‑75 (MQ=255)
                                                               tGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTgcgc          >  2:380898/1‑76 (MQ=255)
                                                                ggCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCt         <  1:1303466/76‑1 (MQ=255)
                                                                    gggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGttt     <  2:193002/76‑1 (MQ=255)
                                                                       aTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAAt  <  1:831922/76‑1 (MQ=255)
                                                                       aTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAAt  <  1:613644/76‑1 (MQ=255)
                                                                       |                                                                           
CTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAAT  >  REL606/1484034‑1484180

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