Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,384,489 G→A 13.3% L358L (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,4890GA13.3% 252.4 / 25.1 90L358L (CTG→TTG) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (39/39);  new base A (5/7);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTC  >  minE/2384423‑2384558
                                                                  |                                                                     
ccTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGc                                                                    >  1:231762/1‑70 (MQ=255)
  tGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAg                                                                 >  1:211605/1‑71 (MQ=255)
  tGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAg                                                                 >  1:2679200/1‑71 (MQ=255)
    aaTCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGt                                                                >  1:1793144/1‑70 (MQ=255)
              cTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGGCAGTTCATGGCAttt                                                     >  1:3048397/1‑71 (MQ=255)
              cTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttt                                                     >  1:3963145/1‑71 (MQ=255)
              cTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttt                                                     >  1:3414020/1‑71 (MQ=255)
                ccACTCTGACAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAtttt                                                    >  1:2709073/1‑70 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:3199369/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:2357012/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:2450233/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  >  1:1663265/1‑71 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:2969347/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:307096/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  >  1:1458930/1‑71 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:3219195/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:3406489/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:1266250/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:3489890/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:1189413/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:4081980/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  <  1:691546/71‑1 (MQ=255)
                 cACTCTGCCAGAGCTGCCGTGTAACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                  >  1:1128089/1‑71 (MQ=255)
                     ctGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAgg                                              >  1:1599981/1‑71 (MQ=255)
                      tGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttttt                                                   >  1:282847/1‑65 (MQ=255)
                      tGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttttt                                                   >  1:549880/1‑65 (MQ=255)
                           gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                        >  1:3786761/1‑71 (MQ=255)
                           gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                        >  1:3858875/1‑71 (MQ=255)
                           gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                        >  1:2436248/1‑71 (MQ=255)
                           gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                        >  1:472278/1‑71 (MQ=255)
                           gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                        >  1:1880534/1‑71 (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:2318466/1‑69 (MQ=255)
                                   gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                                  >  1:3576565/1‑69 (MQ=255)
                                   gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                                  >  1:259854/1‑69 (MQ=255)
                                   gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                                  >  1:1334222/1‑69 (MQ=255)
                                   gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                <  1:3006751/68‑1 (MQ=255)
                                   gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                <  1:2837048/68‑1 (MQ=255)
                                   gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                <  1:2054522/68‑1 (MQ=255)
                                   gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                <  1:354384/68‑1 (MQ=255)
                                   gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                <  1:2547815/68‑1 (MQ=255)
                                   gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                <  1:1143566/68‑1 (MQ=255)
                                   gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                <  1:679723/68‑1 (MQ=255)
                                    tgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttt                                                     <  1:3202546/49‑1 (MQ=255)
                                    tgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttt                                                     <  1:3566495/49‑1 (MQ=255)
                                      tCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCt                             <  1:3268882/71‑1 (MQ=255)
                                      tCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCt                             <  1:3208637/71‑1 (MQ=255)
                                      tCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCt                             <  1:247978/71‑1 (MQ=255)
                                        accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                           >  1:116741/1‑71 (MQ=255)
                                        accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                           <  1:2972727/71‑1 (MQ=255)
                                        accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                           >  1:503024/1‑71 (MQ=255)
                                        accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                           >  1:1407526/1‑71 (MQ=255)
                                        accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                           >  1:2531894/1‑71 (MQ=255)
                                        accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtt                            <  1:777133/70‑1 (MQ=255)
                                        accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtt                            >  1:2934076/1‑70 (MQ=255)
                                        accacTGTTCAGCGCATCGCGCAGCAGTCCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                           <  1:513235/71‑1 (MQ=255)
                                         ccacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt                          <  1:917543/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)
                                          cacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTGATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                           <  1:1661228/69‑1 (MQ=255)
                                            cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGGaa                                    <  1:1769867/58‑1 (MQ=255)
                                            cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaa                                    <  1:552048/58‑1 (MQ=255)
                                            cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaa                                    <  1:570389/58‑1 (MQ=255)
                                            cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaa                                    <  1:3650997/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:1091904/1‑71 (MQ=255)
                                                 cAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg                  >  1:2418600/1‑71 (MQ=255)
                                                 cAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg                  >  1:1541899/1‑71 (MQ=255)
                                                  aGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg                  <  1:3886369/70‑1 (MQ=255)
                                                  aGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg                  <  1:1219284/70‑1 (MQ=255)
                                                       aTCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa            <  1:1516029/71‑1 (MQ=255)
                                                       aTCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa            <  1:4128196/71‑1 (MQ=255)
                                                       aTCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa            <  1:2269720/71‑1 (MQ=255)
                                                        tCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACt                      <  1:4097144/60‑1 (MQ=255)
                                                        tCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACt                      <  1:58292/60‑1 (MQ=255)
                                                          gcgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGa          >  1:1824521/1‑70 (MQ=255)
                                                          gcgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAc         >  1:3129345/1‑71 (MQ=255)
                                                          gcgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAc         >  1:1871080/1‑71 (MQ=255)
                                                          gcgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAc         >  1:2452422/1‑71 (MQ=255)
                                                           cgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAcg        <  1:3668847/71‑1 (MQ=255)
                                                            gcagcaGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAcg        <  1:1446628/70‑1 (MQ=255)
                                                              agcagTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa            <  1:2472017/64‑1 (MQ=255)
                                                              agcagTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa            <  1:595498/64‑1 (MQ=255)
                                                              agcagTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa            <  1:3724605/64‑1 (MQ=255)
                                                                 agTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTc  >  1:2332054/1‑71 (MQ=255)
                                                                  gTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa            <  1:3293700/60‑1 (MQ=255)
                                                                  |                                                                     
CCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTC  >  minE/2384423‑2384558

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: