Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 389,430 T→C 52.5% C316C (TGT→TGC ybdL → methionine aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE389,4300TC52.5% 3.7 / 50.1 61C316C (TGT→TGCybdLmethionine aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (17/12);  new base C (17/15);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTACAGCGCGGTTTCTACCCTGGATGATG  >  minE/389375‑389485
                                                       |                                                       
aGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACt                                            >  1:628187/1‑69 (MQ=255)
aGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACt                                            >  1:1441255/1‑69 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGTCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtg                                    >  1:1553201/1‑69 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTAc                                                 >  1:1079127/1‑56 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtg                                    >  1:2571483/1‑69 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtg                                    >  1:59546/1‑69 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtg                                    >  1:3091629/1‑69 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtg                                    >  1:2731388/1‑69 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtg                                    >  1:2005658/1‑69 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtg                                    >  1:1545323/1‑69 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtg                                    >  1:1344522/1‑69 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtg                                    >  1:1249803/1‑69 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtg                                    >  1:1201324/1‑69 (MQ=255)
        atTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCt                                     >  1:711542/1‑68 (MQ=255)
           cTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAgg                                                    <  1:741025/50‑1 (MQ=255)
           cTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAgg                                                    <  1:3412186/50‑1 (MQ=255)
           cTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAgg                                                    <  1:3281918/50‑1 (MQ=255)
            tGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtggtg                                 <  1:2498958/68‑1 (MQ=255)
            tGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtggtg                                 <  1:2369457/68‑1 (MQ=255)
            tGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtggtg                                 <  1:2344159/68‑1 (MQ=255)
            tGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtggtg                                 <  1:651772/68‑1 (MQ=255)
            tGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCtggtg                                 <  1:933917/68‑1 (MQ=255)
                 aaTGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTAc                           <  1:3043140/69‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGTTACATACTTTTTGCtggtgg                                <  1:2715565/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:1417506/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:2580694/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:1041071/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:2810705/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:3075400/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:3259203/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:2511119/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:1354836/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:123680/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:1592894/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:2186516/57‑1 (MQ=255)
                        tAAATGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:1959119/57‑1 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGGAGGTACATACTTTTTGCTGGTGGAt                              >  1:2565983/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTa                                                  >  1:1213136/1‑38 (MQ=38)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTAGATACTTTTTGCTGGTGGAt                              >  1:1848984/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTg                                       >  1:925150/1‑49 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGa                               >  1:2786762/1‑57 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:968655/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:877345/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:804382/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:1095171/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:551081/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:2053215/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:2080796/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:3170902/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:2193268/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:2500785/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:269874/1‑58 (MQ=255)
                         aaaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCTGGTGGAt                              >  1:2681791/1‑58 (MQ=255)
                          aaTGAAAGTCGTCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtgttgg                                <  1:2909357/55‑1 (MQ=255)
                          aaTGAAAGTCGGCTGGAGATTTTACCGTGCGAAGGTACATACTTTTTGCtggtgg                                <  1:794322/55‑1 (MQ=255)
                          aaTGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTACAGCGCGGt                   <  1:552667/68‑1 (MQ=255)
                          aaTGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTACAGCGCGGt                   <  1:1204237/68‑1 (MQ=255)
                          aaTGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTACAGCACGGt                   <  1:2692978/68‑1 (MQ=255)
                                           gaTTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTACAGCGCGGTTTCTACCCTGGatgatg  >  1:1225151/1‑68 (MQ=255)
                                           gaTTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTACAGCGCGGTTTCTACCCTGGatgatg  >  1:572787/1‑68 (MQ=255)
                                           gaTTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTACAGCGCGGTTTCTACCCTGGatgatg  >  1:2468330/1‑68 (MQ=255)
                                                       |                                                       
AGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTACAGCGCGGTTTCTACCCTGGATGATG  >  minE/389375‑389485

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: