Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,983,7000TG34.2% 64.8 / 23.1 41T397P (ACC→CCC) ECB_03713conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/20);  new base G (14/0);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 4.40e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTG  >  REL606/3983630‑3983774
                                                                      |                                                                           
gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGcccc                                                                        >  1:1746678/1‑74 (MQ=255)
gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTc                                                                       >  1:1679771/1‑76 (MQ=255)
gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTc                                                                       >  1:1294471/1‑76 (MQ=255)
  aaaTGCGGGGTTACTACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTccc                                                                     >  2:118961/1‑76 (MQ=255)
  aaaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTccc                                                                     <  1:158576/76‑1 (MQ=255)
   aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCca                                                                    <  1:1616894/76‑1 (MQ=255)
       cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACaaa                                                                <  1:152113/76‑1 (MQ=255)
       cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCACcaa                                                                >  1:1838118/1‑76 (MQ=255)
        ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCACAAAt                                                               >  2:101194/1‑76 (MQ=255)
         gggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCaaaat                                                              >  2:86489/1‑74 (MQ=255)
            ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCAAAATTGa                                                           >  2:939697/1‑76 (MQ=255)
            ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAtttta                                                           >  2:1408138/1‑74 (MQ=255)
             tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCCAAAtttcc                                                           >  2:306402/1‑74 (MQ=255)
              accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACg                                                         <  1:1436645/76‑1 (MQ=255)
                 accGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCAAAATTTACGATc                                                      >  2:783808/1‑76 (MQ=255)
                 accGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAATTTACGATc                                                      >  2:1132448/1‑76 (MQ=255)
                    gAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCt                                                   <  2:671269/76‑1 (MQ=255)
                     aGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTg                                                  <  1:1432926/76‑1 (MQ=255)
                           aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATg                                            <  1:1172467/76‑1 (MQ=255)
                           aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAAATTTACGATCGCTGGGGATg                                            >  1:422500/1‑76 (MQ=255)
                            ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGa                                           <  2:724065/76‑1 (MQ=255)
                                ttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCCAAATTTTCGATCGCTGGGGATGAAGcc                                       >  1:1144405/1‑76 (MQ=255)
                                 tgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAAATTGACGATCGCCGGGGATGACGCCa                                      >  1:109303/1‑76 (MQ=255)
                                        tACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAATTTACGATCGCCGGGGAAGACGCCATcgcgga                               >  1:1470000/1‑74 (MQ=255)
                                        tACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAAATTTACGATCGCTGGGGATGACGCCCTcgcgga                               >  1:1043619/1‑74 (MQ=255)
                                        tACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCCCAAAATTTTCGGTCGCTGGGGGTGACCCCATtgcgga                               >  1:2538/1‑74 (MQ=255)
                                           ggCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcc                            <  2:570162/76‑1 (MQ=255)
                                           ggCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcc                            <  2:116294/76‑1 (MQ=255)
                                              tGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAATTTACGATCGCTGGAGATGAAGCCATCGCGGAGccacc                         >  2:766235/1‑76 (MQ=255)
                                                ttATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccccccc                       >  2:871422/1‑75 (MQ=255)
                                                   ttGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCacac                    <  2:685699/76‑1 (MQ=255)
                                                   ttGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCacac                    <  1:1452134/76‑1 (MQ=255)
                                                    tggtggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACAcc                   >  1:517213/1‑76 (MQ=255)
                                                     ggtggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACAccc                  >  1:1293207/1‑75 (MQ=255)
                                                       tggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGc                <  1:1247134/76‑1 (MQ=255)
                                                       tggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGc                <  1:1459083/76‑1 (MQ=255)
                                                        ggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCg               <  1:118961/76‑1 (MQ=255)
                                                              ggCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa         <  2:1504957/76‑1 (MQ=255)
                                                              ggCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa         <  1:71590/76‑1 (MQ=255)
                                                                cAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGt       >  2:104083/1‑76 (MQ=255)
                                                                 aGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa      >  1:712893/1‑76 (MQ=255)
                                                                  ggggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTaa     <  1:871422/76‑1 (MQ=255)
                                                                    ggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt   <  1:306402/76‑1 (MQ=255)
                                                                    ggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt   >  1:1321702/1‑76 (MQ=255)
                                                                    ggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt   <  1:783808/76‑1 (MQ=255)
                                                                    ggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCAGCGGTTTCAGTAAtt   >  1:1116147/1‑76 (MQ=255)
                                                                     gTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTg  >  1:860275/1‑76 (MQ=255)
                                                                      |                                                                           
GTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTG  >  REL606/3983630‑3983774

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