Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606430,3180TC46.4% ‑6.5 / 13.5 28D157G (GAC→GGC) insL‑2putative transposase insL for insertion sequence IS186
Reads supporting (aligned to +/- strand):  ref base T (2/13);  new base C (0/13);  total (2/26)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCGGGACGCGAACCGAATCCCCGGTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAGGATCATAT  >  REL606/430253‑430386
                                                                 |                                                                    
atgcGGGCCGCAACCGGAACCCCCGGTCACCAACGCTTATCTCGTCGGCCGTTTGCGCAAATCGGTCCAGCCGTTc                                                            <  2:2798037/73‑1 (MQ=0)
atgcGGGACGCGACCCGAACCCCCGGTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTc                                                            <  1:5967009/73‑1 (MQ=1)
atccGGGCCGCGACCCGAATCCCCGGTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTc                                                            <  1:4201993/73‑1 (MQ=1)
 ttggggCCGCGACCCGAACCCCCGGTCAGCAAGGCTTATCTCGTCTGCCGTTTGCGCAAATCGGCCCAGCCGTTCa                                                           <  2:3020216/73‑1 (MQ=1)
              cctaaTCCCCGGTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGt                                              >  1:6636734/4‑76 (MQ=1)
                 aaccccgGGTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCgg                                           <  2:3007722/73‑1 (MQ=1)
                 aacccccGGTCAGGAGTGCGGCGCGCGTCTGCGGTTTGCCAAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCgg                                           <  2:6138692/73‑1 (MQ=0)
                 aacccccGGTCAGCAATGCGTATCTCGCCTGCCTTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCgg                                           <  2:2538116/73‑1 (MQ=1)
                 aacccccGGTCAGCAATCCTTTCCCCGTCGGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCgg                                           <  1:6944203/73‑1 (MQ=1)
                 aacccccGGTCAGAAATGCTTATCTCGTCTGCGGTTTGCGAAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCgg                                           <  1:6501584/73‑1 (MQ=1)
                 aacccccGGTCAGAAAGGCGTATCTCCTCTGCCGTTTGCCCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCgg                                           <  1:5868493/73‑1 (MQ=1)
                     ccgggTCAGCATTGCGTATCTCGTCTGCCGTTGGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAg                                       <  1:3724628/73‑1 (MQ=1)
                     ccgggTCAGCATGGCGCACCCCGTCTGCCGTTTGCACAAACCTGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAg                                       <  1:691505/73‑1 (MQ=1)
                     ccgggTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAg                                       <  2:282275/73‑1 (MQ=1)
                     ccgggTCAGAAAGCCGTATCCCGTCTGCCGTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAg                                       <  1:156221/73‑1 (MQ=1)
                        gggcAGCATTGCGTATCTCGCCTGCCGTTTGCGCAAACCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTc                                    <  1:4020235/73‑1 (MQ=1)
                        gggcAGAAATGCGTATCTCGTCTGCCGTTGGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTc                                    <  1:1839317/73‑1 (MQ=1)
                        ggccAGAAATGCTTATCTCGTCTGCCTTTGGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTc                                    <  1:8527003/73‑1 (MQ=1)
                            agaaaTGCGTTCCTCGTCTGCCGTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCaaaa                                <  1:9225508/73‑1 (MQ=1)
                              cattGCTTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATc                              <  1:9478431/73‑1 (MQ=0)
                                 tgagTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGt                           >  2:7860471/4‑76 (MQ=1)
                                  gcttATCTCGTCTGCCTTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTg                          <  1:5334096/73‑1 (MQ=1)
                                  gcttATCTCGTCTGCCGTTTGGGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTg                          <  1:4126891/73‑1 (MQ=0)
                                  gcttATCTCGCCTGCCGTTTGCGCAAACCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTg                          <  1:3128762/73‑1 (MQ=1)
                                  gcttATCACGTCTCCCGTTTGCCCAACTCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTg                          <  1:8343155/73‑1 (MQ=1)
                                    gtttctcGTTTGCGGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGaa                        <  1:150571/73‑1 (MQ=1)
                                    gtttctcGCCTGCCGTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGaa                        <  1:6054632/73‑1 (MQ=1)
                                       tcgcGTCTTCCGTTGGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTg                     <  1:8097728/73‑1 (MQ=1)
                                       tcgcGTCTGCGTTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTg                     <  1:5765986/73‑1 (MQ=1)
                                       tcgcGTCTGCCGTTTGCGCAAACCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTg                     <  1:4993036/73‑1 (MQ=1)
                                       tcccGTGTGCCGTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTg                     <  2:6842677/73‑1 (MQ=1)
                                       tcccGTCTGCCTTTTGCGCAAACCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTg                     <  1:2310301/73‑1 (MQ=1)
                                       tcccGTCTGCCGTTGGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTg                     <  1:7504321/73‑1 (MQ=1)
                                          cgccTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACa                  <  2:6125758/73‑1 (MQ=1)
                                          cgccTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACa                  <  2:10052384/73‑1 (MQ=1)
                                          cgccTGCCGTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACa                  <  1:6044593/73‑1 (MQ=0)
                                          cgccTGCCGTTGGCGCAAACCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACa                  <  2:1929798/73‑1 (MQ=1)
                                          cgccTGCCGTTGGCCCAAAGCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACa                  <  1:522666/73‑1 (MQ=1)
                                             ctcccGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGt               <  2:9641367/73‑1 (MQ=1)
                                             ctcccGTTTGCCCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGt               <  2:7143998/73‑1 (MQ=0)
                                              tgacGTGTGCGCAAATGGGCCCAGCCGTGCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTa              <  2:868796/73‑1 (MQ=1)
                                               gcggTTTGCACAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGGACTGACAGGTAt             <  1:5523995/73‑1 (MQ=1)
                                               gcggTTGGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTAt             <  1:6932686/73‑1 (MQ=1)
                                               gcggTTGGCGCAAACCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTAt             <  1:3125164/73‑1 (MQ=1)
                                               gcagTTCGCGCAACGCGTCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTAt             <  2:3577442/73‑1 (MQ=1)
                                                 cggttGCGCAAAGCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGa           <  1:7414085/73‑1 (MQ=0)
                                                   ttggcgcAATCCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAgg         <  1:3612833/73‑1 (MQ=1)
                                                   ttggcgcAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAgg         <  1:475545/73‑1 (MQ=0)
                                                   ttggcgcAAACCGCTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAgg         <  2:5539482/73‑1 (MQ=1)
                                                   ttggcgaAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAgg         <  1:8419232/73‑1 (MQ=1)
                                                   ttggcgaAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAgg         <  1:5269696/73‑1 (MQ=1)
                                                   ttcgcgcAAGTCGTTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAgg         <  1:5404301/73‑1 (MQ=1)
                                                   ttcgcgcAAGTCGGTCCAGCCGTTCAGCGTCTCTGCTGGCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAgg         <  2:6836850/73‑1 (MQ=1)
                                                          aactCGCTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAGGATCatat  <  2:6915874/73‑1 (MQ=0)
                                                                 |                                                                    
ATTCGGGACGCGAACCGAATCCCCGGTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAGGATCATAT  >  REL606/430253‑430386

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