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

CGGATACATTCGGGACGCGAACCGAATCCCCGGTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAGGATCATATC  >  REL606/430246‑430387
                                                                         |                                                                     
cgtaTACATTCGGGACGCGAACCGAATCCCCGGTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCa                                                                    <  1:2110423/73‑1 (MQ=1)
       atgcGGGACGCGAACCGAACCCCCGGTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTc                                                             <  1:2768938/73‑1 (MQ=0)
       atgcGGGACGCGAACCAAATCCCCGGTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTc                                                             <  1:4310304/73‑1 (MQ=0)
        ttggggACGGGAACCGAATCCCGGGTCAGCAAGGCGTATCTCGTCTGCCGTTTGCGCAAATCGGCCCAGCCGTTCa                                                            <  1:4739813/73‑1 (MQ=1)
            ggccgcgAACCAAATCCCCGGTCAGCAATGTGTATCTCGCCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGt                                                        <  1:6338442/73‑1 (MQ=1)
                gcaaaCCGAACCCCCGGCAGCAAAGGCTTGTCCCGCCTGCCGTTGGCGCAAATCGGCCCAGCCGTTCAGCGTctct                                                    <  2:4815798/73‑1 (MQ=0)
                  gatccGAACCCCCGGTCAGCAGTGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTctgc                                                  <  1:4485819/73‑1 (MQ=1)
                  gacccGAACCCCGGGTCAGAAATGCGTTTCCCGTCTGCCGTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTctgc                                                  <  2:4025015/73‑1 (MQ=1)
                            ccgggTCAGCAATGCGTCTCTCGCCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAg                                        <  1:1208834/73‑1 (MQ=1)
                            ccgggTCAGAAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAg                                        <  2:3379794/73‑1 (MQ=1)
                            ccgggTAAGAAATGGTTACCTCGCCTGCCGTTTGCCCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAg                                        <  1:2140571/73‑1 (MQ=1)
                                  caccAATCCTTATCTCGTCTGCCGTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCaaa                                  <  1:3774979/73‑1 (MQ=1)
                                      cattGCGTTCCTCGCCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATc                               <  2:138238/73‑1 (MQ=1)
                                          gcttATCTCGTCTGCCGTTTGCGAAAACCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTg                           <  2:3777848/73‑1 (MQ=1)
                                            gtttctcGTCTGCCGTTTGCGCAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGaa                         <  2:5298404/73‑1 (MQ=0)
                                                  cgccTGCCGTTTGCACAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACa                   <  1:2728198/73‑1 (MQ=1)
                                                     ctcccGTTGCCGCAAACCGCCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGt                <  1:4000700/73‑1 (MQ=1)
                                                       gcggTTTGCGCACATGGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTAt              <  2:1510742/73‑1 (MQ=1)
                                                       gcggTTTGCGAAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTAt              <  2:2407733/73‑1 (MQ=0)
                                                        ccttttGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATg             <  2:5079287/73‑1 (MQ=1)
                                                        ccttttGCGCAAAGGGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATg             <  1:5958792/73‑1 (MQ=1)
                                                        ccttttGCGCAAACCGCTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATg             <  1:725229/73‑1 (MQ=1)
                                                          gtgtGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAg           <  1:1826192/73‑1 (MQ=1)
                                                              gcccAAATCGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAGGATc       <  1:2329582/73‑1 (MQ=0)
                                                               cgaaaaTCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAGGATCa      <  2:2248456/73‑1 (MQ=1)
                                                                   aagcGGCCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAGGATCATATc  <  1:4778871/73‑1 (MQ=0)
                                                                         |                                                                     
CGGATACATTCGGGACGCGAACCGAATCCCCGGTCAGCAATGCGTATCTCGTCTGCCGTTTGCGCAAATCGGTCCAGCCGTTCAGCGTCTCTGCTGTCGGTTAGCTCAAAATCAGTGAACTGACAGGTATGAGGATCATATC  >  REL606/430246‑430387

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