Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176352,321,6420GT27.8% 86.5 / 20.2 37A18A (GCC→GCA) ECW_RS11635glycoside hydrolase family 104 protein
Reads supporting (aligned to +/- strand):  ref base G (13/13);  new base T (8/2);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-02
Rejected as consensus: Frequency below/above cutoff threshold.

CATAGCCCCGGTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGATGTGTGCTAATTAC  >  NC_017635/2321593‑2321689
                                                 |                                               
cATAGCCCCGGTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACg                                                 >  2:518154/1‑50 (MQ=32)
 aTAGCCCCGGTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACg                                                 <  2:1239535/49‑1 (MQ=32)
    gCCCCGGTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACGGcc                                              <  1:502628/49‑1 (MQ=40)
      cccGGTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACGGCCAg                                            >  2:315793/1‑49 (MQ=255)
       ccGGTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACTGCCAGCa                                          >  1:2079101/1‑50 (MQ=255)
        cGGTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACGGCCAGCAt                                         <  1:113608/50‑1 (MQ=255)
        cGGTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACGGCCAGCAt                                         >  2:1621192/1‑50 (MQ=255)
        cGGTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACGGCCAGCAt                                         >  1:1390957/1‑50 (MQ=255)
         ggTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACTGCCAGCATg                                        >  2:53317/1‑50 (MQ=255)
         ggTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACGGCCAGCATg                                        <  2:845851/50‑1 (MQ=255)
         ggTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACGGCCAGCATg                                        >  2:1048429/1‑50 (MQ=255)
         ggTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACCGCCAGCATg                                        >  1:1741280/1‑50 (MQ=255)
            ttttCGTCAGCGGATGATTCGCCGTCCCTTCGGACACTGCCAGCATGTc                                      >  2:540526/1‑49 (MQ=255)
             tttCGTCAGCGGATGATTCGCCGTCCCTTCGGACACTGCCAGCATGTc                                      >  1:688452/1‑48 (MQ=255)
                cGTCAGCGGATGATTCGCCGTCCCTTCGGACACTGCCAGCATGTCgaga                                  >  1:1633357/1‑49 (MQ=255)
                  tCAGCGGATGATTCGCCGTCCCTTCGGACACTGCCAGCATGTCGAGAAAg                               >  1:1221858/1‑50 (MQ=255)
                  tCAGCGGATGATTCGCCGTCCCTTCGGACACGGCCAGCATGTCGAGAAAg                               >  1:67648/1‑50 (MQ=255)
                       ggATGATTCGCCGTCCCTTCGGACACTGCCAGCATGTCGAGAAAggcggc                          >  1:159314/1‑50 (MQ=255)
                       ggATGATTCGCCGTCCCTTCGGACACGGCCAGCATGTCGAGAAAggcggc                          >  2:1158259/1‑50 (MQ=255)
                        gatgatTCGCCGTCCCTTCGGACACGGCCAGCATGTCGAGAAAggcggc                          >  2:485198/1‑49 (MQ=255)
                        gatgatTCGCCGTCCCTTCGGACACGGCCAGCATGTCGAGAAAggcggc                          >  2:1799577/1‑49 (MQ=255)
                          tgatTCGCCGTCCCTTCGGACACTGCCAGCATGTCGAGAAAGGCGGCGa                        >  2:991413/1‑49 (MQ=255)
                            atTCGCCGTCCCTTCGGACACGGCCAGCATGTCGAGAAAGGCGGCGatat                     <  1:514574/50‑1 (MQ=255)
                             ttCGCCGTCCCTTCGGACACTGCCAGCATGTCGAGAAAGGCGGCGATAtt                    <  1:647980/50‑1 (MQ=255)
                                gCCGTCCCTTCGGACACTGCCAGCATGTCGAGAAAGGCGGCGATATTCt                  <  2:1816252/49‑1 (MQ=255)
                                gCCGTCCCTTCGGACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTg                 <  2:431737/50‑1 (MQ=255)
                                  cGTCCCTTCGGACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGAt               <  2:1397804/50‑1 (MQ=255)
                                  cGTCCCTTCGGACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGAt               >  2:57555/1‑50 (MQ=255)
                                   gTCCCTTCGGACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGAtg              <  2:621279/50‑1 (MQ=255)
                                   gTCCCTTCGGACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGAtg              >  2:1555979/1‑50 (MQ=255)
                                    tCCCTTCGGACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGAtgt             >  1:684718/1‑50 (MQ=255)
                                    tCCCTTCGGACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGAtgt             >  1:2004382/1‑50 (MQ=255)
                                       cTTCGGACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGAtgtgtg          <  1:1184834/50‑1 (MQ=255)
                                          cGGACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGATGTGTGCt        <  1:1770979/49‑1 (MQ=255)
                                           ggACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGATGTGTGCTaa      <  1:637963/50‑1 (MQ=255)
                                            gACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGATGTGTGCTAAt     <  1:2040843/50‑1 (MQ=255)
                                            gACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGATGTGTGCTAAt     <  1:583474/50‑1 (MQ=255)
                                              cacGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGATGTGTGCTAATTa   <  1:2041086/50‑1 (MQ=255)
                                               acGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGATGTGTGCTAATTAc  >  1:810229/1‑50 (MQ=255)
                                                 |                                               
CATAGCCCCGGTTTTTCGTCAGCGGATGATTCGCCGTCCCTTCGGACACGGCCAGCATGTCGAGAAAGGCGGCGATATTCTGATGTGTGCTAATTAC  >  NC_017635/2321593‑2321689

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