Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,384,502 G→A 16.2% C353C (TGC→TGT metE ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,384,5020GA16.2% 284.9 / 41.2 111C353C (TGC→TGTmetE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (43/50);  new base A (5/13);  total (48/63)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 114 total aligned reads displayed.

CTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAGG  >  minE/2384437‑2384568
                                                                 |                                                                  
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:3489890/71‑1 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            <  1:3199369/71‑1 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            <  1:307096/71‑1 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            <  1:2969347/71‑1 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            <  1:3219195/71‑1 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            <  1:4081980/71‑1 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            <  1:2450233/71‑1 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            <  1:2357012/71‑1 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            >  1:1663265/1‑71 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            <  1:3406489/71‑1 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            <  1:1189413/71‑1 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            >  1:1458930/1‑71 (MQ=255)
   cACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            <  1:1266250/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:3858875/1‑71 (MQ=255)
             gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                                  >  1:1880534/1‑71 (MQ=255)
             gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                                  >  1:3786761/1‑71 (MQ=255)
             gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                                  >  1:472278/1‑71 (MQ=255)
               gCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg                                                            >  1:2982017/1‑59 (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:1334222/1‑69 (MQ=255)
                     gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                                            >  1:3576565/1‑69 (MQ=255)
                     gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                          <  1:2054522/68‑1 (MQ=255)
                     gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                          <  1:3006751/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)
                     gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                          <  1:354384/68‑1 (MQ=255)
                     gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                          <  1:2837048/68‑1 (MQ=255)
                      tgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttt                                                               <  1:3202546/49‑1 (MQ=255)
                      tgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCAttt                                                               <  1:3566495/49‑1 (MQ=255)
                        tCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCt                                       <  1:3208637/71‑1 (MQ=255)
                        tCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCt                                       <  1:247978/71‑1 (MQ=255)
                          accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                     >  1:2531894/1‑71 (MQ=255)
                          accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                     <  1:2972727/71‑1 (MQ=255)
                          accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                     >  1:116741/1‑71 (MQ=255)
                          accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                     >  1:1407526/1‑71 (MQ=255)
                          accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                     >  1:503024/1‑71 (MQ=255)
                          accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtt                                      <  1:777133/70‑1 (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: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)
                                 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)
                                    aGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg                            <  1:1219284/70‑1 (MQ=255)
                                         aTCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                      <  1:4128196/71‑1 (MQ=255)
                                         aTCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                      <  1:1516029/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:2452422/1‑71 (MQ=255)
                                            gcgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAc                   >  1:3129345/1‑71 (MQ=255)
                                             cgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAcg                  <  1:3668847/71‑1 (MQ=255)
                                              gcagcaGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAcg                  <  1:1446628/70‑1 (MQ=255)
                                                agcagTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                      <  1:3724605/64‑1 (MQ=255)
                                                agcagTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                      <  1:2472017/64‑1 (MQ=255)
                                                   agTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTc            >  1:2332054/1‑71 (MQ=255)
                                                    gTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                      <  1:3293700/60‑1 (MQ=255)
                                                       ccAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCAc         >  1:3628290/1‑70 (MQ=255)
                                                          gTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACt                                <  1:1079116/44‑1 (MQ=255)
                                                          gTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTc     <  1:3546531/71‑1 (MQ=255)
                                                           ttCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                      <  1:2733722/53‑1 (MQ=255)
                                                           ttCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                      <  1:2919233/53‑1 (MQ=255)
                                                           ttCATGACATTTTTGCAGGGCGAAGGGAAACCAGCTTTTCACtt                               <  1:313001/44‑1 (MQ=37)
                                                           ttCATGACATTTTTGCAGGGCGAAGGCAAACCAGCTTTTCACtt                               <  1:4089148/44‑1 (MQ=38)
                                                           ttCATGACATTTTTGCAGGGCGAAGGCAAACCAGCTTTTCACtt                               <  1:3960536/44‑1 (MQ=38)
                                                           ttCATGACATTTTTGCAGGGCGAAGGCAAACCAGCTTTTCACtt                               <  1:1579154/44‑1 (MQ=38)
                                                           ttCATGACATTTTTGCAGGGCGAAGGCAAACCAGCTTTTCACtt                               <  1:826648/44‑1 (MQ=38)
                                                             cATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACt                                <  1:915259/41‑1 (MQ=255)
                                                             cATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAgg  >  1:3582965/1‑71 (MQ=255)
                                                             cATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAgg  >  1:577107/1‑71 (MQ=255)
                                                             cATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAgg  >  1:3260301/1‑71 (MQ=255)
                                                             cATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAgg  >  1:3328887/1‑71 (MQ=255)
                                                             cATGGCATTTTTGTAGGGCGAAGGCAAACAAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAgg  >  1:34398/1‑71 (MQ=255)
                                                                 gCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAg   >  1:3576373/1‑66 (MQ=255)
                                                                 gCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAg   >  1:1334390/1‑66 (MQ=255)
                                                                 gCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAg   >  1:739878/1‑66 (MQ=255)
                                                                 |                                                                  
CTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAGG  >  minE/2384437‑2384568

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: