Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,643,860 A→G 100% R9R (CGT→CGC purC ← phosphoribosylaminoimidazole‑succinocarboxamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,643,8600AG91.7% 30.7 / ‑3.2 12R9R (CGT→CGCpurCphosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (11/0);  total (11/1)
Fisher's exact test for biased strand distribution p-value = 8.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTT  >  minE/1643804‑1643916
                                                        |                                                        
aaTTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGAt                                                      >  1:119149/1‑61 (MQ=255)
aaTTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGAt                                                      >  1:14008/1‑61 (MQ=255)
aaTTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGAt                                                      >  1:141575/1‑61 (MQ=255)
aaTTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGAt                                                      >  1:15754/1‑61 (MQ=255)
aaTTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGAt                                                      >  1:243927/1‑61 (MQ=255)
aaTTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGAt                                                      >  1:244679/1‑61 (MQ=255)
aaTTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGAt                                                      >  1:424577/1‑61 (MQ=255)
        cACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGATACAACTCAg                                             >  1:327943/1‑62 (MQ=255)
        cACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCGCGATACAACTCAg                                             >  1:348631/1‑62 (MQ=255)
                  tCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTg                                   <  1:40890/62‑1 (MQ=255)
                                                    tACCGCGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGtt  >  1:378253/1‑61 (MQ=255)
                                                    tACCGCGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGtt  >  1:389372/1‑61 (MQ=255)
                                                        |                                                        
AATTCGAGCACCAACAGGTCCGGGTTTTCCGTGCTGTATACGGTTTTCGCTTTACCACGATACAACTCAGCTTGCTTTTGCATCTTTATCACTCCTGGGTGTGAATTAACGTT  >  minE/1643804‑1643916

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