Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,223,0290AT52.9% ‑6.5 / 12.3 24L439Q (CTG→CAG) PP_4603transposase
Reads supporting (aligned to +/- strand):  ref base A (0/8);  new base T (0/9);  total (0/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAAC  >  NC_002947/5222945‑5223166
                                                                                     |                                                                                                                                         
gctactgcgcAGCCCGGATGGGCTGGGTGATCTCACACCAAACAAAATGAAACCCTTTGTTTTCTAAGCACTCAAGCGGCCACCATTTGATACCCCTGTTTCGGGGATTTGCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGtt                                                                          <  2:3695839/144‑1 (MQ=1)
             ccggtcgggCTGGGTGATCTCTCATCAAACAAAATTCAACCCTTTTTTTTCTAACCAAACACGCGGCCACCATTTTATACCCCTGTTTCGGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCaa                                                             <  2:843140/144‑1 (MQ=1)
                    gggtGGTTGATCTCTCATAAAACAAAATGCACCCCATTGTTTTCTAAGCACTCACGCGGCCACCAGTTGATACCCCTGTTTCGGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTg                                                      <  2:4657356/147‑1 (MQ=1)
                          gttaTCTCTCATCAAACAAAATCCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCTTTTGATACCCCAGTTTCCGGGCTTTGCGAAGAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAAc                                                <  2:5067557/147‑1 (MQ=1)
                          gttaTCTCTCAACAAAAAAAAGGCAACCATTTTTTTTCTACGCACTCACGCGGCCCCCAGTTGATACCCCAGTTTCCGGGATTTGCGAAGAAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAAc                                                <  1:4932760/147‑1 (MQ=1)
                                tcccatcaAAAAAATTGCAACCCTTTGTTTTCTACGCACTCCCGGGGCCCCCCGTTGATCCCCCTTTTTCGGGGGTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGc                                          <  1:2983572/147‑1 (MQ=1)
                                    aaaaaccaaaATGCAACCCTTTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCGGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGc                                      <  1:5006394/144‑1 (MQ=0)
                                      cacacaaaATCCACCCCCTTTTTTTCAAAGAACTCACGCGGCCCCCAGTGGATACCCCATTTGCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGctt                                    <  1:4951501/147‑1 (MQ=1)
                                         accaaaGGCAACCCCTTGTTTTCTAAGCACTCCCGGGGCCACCTGTTGAAACCCCTGTTTCCGGGCTTTGCGAAGCACCGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGcttctt                                 <  2:1153003/147‑1 (MQ=1)
                                           aaaaaccAACCCATTTTTTTTCTAACCACACACGCGGCCACCTTTTCAACCCCCTGTTGCGGGGTTTTGCCACCCAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTg                                <  1:1723362/144‑1 (MQ=0)
                                                 caccccTTTGTTTTCTAAGCACTCAAGCGGCCACCTGTTGATACCCCTGTTGCCGGGCTTGGCGAAGCAATGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                         <  2:2646873/147‑1 (MQ=1)
                                                 caccccTTTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                         <  2:1318986/147‑1 (MQ=1)
                                                 caccccATTGTTTTCTACGCACCCACGCGGCCACCCGTTGATCCCCCCGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACa                         <  2:2895060/147‑1 (MQ=1)
                                                      ccttttttttCTACCCACCCCCGCGCCCACCCGTTGATCCCCCTGTTTCGGGTCTTGGCGACCAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                     <  1:850857/144‑1 (MQ=1)
                                                      ccttttttttCTACCAACACACGCGGCCACCTTTTTAAACCCCTTTTGCGGGTCTTTGCGAAGCAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                     <  1:4124101/144‑1 (MQ=1)
                                                      cctttGTTTTCTACGCACCCCCGCGGCCCCCCGTTGCTACCCCAGTGGCGGGGCTTTGCGACGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                     <  2:1671363/147‑1 (MQ=1)
                                                      cctttGTTTTCTAAGCACTCCAGGGGCCACCCGTTGAGACCCCTGGTTCCGGGGGTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                     <  2:5329324/147‑1 (MQ=1)
                                                      cctttGTTTTCAACGCACTCACGCGGCCACCTGTTGATACCCCTGTTTCGGGGCTTTGCTACGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                     <  2:4026177/147‑1 (MQ=1)
                                                      cccttGTTTTCTCCGCACCCACGCGGCCACCCGTGGAGACCCCCGTTTCGGGGCTTTGCGAAGCAACGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGtt                     <  2:3246502/147‑1 (MQ=1)
                                                       cttttttttCTATCCACACAGGCGCCCACCCGTCCTTCCCCCTTTTGCGGGGCTTGGCGACCAAACCCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTg                    <  1:4162289/145‑1 (MQ=0)
                                                       cttttttttCTAAGCACTCACGGGGCCACCAGTTGATACCCCCGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTg                    <  1:3482870/145‑1 (MQ=0)
                                                                   agaaCTCAAGCGGCCACCTTTTTATACCCCAGTTTCGGGGGTTTGCTAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGt        <  2:895083/147‑1 (MQ=1)
                                                                      acacacgcgcCCACCTGTTGAAACCCCTGTTTCGGGGCTTTGCGATGAAAACCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGcc     <  1:4761147/144‑1 (MQ=1)
                                                                         cacgCGGCCCCCTTTTGACACCCATGTTGCGGGGCTTTCCGAACAAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAAc  <  1:5099522/147‑1 (MQ=1)
                                                                                     |                                                                                                                                         
GCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGTTTCGAACTCCTCAAGACCTTGTTCAACATAAGCCTGCTTCTTGGTTAACAGGTTGTACACCAGACGTGCCAAC  >  NC_002947/5222945‑5223166

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