Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,036,585 A→G G171G (GGA→GGG appA → bifunctional acid phosphatase/4‑phytase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,036,5850AG92.3% 37.2 / ‑5.8 13G171G (GGA→GGGappAbifunctional acid phosphatase/4‑phytase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (8/4);  minor base T (1/0);  total (9/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCAACTGGATAACGCGAACGTGACTGACGCGATCCTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAACAGGACGAAAGCTGTTCATTAACGCAGGCA  >  NZ_CP009273/1036462‑1036717
                                                                                                                           |                                                                                                                                    
tgCAATAACCGCACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCAACTGGATAACGCGAACGTGACTGACGCGATCTTCAGCTGGTCATGTGGGTCAATTGCTGAc                                                                                                                       >  2:29609/1‑139 (MQ=255)
                  ccAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCAACTGGGTAACGCGAACGTGACTGACGCGATCCTCAGCAGGGCAGGGGGGTCAATTGCTGACTTTACCGGGCATCGGCa                                                                                                      >  2:280004/1‑138 (MQ=255)
                  ccAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCAACTGGATAACGCGAACGTGACTGACGCGATCCTCAGCAGGGCAGGGGGGTCAATTGCTGACTTTACCGGGCATCGGCa                                                                                                      >  1:476138/1‑138 (MQ=255)
                     ggCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCAACTGGATAACGCGAACGTGACTGACGCGATCCTCAGCAGGGCAGGGGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACg                                                                                                  >  1:56132/1‑139 (MQ=255)
                          aTACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCAACTGGATAACGCGAACGTGACTGACGCGATCCTCAGCAGGGCAGGGGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGtt                                                                                             >  2:162524/1‑139 (MQ=255)
                                                    aaTCCTCTAAAAACTGGCGTTTGCCAACTGGATAACGCGAACGTGACTGACGCGATCCTCAGCAGGGCAGGGGGGTCAATTg                                                                                                                            >  1:17370/1‑82 (MQ=255)
                                                    aaTCCTCTAAAAACTGGCGTTTGCCAACTGGATAACGCGAACGTGACTGACGCGATCCTCAGCAGGGCAGGGGGGTCAATTg                                                                                                                            <  2:17370/82‑1 (MQ=255)
                                                      tCCTCTAAAAACTGGCGTTTGCCAACTGGATAACGCGAACGTGACTGACGCGATCCTCAGTAGGGCAGGGGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAAtttt                                                                 >  2:22034/1‑139 (MQ=255)
                                                                             aaCTGGATAACGCGAACGTGACTGACGCGATCCTCAGCAGGGCAGGGGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTaa                                          <  1:162524/139‑1 (MQ=255)
                                                                                                          aTCCTCAGCAGGGCAGGGGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAACAGGACGAAAg                     <  1:214995/131‑1 (MQ=255)
                                                                                                          aTCCTCAGCAGGGCAGGGGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAACAGGACGAAAg                     >  2:214995/1‑131 (MQ=255)
                                                                                                           tCCTCAGTAGGGCAGGGGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAACAGGACGAAAGCTGTTCAtt            <  1:22034/139‑1 (MQ=255)
                                                                                                                     ggcaggGGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAACAGGACGAAAGCTGTTCATTAACGCAGGCa  >  2:188550/1‑139 (MQ=255)
                                                                                                                           |                                                                                                                                    
TGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCAACTGGATAACGCGAACGTGACTGACGCGATCCTCAGCAGGGCAGGAGGGTCAATTGCTGACTTTACCGGGCATCGGCAAACGGCGTTTCGCGAACTGGAACGGGTGCTTAATTTTCCGCAATCAAACTTGTGCCTTAAACGTGAGAAACAGGACGAAAGCTGTTCATTAACGCAGGCA  >  NZ_CP009273/1036462‑1036717

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