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

CTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA  >  REL606/1484034‑1484179
                                                                       |                                                                          
ctgctgCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAttt                                                                         <  2:196589/75‑1 (MQ=255)
 tgctgcTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGAgttt                                                                        >  2:900500/1‑75 (MQ=255)
  gctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTta                                                                       <  1:771580/75‑1 (MQ=255)
  gctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTta                                                                       <  2:465574/75‑1 (MQ=255)
   ctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTtat                                                                      <  2:1396581/75‑1 (MQ=255)
   ctgctTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTtat                                                                      <  2:1554682/75‑1 (MQ=255)
       ttGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTAttt                                                                  >  1:1345633/1‑75 (MQ=255)
       ttGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTAAtt                                                                  >  1:167259/1‑75 (MQ=255)
       ttGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGGTTTATTCtt                                                                  >  1:1936558/1‑75 (MQ=255)
         ggCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACtttt                                                                <  1:40241/75‑1 (MQ=255)
                tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTATTTTTTTCGAtt                                                         >  2:365122/1‑75 (MQ=255)
                tAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGAtt                                                         <  2:1584885/75‑1 (MQ=255)
                 aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTGACGAgtt                                                        >  2:2096905/1‑75 (MQ=255)
                 aaTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGAttt                                                        <  1:1393325/75‑1 (MQ=255)
                  aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTc                                                       <  2:2018897/75‑1 (MQ=255)
                  aTGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTc                                                       <  2:164143/75‑1 (MQ=255)
                      ggTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCa                                                   <  1:792176/75‑1 (MQ=255)
                      ggTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCa                                                   <  1:1393364/75‑1 (MQ=255)
                       gTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTAccct                                                  >  1:582289/1‑73 (MQ=255)
                        tttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTATTTTTTCGATTTTACCCTg                                                 >  1:1250466/1‑75 (MQ=255)
                        tttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTATATTTTTGACGATTTCACCATg                                                 >  2:2123474/1‑75 (MQ=255)
                         ttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGt                                                <  1:1093146/75‑1 (MQ=255)
                         ttGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGt                                                <  2:1225935/75‑1 (MQ=255)
                           ggCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTCACCATGTCg                                              >  2:195289/1‑75 (MQ=255)
                            gCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCgg                                             <  1:1618921/75‑1 (MQ=255)
                               ttttCTATTGCGGCCCGCCGCAGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTaa                                          <  2:688798/72‑1 (MQ=255)
                               ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTGAGGATTTCACCATGTCGGTaa                                          >  2:1806252/1‑75 (MQ=255)
                               ttGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTaa                                          <  2:688792/75‑1 (MQ=255)
                                tGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGGTTTTACCATGTCGGGaaa                                         >  1:1629300/1‑74 (MQ=255)
                                  tCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTTCGATTTTACCATGTCGGGAACGc                                       >  2:66296/1‑75 (MQ=255)
                                   cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTACGATTTTACCCTGTCGGGAAAGcc                                      >  1:1102460/1‑75 (MQ=255)
                                   cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTTAAGATTTTACCATGTCGGTAACGcc                                      >  1:1968634/1‑75 (MQ=255)
                                   cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGcc                                      <  1:1401247/75‑1 (MQ=255)
                                   cTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGcc                                      <  1:898760/75‑1 (MQ=255)
                                    tATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTTCTTTTTACGATTTTACCCTGTCGGGAAAGCCg                                     >  1:843695/1‑75 (MQ=255)
                                     aTTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTGACGATTTCACCATGTCGGGAACGCCgg                                    >  1:250294/1‑75 (MQ=255)
                                     aTTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCgg                                    <  1:681160/75‑1 (MQ=255)
                                       ttCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGaa                                  <  1:1134446/75‑1 (MQ=255)
                                       ttCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGaa                                  <  2:1192331/75‑1 (MQ=255)
                                       ttCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGaa                                  <  2:1789199/75‑1 (MQ=255)
                                        tCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAAc                                 <  2:1110514/75‑1 (MQ=255)
                                         cGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACt                                <  2:2152668/75‑1 (MQ=255)
                                         cGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACt                                <  1:2019049/75‑1 (MQ=255)
                                          ggCTCGCTGCTGACACTGTTATGGCAGGGGTTTTTTAATTTTTACGATTTTACCATGTCGGGAACCCCGGGACTg                               >  1:1585680/1‑75 (MQ=255)
                                           gCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGa                              <  1:1766206/75‑1 (MQ=255)
                                            cTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGAc                             <  1:916577/75‑1 (MQ=255)
                                             tCGCTGCTGACACTGTTATGGCAGGGGTTTTTTACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACg                            >  2:2093641/1‑75 (MQ=255)
                                             tCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACg                            >  1:365276/1‑75 (MQ=255)
                                             tCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACg                            <  1:2082048/75‑1 (MQ=255)
                                              cGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGc                           <  2:1345633/75‑1 (MQ=255)
                                              cGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGc                           <  2:857051/75‑1 (MQ=255)
                                                 tgctgACACTGTTATGGCAGGGGTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTgg                        >  1:479191/1‑75 (MQ=255)
                                                 tgctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTgg                        <  1:637389/75‑1 (MQ=255)
                                                  gctgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGc                       >  2:373806/1‑75 (MQ=255)
                                                    tgACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGa                     >  1:717110/1‑75 (MQ=255)
                                                         cTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATc                >  2:884368/1‑75 (MQ=255)
                                                         cTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATc                <  2:1002657/75‑1 (MQ=255)
                                                           gTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCg              >  2:1564202/1‑75 (MQ=255)
                                                           gTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCg              <  2:2124479/75‑1 (MQ=255)
                                                            ttATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGt             <  1:715379/75‑1 (MQ=255)
                                                                  cAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTg       >  1:997405/1‑75 (MQ=255)
                                                                  cAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTg       >  1:1112395/1‑75 (MQ=255)
                                                                   aGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGt      <  2:1968634/75‑1 (MQ=255)
                                                                   aGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGt      >  2:864089/1‑75 (MQ=255)
                                                                    gggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGtt     <  2:1146296/75‑1 (MQ=255)
                                                                     ggATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGttt    <  2:2003103/75‑1 (MQ=255)
                                                                      gATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTa   >  1:1947063/1‑75 (MQ=255)
                                                                      gATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTa   <  1:8987/75‑1 (MQ=255)
                                                                      gATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTa   <  1:416069/75‑1 (MQ=255)
                                                                       aTTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTaa  <  1:188378/75‑1 (MQ=255)
                                                                       |                                                                          
CTGCTGCTTGGCCCACTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGATTTCACCATGTCGGTAACGCCGGAACTGACGCTGGCGAATATCCGTGCGCTGTTTAA  >  REL606/1484034‑1484179

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