Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273276,8000AC52.9% 14.9 / 19.7 17Y299D (TAT→GAT) yagAIS481‑like element ISEc19 family transposase
Reads supporting (aligned to +/- strand):  ref base A (6/2);  new base C (0/9);  total (6/11)
Fisher's exact test for biased strand distribution p-value = 2.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGG  >  NZ_CP009273/276717‑276877
                                                                                   |                                                                             
gATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGa                                                                         >  1:186028/1‑90 (MQ=255)
gATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGa                                                                         >  2:184440/1‑90 (MQ=255)
gATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGa                                                                         >  2:234794/1‑90 (MQ=255)
    tCCCCTTTCCTGACCACCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGccc                                                                     <  2:209180/90‑1 (MQ=255)
         tttCCTGACCCCCCCCCCTTCATCGTATCCGGGGGGCGTTGTGTTGCCGTTGTCCTGCCGCCCGGACGGCTGATCCCGCGAGCCGGGtcc                                                                <  1:210918/90‑3 (MQ=255)
         ttcccTGCCCACCCCCCCTTCATCGTTTTCCGGGGGCGTTGTTTTGCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGAGCCCGGTAc                                                                <  1:101699/87‑1 (MQ=255)
              ttaccaccCCCCCTCCATCGTTTTCCGGGTGCGTTGTTTTCCCCTTCTCCTGCCGCGCGGAAGGCTGATACCGCGAGCCCGGTACCGCCa                                                           <  1:26372/88‑1 (MQ=255)
                                       cGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTc                                  <  2:39669/90‑1 (MQ=255)
                                        gggggCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                 >  2:203079/1‑90 (MQ=255)
                                         ggggTGTTTTTTTGCCGGGGTCCGCCCGCTCGGACGGCTGATCCCGCGACCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa                                <  2:207853/90‑1 (MQ=255)
                                         ggggCTTTGTTTCCCCGCTGTCCCCCCGCGCGGAGGGCAGATCCCGCGAGCCGGGCCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa                                <  1:184440/90‑1 (MQ=255)
                                         ggggCTTTGTGTTCCCGCTGTACGGCCGCGCGGACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa                                <  2:106498/90‑1 (MQ=255)
                                         ggggCGTTGTGTTCCCGCTGTACTGCCGCGCTGCCGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa                                <  2:74255/90‑1 (MQ=255)
                                             cGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGtt                            >  2:8060/1‑90 (MQ=255)
                                                           tGTCCTGCCGCGCTGACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGcc              <  2:213798/90‑1 (MQ=255)
                                                                      gcTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAg   >  2:39292/1‑90 (MQ=255)
                                                                       cTGCCGGCTGATCCCGCGGGCCCGGTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAgg  <  2:309414/90‑1 (MQ=255)
                                                                                   |                                                                             
GATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGG  >  NZ_CP009273/276717‑276877

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