Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,185,010 C→T 10.6% Y325Y (TAC→TAT yeaG → conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,185,0100CT10.6% 347.7 / 24.0 122Y325Y (TAC→TATyeaGconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base C (41/68);  new base T (11/2);  total (52/70)
Fisher's exact test for biased strand distribution p-value = 1.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 126 total aligned reads displayed.

GAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCA  >  minE/1184944‑1185078
                                                                  |                                                                    
gAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCg                                                                  >  1:1649587/1‑71 (MQ=255)
   cGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:2610481/71‑1 (MQ=255)
   cGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:218075/71‑1 (MQ=255)
   cGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:899725/71‑1 (MQ=255)
   cGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:3669232/71‑1 (MQ=255)
   cGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:5153118/71‑1 (MQ=255)
    gAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGaa                                                              <  1:4863717/71‑1 (MQ=255)
              gggTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGaa                                                              <  1:2938971/61‑1 (MQ=255)
               tgTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAtt                                                   <  1:3752969/70‑1 (MQ=255)
               ggTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAtt                                                   <  1:1239636/71‑1 (MQ=255)
                    cTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAttgcttg                                              <  1:4523425/71‑1 (MQ=255)
                    cTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAttgcttg                                              <  1:3220125/71‑1 (MQ=255)
                    cTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTAttgcttg                                              <  1:3882260/71‑1 (MQ=255)
                     tttCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:5108633/53‑1 (MQ=255)
                     tttCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:2367878/53‑1 (MQ=255)
                        ccGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCa                                          >  1:280797/1‑71 (MQ=255)
                          gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:2545427/48‑1 (MQ=255)
                          gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:4829017/48‑1 (MQ=255)
                          gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:919603/48‑1 (MQ=255)
                          gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:2128243/48‑1 (MQ=255)
                          gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:4089734/48‑1 (MQ=255)
                          gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:4964384/48‑1 (MQ=255)
                          gTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:3250381/48‑1 (MQ=255)
                           taataaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:4254376/47‑1 (MQ=255)
                           taataaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGa                                                               <  1:4002593/47‑1 (MQ=255)
                              taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    <  1:1346692/71‑1 (MQ=255)
                              taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    <  1:2230309/71‑1 (MQ=255)
                              taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    <  1:2328440/71‑1 (MQ=255)
                              taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    <  1:4328090/71‑1 (MQ=255)
                              taaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    <  1:1500764/71‑1 (MQ=255)
                               aaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTgcgc                                           <  1:3312983/63‑1 (MQ=255)
                               aaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTgcgc                                           <  1:4457892/63‑1 (MQ=255)
                               aaCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTgcgc                                           <  1:3078019/63‑1 (MQ=255)
                                aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    >  1:106081/1‑69 (MQ=255)
                                aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    >  1:2724458/1‑69 (MQ=255)
                                aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    >  1:1273432/1‑69 (MQ=255)
                                aCAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    >  1:1417489/1‑69 (MQ=255)
                                 cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGt                                                      >  1:3336316/1‑50 (MQ=255)
                                 cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTcag                                    >  1:3904249/1‑66 (MQ=255)
                                 cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTc                                      >  1:4035432/1‑66 (MQ=255)
                                 cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    >  1:1907883/1‑68 (MQ=255)
                                 cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    >  1:3690936/1‑68 (MQ=255)
                                 cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    >  1:4716675/1‑68 (MQ=255)
                                 cAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    >  1:2544656/1‑68 (MQ=255)
                                 cAAAAACAAAGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCg                                    >  1:4307666/1‑68 (MQ=255)
                                  aaaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:4189999/1‑69 (MQ=255)
                                  aaaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:2582564/1‑69 (MQ=255)
                                   aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCa                                          >  1:3467636/1‑60 (MQ=255)
                                   aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:3188188/1‑68 (MQ=255)
                                   aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:1947465/1‑68 (MQ=255)
                                   aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:2166856/1‑68 (MQ=255)
                                   aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:4838604/1‑68 (MQ=255)
                                   aaaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:3909642/1‑68 (MQ=255)
                                    aaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTg                                                                <  1:2529728/37‑1 (MQ=255)
                                    aaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTg                                                                <  1:3783730/37‑1 (MQ=255)
                                    aaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTg                                                                <  1:1701134/37‑1 (MQ=255)
                                    aaaCAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:1081638/1‑67 (MQ=255)
                                       caacGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATc                            >  1:2834938/1‑70 (MQ=255)
                                         acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTAt                                                    >  1:3568936/1‑44 (MQ=16)
                                         acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTgcg                                            >  1:3466443/1‑52 (MQ=19)
                                         acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                            >  1:408721/1‑68 (MQ=255)
                                         acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                            >  1:658818/1‑68 (MQ=255)
                                         acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                            >  1:531556/1‑68 (MQ=255)
                                         acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                            >  1:1061672/1‑68 (MQ=255)
                                         acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                            >  1:4423473/1‑68 (MQ=255)
                                         acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                            >  1:3615779/1‑68 (MQ=255)
                                         acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                            >  1:1582473/1‑68 (MQ=255)
                                         acGAGGCGTTCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGATGAAATc                            >  1:3124410/1‑68 (MQ=255)
                                           gAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAAt                       <  1:2536838/71‑1 (MQ=255)
                                             aGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCt                     <  1:4073208/71‑1 (MQ=255)
                                                 ttCCTCGACCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGAtga                                <  1:4543288/56‑1 (MQ=255)
                                                   ccTCGATCGTATCTATATCGTGAAGGTGCCGTATTGCCTGCGCGTGTCCGAtga                                <  1:1753935/54‑1 (MQ=255)
                                                    cTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga               <  1:4069448/70‑1 (MQ=255)
                                                    cTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga               <  1:2976790/70‑1 (MQ=255)
                                                    cTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGaa              <  1:736434/71‑1 (MQ=255)
                                                     tCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACgaga               <  1:3825626/69‑1 (MQ=255)
                                                     tCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGaa              <  1:2424165/70‑1 (MQ=255)
                                                      cGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGaa              <  1:3375686/69‑1 (MQ=255)
                                                              tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                         <  1:1879344/50‑1 (MQ=255)
                                                              tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                         <  1:1065409/50‑1 (MQ=255)
                                                              tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                         <  1:3965795/50‑1 (MQ=255)
                                                              tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                         <  1:709103/50‑1 (MQ=255)
                                                              tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                         <  1:337610/50‑1 (MQ=255)
                                                              tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                         <  1:1884815/50‑1 (MQ=255)
                                                              tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAAtt           <  1:2767203/64‑1 (MQ=255)
                                                              tttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAAt    >  1:3043214/1‑71 (MQ=255)
                                                               ttACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCaaa                         <  1:4195978/49‑1 (MQ=255)
                                                                tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGc         >  1:2727891/1‑64 (MQ=255)
                                                                tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGc         >  1:321991/1‑64 (MQ=255)
                                                                tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGc         >  1:1378579/1‑64 (MQ=255)
                                                                tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGc         >  1:4222585/1‑64 (MQ=255)
                                                                tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGc         >  1:127299/1‑64 (MQ=255)
                                                                tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATca  >  1:728910/1‑71 (MQ=255)
                                                                tACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATca  >  1:874782/1‑71 (MQ=255)
                                                                 aCATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:4687909/1‑38 (MQ=255)
                                                                 aCATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:2762159/1‑38 (MQ=255)
                                                                 aCATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:1406821/1‑38 (MQ=255)
                                                                 aCATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGaa                                  >  1:4247604/1‑38 (MQ=255)
                                                                  cATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAACCTACGAGAAATTGc         >  1:5107649/1‑62 (MQ=255)
                                                                  cATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAACCTACGAGAAATTGc         >  1:3537176/1‑62 (MQ=255)
                                                                  |                                                                    
GAACGAGTCCGAATGGGTCACTTTCCGTAATAACAAAAACAACGAAGCCTTCCTCGATCGTGTTTACATCGTGAAGGTGCCGTATTGCTTGCGCATTTCCGAAGAGATCAAAATCTACGAGAAATTGCTTAATCA  >  minE/1184944‑1185078

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