Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,038,1280CA46.7% ‑5.1 / 15.8 15R303S (AGG→AGT) ‡PP_5290transposase
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (3/4);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GAGGTACGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATGGGCGTAATACTGCTTCCAGGTCGCGGTTCGACTGCCGCTCATCGAGGCGTTGTAC  >  NC_002947/6037985‑6038269
                                                                                                                                               |                                                                                                                                             
gatgTACGCGTCAGCGTGCCAGAGAATTGCCCTCCTCAAGAGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGCCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCtcatc                                                                                                                                         <  1:197302/147‑1 (MQ=1)
    cacgcggcAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGGATCAGCTGGCCGGCCTGGCCCCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATGAGGc                                                                                                                                     <  2:1132000/143‑1 (MQ=1)
     aCGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                    <  1:1263459/150‑1 (MQ=0)
     aCGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                    <  2:2735105/150‑1 (MQ=0)
                                cctaGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGc                                                                                                         >  1:3827732/4‑150 (MQ=1)
                                 ccAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                    >  1:1997527/1‑122 (MQ=0)
                                 ccAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                    <  2:1997527/122‑1 (MQ=0)
                                                     tcggTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGCaa                                                                                                       >  1:3829176/4‑131 (MQ=1)
                                                                    tGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                    <  1:1801496/87‑1 (MQ=1)
                                                                    tGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                    >  2:1801496/1‑87 (MQ=1)
                                                                            aggcccAGCCTCAGCTGGCCGGCTTGGACCCTCCGGTATATATTTCGGGCTTCCAGTCGGCCTGATGCCTCATCACGCAGAAGGCCAGCCGAGCCAAGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGt                                                             <  2:2536718/147‑1 (MQ=1)
                                                                                    ccgcAGCGGGCCGGCTGGGCCCCTCCGGTATCTATCTCGGGCTTCCAGCCGGCCTGATGCCTCATCAGGCCGAAGGCCATCCGAGCCAAGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATGGGCg                                                     <  2:231184/147‑1 (MQ=1)
                                                                                     cTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGAGGCATCAGGcc                                                                                                                                    >  2:1460863/1‑70 (MQ=1)
                                                                                                     ggtaCCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGCAAGGATc                                                                                                  >  2:302427/4‑88 (MQ=1)
                                                                                                                                       gcatGATGCCTCATCAGGCCGAAGGCCAGCCGAGACAAGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATGGGCGTAATACTGCTTCCAGTTCGCGGTTCGACTCCCGCTGATCGAGGCGTTGtca  >  2:3065229/4‑148 (MQ=1)
                                                                                                                                               |                                                                                                                                             
GAGGTACGAGTCAGCGTGCCAGAGAATTTCTCCCCAGCATTGGATACTGCCCGTCAGTTTGGGCGGGCTGCGTTGCAGCCCCAGCCTCAGCTGGCCGGCTTGGCACCTCCGGTATATATCTCGGGCTTCCAGTCGGCCTGATGCCTCATCAGGCCGAAGGCCAGCCGAGCCAAGCGACGTGCAAGGATCACCAAGGCTTGAGTGGTTTTCTTCCCCCGGGCTTGGTGATGGGCGTAATACTGCTTCCAGGTCGCGGTTCGACTGCCGCTCATCGAGGCGTTGTAC  >  NC_002947/6037985‑6038269

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