Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176351,637,7570GA41.7% 33.2 / 29.4 24E599E (GAG→GAA) ECW_RS08125RHS element core protein
Reads supporting (aligned to +/- strand):  ref base G (8/6);  new base A (4/6);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCAGTATGAGTACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGG  >  NC_017635/1637710‑1637801
                                               |                                            
cGGTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCAGTATGAGTa                                            >  2:127549/1‑50 (MQ=255)
 ggTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCGGTATGAATa                                            >  2:1953397/1‑49 (MQ=11)
 ggTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCGGTATGAATa                                            >  1:1334104/1‑49 (MQ=11)
   ttAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCAGTATGAGTACaa                                         >  2:844644/1‑50 (MQ=255)
      aCCTCGGTGAAAGACGCACAGGGCCGTGAAACGCAGTATGAGTACAAcg                                       <  1:817383/49‑1 (MQ=255)
          cGGTGAAAGACGCACAGGGCCGTGAAACGCGGTATGAATACAAcgccg                                    >  2:388735/1‑48 (MQ=11)
           ggTGAAAGACGCACAGGGCCGTGAAACGCAGTATGAGTACAACGCCGCAg                                 >  1:1749365/1‑50 (MQ=255)
            gTGAAAGACGCACAGGGCCGTGAAACGCAGTATGAGTACAACGCCGCAgg                                >  1:58450/1‑50 (MQ=255)
                  gACGCACAGGGCCGTGAAACGCGGTATGAATACAACGCCGCAGGCGACCt                          <  2:1370353/50‑1 (MQ=14)
                   aCGCACAGGGCCGTGAAACGCAGTATGAGTACAACGCCGCAGGCGACCTg                         <  2:856716/50‑1 (MQ=255)
                    cGCACAGGGCCGTGAAACGCAGTATGAGTACAACGCCGCAGGCGACCTGa                        >  1:1526129/1‑50 (MQ=255)
                      cacaGGGCCGTGAAACGCGGTATGAATACAACGCCGCAGGCGACCTGACt                      <  2:1569035/50‑1 (MQ=14)
                       acaGGGCCGTGAAACGCAGTATGAGTACAACGCCGCAGGCGACCTGACTg                     <  1:1783426/50‑1 (MQ=255)
                       acaGGGCCGTGAAACGCAGTATGAGTACAACGCCGCAGGCGACCTGACTg                     <  2:1196764/50‑1 (MQ=255)
                         aGGGCCGTGAAACGCAGTATGAGAACAACGCCGCAGGCGACCTGACTGcc                   <  2:1232589/50‑1 (MQ=255)
                              cGTGAAACGCGGTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTAt              >  1:1880293/1‑50 (MQ=14)
                              cGTGAAACGCAGTATGAGTACAACGCCGCAGGCGACCTGACTGCCGTTAt              >  2:152869/1‑50 (MQ=255)
                                tGAAACGCGGTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCa            <  1:454543/50‑1 (MQ=14)
                                tGAAACGCGGTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCa            <  2:1262126/50‑1 (MQ=14)
                                tGAAACGCAGTATGAGTACATCGCCGCAGGCGACCTGACTGCCGTTATCa            <  1:1401637/50‑1 (MQ=255)
                                  aaaCGCGGTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCAcc          <  1:132596/50‑1 (MQ=14)
                                     cgcgGTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCg       <  1:759972/50‑1 (MQ=14)
                                      gCAGTATGAGTACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCg       >  1:960689/1‑49 (MQ=255)
                                          tATGAGTACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACgg  >  1:1153113/1‑50 (MQ=255)
                                               |                                            
CGGTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCAGTATGAGTACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGG  >  NC_017635/1637710‑1637801

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