Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,484,1050AG32.0% 32.3 / 11.6 25G55G (GGA→GGGydcUpredicted spermidine/putrescine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (4/13);  new base G (8/0);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 4.58e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA  >  REL606/1484031‑1484179
                                                                          |                                                                          
ttACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAt                                                                           <  1:1114683/76‑1 (MQ=255)
ttACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAt                                                                           <  1:970918/76‑1 (MQ=255)
     gctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGtttttt                                                                      >  1:602562/1‑74 (MQ=255)
     gctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTtat                                                                      <  2:876368/76‑1 (MQ=255)
         ctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACtt                                                                  <  1:644683/76‑1 (MQ=255)
           tGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTTCtttt                                                                >  2:775766/1‑76 (MQ=255)
            ggCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTg                                                               <  1:219862/76‑1 (MQ=255)
              cccACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTAAttttttc                                                             >  1:109719/1‑73 (MQ=255)
                   tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTTTTCTTTTTTCGAttt                                                        >  2:59511/1‑76 (MQ=255)
                    aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTACGGtttt                                                       >  1:736095/1‑75 (MQ=255)
                       gtgGTTTGGCATTGTCTATTTCGGCTCGCTGCTGATACTGTTATGGCAGGGGTTTTTTTCTTTTTTCGATTTTAcc                                                    >  1:1036646/1‑76 (MQ=255)
                         ggTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCAt                                                  <  2:58979/76‑1 (MQ=255)
                          gTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTCCCATg                                                 >  2:247738/1‑76 (MQ=255)
                                 aTTGTCTATGTCGGCTCGCTGCTGACACTGTTATGGCTGGGATTTTATACTTTTGACGATTTCACCATGTCGGTaa                                          <  1:225796/76‑1 (MQ=255)
                                      cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTTCGATTTTACCCTGTCGGGAAAGccc                                     >  1:829138/1‑75 (MQ=255)
                                           tCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACt                                <  2:821219/76‑1 (MQ=255)
                                             ggCTCGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACCCCGGAACTGa                              >  2:991135/1‑76 (MQ=255)
                                               cTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACg                            <  1:824898/76‑1 (MQ=255)
                                                    tgctgACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCCGGc                       >  2:587820/1‑76 (MQ=255)
                                                     gctgACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTCACGCCGGAACTGACGCTGGCg                      >  2:718424/1‑76 (MQ=255)
                                                        gACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAat                   <  2:691991/76‑1 (MQ=255)
                                                         acacTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAata                  <  1:89210/76‑1 (MQ=255)
                                                          cacTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAatat                 <  2:479131/76‑1 (MQ=255)
                                                            cTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCCTGTCAGTAACGCCGGAACTGACGCTGGCGAATATcc               >  2:78776/1‑76 (MQ=255)
                                                               ttATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGGAACGCCGGAACTGACGCCGGCGAATATCCGTg            >  2:352532/1‑76 (MQ=255)
                                                                     cAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGt      >  1:576630/1‑76 (MQ=255)
                                                                     cAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCACTGCGCTGt      >  1:295005/1‑76 (MQ=255)
                                                                         gATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTaa  <  2:89750/76‑1 (MQ=255)
                                                                          |                                                                          
TTACTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA  >  REL606/1484031‑1484179

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