Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,775,1250AC53.3% ‑6.0 / 11.7 15V63G (GTA→GGA) insK‑3IS150 putative transposase
Reads supporting (aligned to +/- strand):  ref base A (1/6);  new base C (0/8);  total (1/14)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTTCCCATCAGGCGCTGAACAGCTTTATGGTTAATCTGTTTCCCTTCTCGATGAAGAGACAGCGTTACCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCTTTTTAACGTCCGC  >  REL606/2775057‑2775192
                                                                    |                                                                   
ggtttCCCATCAGGCGCTGAACAGCTTTATGGTTAATCTGTTTCCCTTCTCGATGAAGAGACAGCGTTACCCTACg                                                              <  1:5092410/73‑1 (MQ=1)
       cagcAGGCGCTGAACAGCTTTATGGTTAATCGGTTGCCCTTCTCGATGAAGAGGCAGCGTTACCCTACGGTATCCg                                                       <  2:2450770/73‑1 (MQ=1)
         tccggCGCTAAACAGCTTTATGGTTAATCTGTTTCCCTTCCCGATGAAGAACCGGCGTTCCCCTACGGTATCCGTa                                                     <  2:8118704/73‑1 (MQ=1)
           agtcgcTGAACAGTTTTATGGTTAATCTGTTTCCCTTCTCGATGAAGAGACAGCGTTCCCCTACGTTATCCGTATc                                                   <  2:1617471/73‑1 (MQ=1)
                  gagcAGCTTTATGTTTAATCTGTTTCCTTTCTCGATGAAGAGACAGCGTTACCCTACGGTATCCGTATCGGCctct                                            <  1:6665093/73‑1 (MQ=1)
                  gaccAGTTTTGTGGTTATTTTTGTCCCCTTCTCGGTGAGGAGACAGTGTTGCCCTCGGGTATCCGTATCGGCctct                                            <  2:7804906/73‑1 (MQ=0)
                  gaccAGCTTTATGGTTAATCTGTTCCCCTTCTTGATGAAGAGACAGCGTTCCCCTCCGGTATCCGTATCGGCctct                                            <  1:2157577/73‑1 (MQ=1)
                     catcTTTATGGTTAATCTTTTTCCCTTCTCGATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTAtt                                         <  2:5033686/73‑1 (MQ=1)
                     catcTTTATGGTTAATCTGTTTCCCTTCTCGATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTAtt                                         <  1:4712594/73‑1 (MQ=0)
                     catcTTTATGGTTAATCTGTTTCCCTTCTCGATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTAtt                                         <  2:291679/73‑1 (MQ=0)
                     catcTTTATGGTTAATCTGTTTCCCTTCTCGATAAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTAtt                                         <  2:4645661/73‑1 (MQ=1)
                     catcTTTAGGGTTAATCTGTTTCCCTTCTCGATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTAtt                                         <  2:6638196/73‑1 (MQ=1)
                     catcTTTAGGGTTAATCTCTTGCCCTTCACGATGAGGAGACAGCGTTACCCTACGGTATCCGTATCGGCCTCTAtt                                         <  1:6529114/73‑1 (MQ=1)
                     caccTTTATGGTTAATCTGTTTCCCTCCTCGATGAAGAGACAGCGTTCCCCTCCGGTATCCGTATCGGCCTCTAtt                                         <  1:2160932/73‑1 (MQ=1)
                      agttttATGTTTAATCTGTTTCCCTTCTCGATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATtc                                        <  2:6149955/73‑1 (MQ=1)
                      agttttAGGTTTATTCTTTTTCCCTTCCCGAAGAAAAGACAGCGTTACCCTACGGTATCCGTATCGGCCTCTATtc                                        <  2:6630574/73‑1 (MQ=1)
                          ttttGGTTAATCTGTTCCCCTTCTCGATAAAGAGACAGCGTTCCCCTACGTTATCCGTATCGGCCTCTATTCTCGt                                    <  2:865235/73‑1 (MQ=1)
                           tagggTTAATCTTTTTCCCTCCTCGATGAAGAGCCAGCGTTCCCCTGCGTTATCCGTATCGGCCTCTATTCTCGTg                                   <  2:312719/73‑1 (MQ=1)
                           tagggTTAATCTGTTTCCCTTCTGGATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTg                                   <  2:2176804/73‑1 (MQ=1)
                           tagggTTAATCTGTTTCCCTTCTCGATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTg                                   <  1:1642759/73‑1 (MQ=0)
                           tagggTTAATCTGTTTCCCTTCTCGAGGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTg                                   <  2:1989785/73‑1 (MQ=1)
                             tgtttACTCTGTTCCCCTTCCCGATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGAt                                 <  1:6255839/73‑1 (MQ=1)
                             tgtttAATCTTTTTCCCTTCTCGATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGAt                                 <  1:4080527/73‑1 (MQ=1)
                             tgtttAATCTTTTTCCCTTCTCGATGAAGAGACAGCGTTACCCTACGGTATCCGTATCGGCCTCTATTCTCGTGAt                                 <  1:6242189/73‑1 (MQ=0)
                                     ctttttCCCTTCTCGATGAAGAGACAGCTTTACCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCa                         <  2:4071822/73‑1 (MQ=0)
                                     ctttttCCCTTCTCGATGAAGAGACAGCGTTACCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCa                         <  1:7662706/73‑1 (MQ=1)
                                        ttccccTTCTGGATGAAGAGACAGCTTTCCCCTCCGTTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTa                      <  2:2846573/73‑1 (MQ=1)
                                        ttccccTTCTCGATGAGGAGACAGCGTTACCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTa                      <  2:56504/73‑1 (MQ=0)
                                        ttccccTTCTCGATGAAGAGACAGCGTTCCCCTACGTTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTa                      <  1:1071207/73‑1 (MQ=1)
                                        ttccccTTCTCGATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTa                      <  2:878227/73‑1 (MQ=0)
                                        ttccccTTCTCGATGAAGAGACAGCGTTACCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTa                      <  1:2867862/73‑1 (MQ=1)
                                           cctttCTCGGTAAAGAAACAGCTTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATa                   <  1:7926430/73‑1 (MQ=1)
                                           cctttCTCGATGAAGAGCCAGCGTTACCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATa                   <  2:6782645/73‑1 (MQ=0)
                                           cctttCTCGATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATa                   <  2:4205184/73‑1 (MQ=0)
                                              ttatcGATGAGGACACGCCGTATCCCGACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGc                <  1:7647282/73‑1 (MQ=1)
                                                ctggATGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCtt              <  2:2349606/73‑1 (MQ=0)
                                                ctagATGAAGAGACAGCGTTACCCCACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTaa                     <  2:2793101/66‑1 (MQ=0)
                                                 tcaaTGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCttt             <  2:6379708/73‑1 (MQ=0)
                                                  cgctGAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCtttt            <  1:6124923/73‑1 (MQ=0)
                                                   gaggAAGAGACAGCTTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCttttt           <  1:6501085/73‑1 (MQ=1)
                                                   gaggAAGAGACAGCTTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCttttt           <  1:2701110/73‑1 (MQ=1)
                                                   gaggAAGAGACAGCGTTCCCCTAGGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCttttt           <  2:7486320/73‑1 (MQ=1)
                                                   gaggAAGAGACAGCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCttttt           <  2:4324955/73‑1 (MQ=0)
                                                   gaggAAGAGACACTGTTCCCCTACGTTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCttttt           <  1:1386113/73‑1 (MQ=1)
                                                    ataaaGAGACGGCTTTCCCCTACGTTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCTTTTTa          <  2:6629964/73‑1 (MQ=1)
                                                      gaggagaCCGAGTAACCGTAGGTTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCTTTTTAAc        <  2:561936/73‑1 (MQ=1)
                                                      gaggagaCAGCGTTACCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATTTCACTAATACGCTTTTTAAc        >  1:623117/4‑76 (MQ=0)
                                                            acggCGTTCCCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCTTTTTAACGTCCGc  <  1:3039839/73‑1 (MQ=0)
                                                                    |                                                                   
GGGTTCCCATCAGGCGCTGAACAGCTTTATGGTTAATCTGTTTCCCTTCTCGATGAAGAGACAGCGTTACCCTACGGTATCCGTATCGGCCTCTATTCTCGTGATAAATCTCACTAATACGCTTTTTAACGTCCGC  >  REL606/2775057‑2775192

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