Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,188,733 T→C 46.2% L42L (TTA→TTG nanA ← N‑acetylneuraminate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,188,7330TC46.2% 0.1 / 30.4 39L42L (TTA→TTGnanAN‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base T (9/12);  new base C (9/9);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCT  >  minE/2188673‑2188799
                                                            |                                                                  
tACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGta                                                            <  1:1064843/69‑2 (MQ=255)
tACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTc                                                            <  1:1570276/69‑1 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACGCACGTATAAACCGTCGATGccc                                                     <  1:197291/68‑1 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGccc                                                     <  1:1961253/68‑1 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGccc                                                     <  1:3387605/68‑1 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGccc                                                     <  1:2564508/68‑1 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGccc                                                     <  1:656908/68‑1 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGccc                                                     <  1:1954099/68‑1 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGccc                                                     <  1:1869254/68‑1 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGccc                                                     <  1:1761690/68‑1 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGccc                                                     <  1:1664990/68‑1 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGccc                                                     <  1:1374054/68‑1 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCct                                                    >  1:1932953/1‑69 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCct                                                    >  1:1940893/1‑69 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCct                                                    >  1:2775173/1‑69 (MQ=255)
        cACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCct                                                    >  1:2739998/1‑69 (MQ=255)
           gCTCTGAAAGGCTCTGTACAAAGGCTTCGCCTGTCGAACCCCCCACGTACAAACCGTCGATGCCctgct                                                 <  1:3110317/69‑1 (MQ=255)
           gCTCTGAAAGGCTCTGTACAAAGGCTTCGCCTGTCGAACCCCCCACGTACAAACCGTCGATGCCctgct                                                 <  1:1769255/69‑1 (MQ=255)
           gCTCTGAAAGGCTCTGTACAAAGGCTTCGCCTGTCGAACCCCCCACGTACAAACCGTCGATGCCctgct                                                 <  1:2860584/69‑1 (MQ=255)
           gCTCTGAAAGGCTCTGTACAAAGGCTTCGCCTGTCGAACCCCCCACGTACAAACCGTCGATGCCctgct                                                 <  1:1133991/69‑1 (MQ=255)
           gCTCTGAAAGGCTCTGTACAAAGGCTTCGCCTGTCGAACCCCCCACGTACAAACCGTCGATGCCctgct                                                 <  1:2319035/69‑1 (MQ=255)
           gCTCTGAAAGGCTCTGTACAAAGGCTTCGCCTGTCGAACCCCCCACGTACAAACCGTCGATGCCctgct                                                 <  1:1267334/69‑1 (MQ=255)
           gCTCTGAAAGGCTCTGTACAAAGGCTTCGCCTGTCGAACCCCCCACGTACAAACCGTCGATGCCctgct                                                 <  1:2070608/69‑1 (MQ=255)
           gCTCTGAAAGGCTCTGTACAAAGGCTTCGCCTGTCGAACCCCCCACGTACAAACCGTCGATGCCctgct                                                 <  1:206424/69‑1 (MQ=255)
           gCTCTGAAAGGCTCTGTACAAAGGCTTCGCCTGTCGAACCCCCCACGTACAAACCGTCGATGCCctgct                                                 <  1:1055692/69‑1 (MQ=255)
                    ggCTTTGTACAAAGGCATCGCCGGTCGAACCACCCACGTATAAACCGTCGATg                                                        >  1:2659932/1‑53 (MQ=255)
                             cAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCCTGCTGaa                                              >  1:1127828/1‑54 (MQ=255)
                             cAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCCTGCTGaa                                              >  1:3007947/1‑54 (MQ=255)
                             cAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCCTGCTGaa                                              >  1:3343028/1‑54 (MQ=255)
                             cAAAGGCCTCGCCGGTCGAACCACACACGTATAAACCGTCGATGCCCTGCTGa                                               >  1:2600904/1‑53 (MQ=255)
                                                     ccacGTACAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGGCGCAGa                    >  1:3404817/1‑56 (MQ=255)
                                                     ccacGTACAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGGCGCAGa                    >  1:789998/1‑56 (MQ=255)
                                                     ccacGTACAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGGCGCAGa                    >  1:756946/1‑56 (MQ=255)
                                                     ccacGTACAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGGCGCAGa                    >  1:3083644/1‑56 (MQ=255)
                                                     ccacGTACAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGGCGCAGa                    >  1:2870645/1‑56 (MQ=255)
                                                     ccacGTACAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGGCGCAGa                    >  1:2366887/1‑56 (MQ=255)
                                                     ccacGTACAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGGCGCAGa                    >  1:2132783/1‑56 (MQ=255)
                                                     ccacGTACAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGGCGCAGa                    >  1:1911746/1‑56 (MQ=255)
                                                     ccacGTACAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGGCGCAGa                    >  1:189022/1‑56 (MQ=255)
                                                          tataAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGACGTAGACTCGCTTTATCCAGTGCt  >  1:692610/1‑69 (MQ=255)
                                                          tataAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCt  >  1:527941/1‑69 (MQ=255)
                                                            |                                                                  
TACCTGTTCACGCTCGGAAAGGCTTTGTACAAAGGCCTCGCCGGTCGAACCACCCACGTATAAACCGTCGATGCCCTGCTGAATATTGAACTGAACCAGGCGACGCAGACTCGCTTTATCCAGTGCT  >  minE/2188673‑2188799

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