Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP025268 390,887 G→A A356V (GCG→GTG)  CXP41_RS02035 ← IS3‑like element IS3 family transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP025268390,8870GA81.2% 12.3 / ‑3.9 16A356V (GCG→GTG) CXP41_RS02035IS3‑like element IS3 family transposase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (12/1);  total (13/3)
Fisher's exact test for biased strand distribution p-value = 7.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACGCCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACG  >  NZ_CP025268/390751‑391014
                                                                                                                                        |                                                                                                                               
ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGTCTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCg                                                                                                                              >  1:1129614/1‑140 (MQ=255)
ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCg                                                                                                                              >  1:1888743/1‑140 (MQ=9)
ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCg                                                                                                                              >  1:1953485/1‑140 (MQ=18)
ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCg                                                                                                                              >  1:3538698/1‑140 (MQ=12)
ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCg                                                                                                                              >  1:4181596/1‑140 (MQ=12)
ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCg                                                                                                                              >  1:5341277/1‑140 (MQ=12)
ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCg                                                                                                                              >  2:1111497/1‑140 (MQ=9)
ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCg                                                                                                                              >  2:1112598/1‑140 (MQ=12)
ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGCTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCg                                                                                                                              >  2:1109759/1‑140 (MQ=255)
  ggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGtt                                                                                                                            <  1:4263188/140‑1 (MQ=2)
  ggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGtt                                                                                                                            >  2:1324427/1‑140 (MQ=0)
  ggCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTGCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGtt                                                                                                                            >  2:11203/1‑140 (MQ=9)
  ggCTGATTGATCATACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGtt                                                                                                                            >  2:5515345/1‑140 (MQ=0)
                                     gcctAAGCGAGGTTATTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACCCTAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACGCCGTTGCCCGCATTATTTCCCGGCTGCTAGAGTGTtcgtc                                                                                         >  2:4080567/4‑137 (MQ=11)
                                                                                agtcatacCAACTGTGCCGCCCCAACCGATTGTAATCACATTCGATATAATTAAACGCCGTTGCCCGCATTATTTCCCGGCTGATAAAGTTTTCTCCATGGATTCATTCCACTTTCAGCATATGAATGAAGCTTTCCACg                                              <  2:644016/134‑1 (MQ=11)
                                                                                                                            agataATTAAAGGCCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACAGTCCACTTTCAGCGAATGAAAGAAGCGTTCCACGCAGGCAGTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACg  <  2:5345041/138‑1 (MQ=11)
                                                                                                                                        |                                                                                                                               
GGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACGCCGTTGCCCGCATTATTTCCCGGCTGATAAAGTGTTCTCCATGGATACATTCCACTTTCAGCGAATGAAAGAAGCTTTCCACGCAGGCATTATCGTAGCAGCAACCTTTTGCGCTCATACTTCCACG  >  NZ_CP025268/390751‑391014

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