Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,005,5610CT39.6% 35.5 / 58.9 48S146S (TCC→TCT) ECB_02802conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/11);  new base T (10/9);  total (28/20)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACGGCCATCGTGTCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGCAGAGGCTTTT  >  REL606/3005488‑3005635
                                                                         |                                                                          
aCGGCCATCGTGTCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCgg                                                                          <  1:242479/76‑1 (MQ=255)
   ttcATCGTGTCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCTGGaaa                                                                       <  1:930305/74‑1 (MQ=255)
      aTCGTGTCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCTGGaaaaan                                                                    <  1:1366061/76‑2 (MQ=255)
      aTCGTGTCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGaaaaaa                                                                    >  2:119282/1‑76 (MQ=255)
        cGTGTCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACg                                                                  <  1:616732/76‑1 (MQ=255)
        cGTGTCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGCTTCCGGAAAAAACg                                                                  >  1:923114/1‑76 (MQ=255)
          tgtCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCTGGAAAAAACGtt                                                                <  1:47104/76‑1 (MQ=255)
            tCAGTCTGCTGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTaa                                                              <  1:410686/76‑1 (MQ=255)
            tCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCTGGAAAAAACGTTaa                                                              <  2:1244275/76‑1 (MQ=255)
            tCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTaa                                                              <  1:337103/76‑1 (MQ=255)
            tCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCCTTTCTGGAAAAAACGATaa                                                              >  1:874133/1‑76 (MQ=255)
             cAGTCTGTCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCTGGAAAAAACGTTAAc                                                             >  1:1142185/1‑76 (MQ=255)
             cAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAAc                                                             >  2:793394/1‑76 (MQ=255)
             cAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAAc                                                             >  2:793387/1‑76 (MQ=255)
              aGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCCTTTCCGGAAAAAACCTTAAcc                                                            >  1:830190/1‑75 (MQ=255)
                tCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCTGGAAAAAACGTTAACTGa                                                          >  1:239141/1‑76 (MQ=255)
                tCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGa                                                          <  2:1596201/76‑1 (MQ=255)
                 cTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCTGGAAAAAACGTTAACTGAt                                                         >  1:408971/1‑76 (MQ=255)
                  tGCCGTTATCGGTTTTCGACCCCCAGGACTACCATCCACTGCTAATTACCGTTTCTGGAAAAAACGATAACTCATc                                                        >  1:1506301/1‑76 (MQ=255)
                      gTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCTGGAAAAAACGTTAACTGATCCCTg                                                    >  1:13471/1‑76 (MQ=255)
                       ttATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGa                                                   >  2:364648/1‑76 (MQ=255)
                       ttATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGa                                                   >  2:1173959/1‑76 (MQ=255)
                          tCGGTTTGCGATCCCCAGGTCTCCCATCCACTGCTAATTTCCGTTTCTGGAAAAAACGTTAACTGATCCCTGAAGt                                                <  2:1090331/76‑1 (MQ=255)
                           cGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTc                                               >  2:153873/1‑76 (MQ=255)
                                  cGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCa                                        <  1:1738068/76‑1 (MQ=255)
                                     tCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGc                                     >  1:312477/1‑76 (MQ=255)
                                       cccAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCtt                                   >  2:1062572/1‑76 (MQ=255)
                                         cAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCa                                 >  1:357763/1‑76 (MQ=255)
                                            gACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGcc                              >  1:618129/1‑76 (MQ=255)
                                               tACCATCCACTGCTAATTACCGTTTCTGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCctct                           >  1:610420/1‑76 (MQ=255)
                                               tACCATCCACTGCTAATTACCGTTTCTGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCctct                           <  2:275715/76‑1 (MQ=255)
                                                 ccATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGt                         >  1:681619/1‑76 (MQ=255)
                                                      cACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCGGCCTCTGTCACaa                    >  2:1320968/1‑76 (MQ=255)
                                                        cTGCTAATTACCGTTTCTGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAAcg                  >  1:1236350/1‑76 (MQ=255)
                                                        cTGCTAATTACCGTTTCTGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAAcg                  >  1:1430571/1‑76 (MQ=255)
                                                        cTGCTAATTACCGTTTCTGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAAcg                  <  2:662039/76‑1 (MQ=255)
                                                        cTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAAcg                  <  1:19747/76‑1 (MQ=255)
                                                         tGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAAcgc                 >  1:853103/1‑76 (MQ=255)
                                                           cTAATTACCGTTTCTGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGt               >  2:928939/1‑76 (MQ=255)
                                                           cTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACCACGCGt               >  1:613503/1‑76 (MQ=255)
                                                           cTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGt               >  1:532513/1‑76 (MQ=255)
                                                            tAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTg              <  2:797016/76‑1 (MQ=255)
                                                              tttACCGTTTCTGGAGAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAGCGCGTGGc            <  2:1506301/75‑1 (MQ=255)
                                                              aTTACCGTTTCTGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGc            >  1:1119484/1‑76 (MQ=255)
                                                              aTTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGc            <  2:195126/76‑1 (MQ=255)
                                                               ttACCGTTTCTGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGCa           <  2:906766/76‑1 (MQ=255)
                                                               ttACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGCa           >  1:1459798/1‑76 (MQ=255)
                                                                tACCGTTTCTGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGCag          >  2:691725/1‑76 (MQ=255)
                                                                tACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGCag          <  2:125324/76‑1 (MQ=255)
                                                                     tttCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGCAGAGGCt     >  1:298714/1‑76 (MQ=255)
                                                                     tttCCGAAAAAAACGTGAACTGAGCCGTGAAGTCAGTTTCATGCTTCGGCCTCTGTCACAACGCGTGGCAGAGGCt     <  2:1406383/76‑1 (MQ=255)
                                                                      ttCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGCAGAGGCtt    <  1:1685844/76‑1 (MQ=255)
                                                                      ttCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTGCAGCCTCTGTCACAACGCGTGGCAGAGGCtt    <  2:1619826/76‑1 (MQ=255)
                                                                        ccGGACAGAACGTTAACTGATCCCTGAAGTCAGTTCTATGCTTCAGCCTCTGTCACAACGCGTGGCAGAGGCtttt  >  2:1581381/1‑76 (MQ=255)
                                                                        ccGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGCAGAGGCtttt  <  2:745636/76‑1 (MQ=255)
                                                                        ccGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGCAGAGGCtttt  <  2:427248/76‑1 (MQ=255)
                                                                        ccGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGCAGAGGCtttt  >  1:1204529/1‑76 (MQ=255)
                                                                         |                                                                          
ACGGCCATCGTGTCAGTCTGCCGTTATCGGTTTTCGATCCCCAGGACTACCATCCACTGCTAATTACCGTTTCCGGAAAAAACGTTAACTGATCCCTGAAGTCAGTTTCATGCTTCAGCCTCTGTCACAACGCGTGGCAGAGGCTTTT  >  REL606/3005488‑3005635

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