Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,983,7000TG37.9% 40.0 / 14.0 29T397P (ACC→CCC) ECB_03713conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/15);  new base G (10/1);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 2.61e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTG  >  REL606/3983629‑3983774
                                                                       |                                                                          
tgtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCt                                                                        <  2:547834/76‑1 (MQ=255)
    aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCca                                                                    <  2:909078/76‑1 (MQ=255)
    aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGccccccc                                                                    >  1:578127/1‑75 (MQ=255)
       gCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACaa                                                                 <  1:1403314/76‑1 (MQ=255)
       gCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACaa                                                                 <  1:1587080/76‑1 (MQ=255)
         ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCaaaa                                                               >  1:118198/1‑75 (MQ=255)
            gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCCCCAAATTg                                                            >  1:610584/1‑76 (MQ=255)
            gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCACAATTTg                                                            >  1:222351/1‑76 (MQ=255)
            gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGACTCCCCCAAAtttt                                                            >  2:1689803/1‑75 (MQ=255)
             ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCAAAattta                                                           >  1:1659151/1‑74 (MQ=255)
                 ccccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGAt                                                       <  1:1212541/74‑1 (MQ=255)
                 caccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCACAAAATTGACGAt                                                       >  2:1021484/1‑76 (MQ=255)
                     gAGGGTGATTGCTTGCATCCTACGGCTGTTATTTGTGGCATCGGCAGGGCTGCCTCCCACAAATTTGACGATCGCt                                                   <  2:1390836/76‑1 (MQ=255)
                     gAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAATTTACGATCGcc                                                   >  2:1370947/1‑75 (MQ=255)
                       gggTGATTGCTCGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCACACATTTGACGATCGCTgg                                                 >  2:515627/1‑76 (MQ=255)
                        ggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGga                                                <  1:480535/76‑1 (MQ=255)
                         gTGATTGCTTGCATCCTACGGCNGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAATTTGACGATCGCTgggg                                               >  2:237083/1‑74 (MQ=255)
                           gATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAATTGACGATCGCTGGAGAt                                             >  1:1653363/1‑76 (MQ=255)
                                    cATCCTACGGCTGTTATTGGTGGCATCGGGAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATc                                    <  2:289691/76‑1 (MQ=255)
                                     aTCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAAATTGACGACCGCTGGAGATGCAGCCacgg                                   >  2:391080/1‑73 (MQ=255)
                                       ccTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcgcg                                 <  1:608779/76‑1 (MQ=255)
                                       ccTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCAAAATTTACGATCGCTGGGGATGACGCCATcgcg                                 >  1:692703/1‑76 (MQ=255)
                                        cTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCAAATTTTTCGATCGCTGGGGATGACGCCATcgcgg                                >  1:118214/1‑75 (MQ=255)
                                          aCGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAg                              <  1:1332414/76‑1 (MQ=255)
                                            ggCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcc                            <  1:1552032/76‑1 (MQ=255)
                                            ggCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCCCAAATTTGACGCTCGCCGGGGCTGCCGCCCTCCCGAAGcc                            >  1:841564/1‑76 (MQ=255)
                                             gCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcca                           <  2:1503737/76‑1 (MQ=255)
                                             gCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcca                           <  2:217360/76‑1 (MQ=255)
                                                 ttATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTGGGGATGACGCCATCGCGAAGccaccca                       >  1:532732/1‑76 (MQ=255)
                                                     tggtggCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACAcc                   >  1:487824/1‑76 (MQ=255)
                                                              cGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGAGGTTTc          >  2:1011726/1‑76 (MQ=255)
                                                               ggCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa         <  1:1359306/76‑1 (MQ=255)
                                                                gCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAg        <  1:200597/76‑1 (MQ=255)
                                                                 cAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGt       <  2:925657/76‑1 (MQ=255)
                                                                  aGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa      >  2:469154/1‑76 (MQ=255)
                                                                      gTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTg  >  2:1125813/1‑76 (MQ=255)
                                                                       |                                                                          
TGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTG  >  REL606/3983629‑3983774

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