Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176351,364,7860AT21.1% 83.0 / 12.2 38V8E (GTA→GAA) chaAsodium‑potassium/proton antiporter ChaA
Reads supporting (aligned to +/- strand):  ref base A (5/25);  new base T (8/0);  total (13/25)
Fisher's exact test for biased strand distribution p-value = 2.63e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATT  >  NC_017635/1364739‑1364834
                                               |                                                
gcgcCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTAcc                                                <  2:631575/50‑1 (MQ=255)
gcgcCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTcc                                                >  1:1996868/1‑50 (MQ=255)
 cgcCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCg                                               <  1:372752/50‑1 (MQ=255)
 cgcCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCg                                               <  2:1452521/50‑1 (MQ=255)
    cAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCct                                            <  1:1247128/50‑1 (MQ=255)
    cAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGCct                                            >  1:1345676/1‑50 (MQ=255)
     aGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCcccccc                                           >  1:1740679/1‑48 (MQ=255)
        aCCGGGAAGATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTg                                        <  1:708860/50‑1 (MQ=255)
           gggAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGc                                     <  1:1544268/50‑1 (MQ=255)
             gAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCCCCTTTAGCAt                                   >  1:1681895/1‑50 (MQ=255)
             gAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCCCCCTGAGCaa                                   >  1:79025/1‑49 (MQ=255)
              aaaaTAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCAtt                                  <  2:306253/50‑1 (MQ=255)
               aaaTAAGGGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGAAttt                                 <  1:165541/50‑1 (MQ=39)
               aaaTAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCAttt                                 <  2:786248/50‑1 (MQ=255)
               aaaTAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCAttt                                 <  1:325776/50‑1 (MQ=255)
               aaaTAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCAttt                                 <  2:1191040/50‑1 (MQ=255)
               aaaTAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCAttt                                 <  2:1554681/50‑1 (MQ=255)
                aaTAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCCCCTTTTGCAtttt                                >  1:1914896/1‑49 (MQ=255)
                 aTAAGCGAAGTCTCCTTGTGGCGGGTTTTTTTCCCCCCTTGAGCAtttta                               >  1:457899/1‑48 (MQ=255)
                 aTAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCCCTTTAGCATTTGa                               >  1:326701/1‑50 (MQ=255)
                  tAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGAc                              <  1:1376672/50‑1 (MQ=255)
                    aGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACa                             <  1:1770797/49‑1 (MQ=255)
                    aGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACAt                            <  2:974593/50‑1 (MQ=255)
                     gCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACAtt                           <  2:1507414/50‑1 (MQ=255)
                     gCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACAtt                           <  2:29213/50‑1 (MQ=255)
                        aaGTCTCCTTTTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATg                        >  1:320214/1‑50 (MQ=255)
                           tctcCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGtt                     <  1:1394366/50‑1 (MQ=255)
                                   tgGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCttt             <  2:1822701/50‑1 (MQ=255)
                                   tgGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCttt             >  1:1046247/1‑50 (MQ=255)
                                    ggCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTg            >  1:802649/1‑50 (MQ=255)
                                     gCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTg            <  1:809632/49‑1 (MQ=255)
                                       gggTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTg            >  1:2079208/1‑47 (MQ=255)
                                       gggTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAg         <  1:2039864/50‑1 (MQ=255)
                                        ggTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGa        <  2:184005/50‑1 (MQ=255)
                                            tttACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGa     <  2:1781697/49‑1 (MQ=255)
                                            tttACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGa     <  2:398460/49‑1 (MQ=255)
                                              tACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAAtt  >  2:1200070/1‑50 (MQ=255)
                                              tACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAAt   <  1:1047920/49‑1 (MQ=255)
                                               aCCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAAtt  <  1:2007207/49‑1 (MQ=255)
                                               |                                                
GCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTATGGTTATCCCTTTGCAGATGAATT  >  NC_017635/1364739‑1364834

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