Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,384,474 G→A 13.3% L363L (CTG→TTG)  metE ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,384,4740GA13.3% 251.5 / 25.3 90L363L (CTG→TTG) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (40/38);  new base A (5/7);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTG  >  minE/2384406‑2384542
                                                                    |                                                                    
ggTAGAGTGACGACGTGCCTGAATAGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAgcg                                                                    <  1:924524/71‑1 (MQ=255)
  tAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCa                                                                  >  1:992301/1‑71 (MQ=255)
    gagTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAg                                                                      <  1:3463208/65‑1 (MQ=255)
    gagTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAg                                                                      <  1:3493102/65‑1 (MQ=255)
    gagTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAg                                                                      <  1:36010/65‑1 (MQ=255)
    gagTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAg                                                                      <  1:1672229/65‑1 (MQ=255)
    gagTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAg                                                                      <  1:1532096/65‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                         <  1:1030434/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                         <  1:1573472/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                         <  1:1553814/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                         <  1:144022/71‑1 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcag                                                          >  1:580502/1‑70 (MQ=255)
           gacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGAgcag                                                          >  1:744520/1‑70 (MQ=255)
            acgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                         <  1:3438990/70‑1 (MQ=255)
                 ccTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGc                                                    >  1:231762/1‑70 (MQ=255)
                   tGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAg                                                 >  1:2679200/1‑71 (MQ=255)
                   tGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAg                                                 >  1:211605/1‑71 (MQ=255)
                     aaTCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGt                                                >  1:1793144/1‑70 (MQ=255)
                         ggggCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCAt                                                                 >  1:521600/1‑49 (MQ=255)
                         ggggCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCAt                                                                 >  1:2143992/1‑49 (MQ=255)
                               cTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGGCAGTTCATGGCAttt                                     >  1:3048397/1‑71 (MQ=255)
                               cTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttt                                     >  1:3414020/1‑71 (MQ=255)
                               cTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttt                                     >  1:3963145/1‑71 (MQ=255)
                                 ccACTCTGACAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAtttt                                    >  1:2709073/1‑70 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:1189413/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:2357012/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:3219195/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:2450233/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:3489890/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:3406489/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:3199369/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:307096/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:691546/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  >  1:1663265/1‑71 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:2969347/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:4081980/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  <  1:1266250/71‑1 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  >  1:1458930/1‑71 (MQ=255)
                                  cACTCTGCCAGAGCTGCCGTGTAACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  >  1:1128089/1‑71 (MQ=255)
                                      ctGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAgg                              >  1:1599981/1‑71 (MQ=255)
                                       tGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttttt                                   >  1:549880/1‑65 (MQ=255)
                                       tGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttttt                                   >  1:282847/1‑65 (MQ=255)
                                            gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                        >  1:3786761/1‑71 (MQ=255)
                                            gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                        >  1:2436248/1‑71 (MQ=255)
                                            gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                        >  1:3858875/1‑71 (MQ=255)
                                            gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                        >  1:1880534/1‑71 (MQ=255)
                                            gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                        >  1:472278/1‑71 (MQ=255)
                                             agCTGCCGTGTCACCACTGTTCAGCGCATCGCgcagc                                                         <  1:4071566/37‑1 (MQ=255)
                                              gCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                  >  1:2982017/1‑59 (MQ=255)
                                                 gCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCAt                                            >  1:56710/1‑46 (MQ=255)
                                                 gCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaaa                   >  1:1520747/1‑71 (MQ=255)
                                                   cGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAcc                 <  1:291584/71‑1 (MQ=255)
                                                    gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaaa                   >  1:3581046/1‑68 (MQ=255)
                                                    gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                  >  1:1334222/1‑69 (MQ=255)
                                                    gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                  >  1:2318466/1‑69 (MQ=255)
                                                    gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                  >  1:3576565/1‑69 (MQ=255)
                                                    gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                  >  1:259854/1‑69 (MQ=255)
                                                    gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                <  1:3006751/68‑1 (MQ=255)
                                                    gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                <  1:679723/68‑1 (MQ=255)
                                                    gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                <  1:2054522/68‑1 (MQ=255)
                                                    gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                <  1:1143566/68‑1 (MQ=255)
                                                    gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                <  1:354384/68‑1 (MQ=255)
                                                    gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                <  1:2547815/68‑1 (MQ=255)
                                                    gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                <  1:2837048/68‑1 (MQ=255)
                                                     tgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttt                                     <  1:3566495/49‑1 (MQ=255)
                                                     tgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttt                                     <  1:3202546/49‑1 (MQ=255)
                                                       tCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCt             <  1:247978/71‑1 (MQ=255)
                                                       tCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCt             <  1:3268882/71‑1 (MQ=255)
                                                       tCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCt             <  1:3208637/71‑1 (MQ=255)
                                                         accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt           >  1:2531894/1‑71 (MQ=255)
                                                         accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt           >  1:116741/1‑71 (MQ=255)
                                                         accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt           <  1:2972727/71‑1 (MQ=255)
                                                         accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt           >  1:1407526/1‑71 (MQ=255)
                                                         accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt           >  1:503024/1‑71 (MQ=255)
                                                         accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtt            >  1:2934076/1‑70 (MQ=255)
                                                         accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtt            <  1:777133/70‑1 (MQ=255)
                                                         accacTGTTCAGCGCATCGCGCAGCAGTCCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt           <  1:513235/71‑1 (MQ=255)
                                                          ccacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt          <  1:3164271/71‑1 (MQ=255)
                                                          ccacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt          <  1:2835817/71‑1 (MQ=255)
                                                          ccacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt          <  1:1608458/71‑1 (MQ=255)
                                                          ccacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt          <  1:917543/71‑1 (MQ=255)
                                                           cacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTGATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt           <  1:1661228/69‑1 (MQ=255)
                                                             cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGGaa                    <  1:1769867/58‑1 (MQ=255)
                                                             cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaa                    <  1:3650997/58‑1 (MQ=255)
                                                             cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaa                    <  1:552048/58‑1 (MQ=255)
                                                             cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaa                    <  1:570389/58‑1 (MQ=255)
                                                               gTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCa                >  1:3826583/1‑60 (MQ=255)
                                                               gTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCAttt     >  1:2579579/1‑71 (MQ=255)
                                                               gTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACtt     >  1:3077669/1‑71 (MQ=255)
                                                                ttCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGc                      >  1:3507569/1‑53 (MQ=255)
                                                                ttCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTc    >  1:3027212/1‑71 (MQ=255)
                                                                  cAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg  >  1:1541899/1‑71 (MQ=255)
                                                                  cAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg  >  1:2418600/1‑71 (MQ=255)
                                                                  cAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg  >  1:1091904/1‑71 (MQ=255)
                                                                   aGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg  <  1:1219284/70‑1 (MQ=255)
                                                                   aGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg  <  1:3886369/70‑1 (MQ=255)
                                                                    |                                                                    
GGTAGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTG  >  minE/2384406‑2384542

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: