Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,974,119 C→T 64.0% T140T (ACC→ACT metK → methionine adenosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,974,1190CT64.0% 35.4 / 35.5 61T140T (ACC→ACTmetKmethionine adenosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base T (25/14);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ACCAGGGTCTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGC  >  minE/1974055‑1974185
                                                                |                                                                  
aCCAGGGTCTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACCTATg                                                                <  1:1730067/69‑1 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCCGCACCTATCACTTATGCacac                                                           >  1:1018211/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGcaca                                                            >  1:236908/4‑68 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:1628525/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:2267860/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:675385/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:2067339/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:577321/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:677312/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:518223/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:1797467/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:2388598/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:168460/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:2509684/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:348056/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:812040/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:1477016/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:1346323/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:3431204/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:127051/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:1264759/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:1191973/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:3036823/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:3106337/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:3118770/4‑69 (MQ=255)
     gggcTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACTTATGCacac                                                           >  1:3196129/4‑69 (MQ=255)
      gTCTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACCTATg                                                                <  1:1317703/63‑1 (MQ=255)
             gTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTgg                                                    <  1:1933265/68‑1 (MQ=255)
                                 aaaCCGACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTg                                >  1:2051236/1‑68 (MQ=255)
                                 aaaCCGACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTg                                >  1:2751071/1‑68 (MQ=255)
                                 aaaCCGACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTg                                >  1:1175665/1‑68 (MQ=255)
                                 aaaCCGACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTg                                >  1:698513/1‑68 (MQ=255)
                                      gACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGTTACAGCGTCAGGCTGAAGTGCGTaaa                          <  1:1892288/69‑1 (MQ=255)
                                      gACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTaaa                          <  1:2201153/69‑1 (MQ=255)
                                      gACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTaaa                          <  1:3381188/69‑1 (MQ=255)
                                      gACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTaaa                          <  1:1229350/69‑1 (MQ=255)
                                      gACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTaaa                          >  1:3382280/1‑69 (MQ=255)
                                      gACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTaaa                          <  1:897671/69‑1 (MQ=255)
                                      gACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTaaa                          <  1:1640693/69‑1 (MQ=255)
                                      gACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTaaa                          <  1:1847896/69‑1 (MQ=255)
                                      gACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTaaa                          >  1:2130204/1‑69 (MQ=255)
                                            cTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGt                                            >  1:2653797/1‑45 (MQ=255)
                                            cTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGc                    >  1:1595702/1‑69 (MQ=255)
                                            cTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGc                    >  1:2875436/1‑69 (MQ=255)
                                            cTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGc                    >  1:2134290/1‑69 (MQ=255)
                                            cTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGc                    >  1:111870/1‑69 (MQ=255)
                                            cTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGc                    >  1:105660/1‑69 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:407663/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:710704/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:2341074/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:3304554/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:90732/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:1037600/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:287378/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:2660799/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:2592136/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:2555529/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:2293269/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:3226219/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:2092897/68‑1 (MQ=255)
                                                    aGCACCTATCACTTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGc             <  1:1986149/68‑1 (MQ=255)
                                                            tcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTgc    >  1:1754579/1‑69 (MQ=255)
                                                              acctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTtcg   >  1:2301426/1‑68 (MQ=255)
                                                              acctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTgcgc  >  1:622164/1‑69 (MQ=255)
                                                              acctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTgcgc  >  1:1307841/1‑69 (MQ=255)
                                                              acctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTgcgc  >  1:1576733/1‑69 (MQ=255)
                                                              acctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTgcgc  >  1:780191/1‑69 (MQ=255)
                                                              acctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTgcgc  >  1:2512282/1‑69 (MQ=255)
                                                              acctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTgcgc  >  1:3030382/1‑69 (MQ=255)
                                                              acctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTgcg   >  1:1905706/1‑68 (MQ=255)
                                                                |                                                                  
ACCAGGGTCTGATGTTTGGCTACGCAACTAATGAAACCGACGTGCTGATGCCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGC  >  minE/1974055‑1974185

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