Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 232,790 G→T 5.6% intergenic (‑29/+141) hemB ← / ← ykiB porphobilinogen synthase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE232,7900GT 5.6% 349.0 / 2.8 102intergenic (‑29/+141)hemB/ykiBporphobilinogen synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (55/41);  new base T (4/2);  total (59/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

AGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGA  >  minE/232720‑232856
                                                                      |                                                                  
agaTTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAggg                                                                    <  1:4372453/71‑1 (MQ=255)
 gaTTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAgggg                                                                   >  1:252135/1‑71 (MQ=255)
   tttGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCa                                                                 <  1:3299892/71‑1 (MQ=255)
    ttGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAgggg                                                                   >  1:721864/1‑68 (MQ=255)
    ttGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCa                                                                 >  1:1408909/1‑70 (MQ=255)
    ttGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCa                                                                 >  1:1575050/1‑70 (MQ=255)
     tGCGCAGGCGACGAGGGCGTTGGATTAAGTCAGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCa                                                                 >  1:3778778/1‑69 (MQ=255)
        gcAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCAt                                                                >  1:3362738/1‑67 (MQ=255)
         cAGGCGACGAGGGCGTTTGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGtat                                                           <  1:3653228/71‑1 (MQ=255)
          aGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGtata                                                          >  1:2583104/1‑71 (MQ=255)
                  aGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCAt                                                                <  1:2820226/57‑1 (MQ=255)
                  aGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCAt                                                                <  1:338814/57‑1 (MQ=255)
                    ggCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAg                                                   >  1:96774/1‑68 (MQ=255)
                    ggCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAg                                                   >  1:1443165/1‑68 (MQ=255)
                           gATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGta                                                            >  1:3599018/1‑52 (MQ=255)
                           gATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGta                                                            >  1:4635853/1‑52 (MQ=255)
                           gATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGta                                                            >  1:2663139/1‑52 (MQ=255)
                           gATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGta                                                            >  1:2550627/1‑52 (MQ=255)
                           gATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGta                                                            >  1:2365672/1‑52 (MQ=255)
                           gATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGt                                                             >  1:1202689/1‑51 (MQ=255)
                           gATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAggg                                         <  1:4490560/71‑1 (MQ=255)
                              tAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATTAAGGTGCACAGTATAcc                                                        >  1:1522775/1‑53 (MQ=255)
                              tAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATTAAGGTGCACAGTATAcc                                                        >  1:1002260/1‑53 (MQ=255)
                              tAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATTAAGGTGCACAGTATAcc                                                        >  1:2110716/1‑53 (MQ=255)
                              tAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATTAAGGTGCACAGTATAcc                                                        >  1:4113251/1‑53 (MQ=255)
                              tAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATg                                      >  1:3559308/1‑71 (MQ=255)
                              tAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATg                                      >  1:2864531/1‑71 (MQ=255)
                              tAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATg                                      >  1:2784385/1‑71 (MQ=255)
                              tAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATg                                      >  1:2242289/1‑71 (MQ=255)
                               aaGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATg                                      <  1:2757172/70‑1 (MQ=255)
                               aaGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGt                                     >  1:257502/1‑71 (MQ=255)
                               aaGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGt                                     >  1:542502/1‑71 (MQ=255)
                               aaGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGt                                     >  1:759860/1‑71 (MQ=255)
                               aaGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGt                                     >  1:4366952/1‑71 (MQ=255)
                               aaGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGt                                     >  1:48067/1‑71 (MQ=255)
                               aaGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGt                                     >  1:1113323/1‑71 (MQ=255)
                               aaGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGAGGt                                     >  1:3705879/1‑71 (MQ=255)
                               aaGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAAGGATGt                                     >  1:2254463/1‑71 (MQ=255)
                                aGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGt                                     <  1:1426285/70‑1 (MQ=255)
                                 gTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGTGCATAGTATACCTGAAGCAGGGTAGGGATGttt                                   <  1:2206511/71‑1 (MQ=255)
                                 gTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGttt                                   <  1:4515843/71‑1 (MQ=255)
                                  tCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGtttt                                  <  1:693487/71‑1 (MQ=255)
                                  tCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGtttt                                  <  1:456355/71‑1 (MQ=255)
                                  tCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGtttt                                  <  1:1269927/71‑1 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGCAGCAGGGTAGGGATGTTTTAc                                >  1:701344/1‑69 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGttt                                   <  1:3695294/66‑1 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTAc                                <  1:3549365/69‑1 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACg                               >  1:2531425/1‑70 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACg                               >  1:4467977/1‑70 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACg                               >  1:217451/1‑70 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACg                               >  1:4462430/1‑70 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACg                               >  1:883900/1‑70 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGa                              <  1:4033688/71‑1 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGa                              <  1:4233431/71‑1 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGa                              <  1:4450116/71‑1 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGa                              <  1:144013/71‑1 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGa                              <  1:1405090/71‑1 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGa                              <  1:4715860/71‑1 (MQ=255)
                                      tCATGGTCTGCCTGATGTTTGTGGAATCCAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGa                              <  1:3847946/71‑1 (MQ=255)
                                        aTGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGa                              >  1:404155/1‑69 (MQ=255)
                                        aTGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGa                              >  1:2841829/1‑69 (MQ=255)
                                        aTGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGa                              >  1:1639764/1‑69 (MQ=255)
                                         tGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGAt                                       <  1:2981103/59‑1 (MQ=255)
                                         tGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGAt                                       <  1:583889/59‑1 (MQ=255)
                                         tGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAg                           <  1:4199204/71‑1 (MQ=255)
                                          ggTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTa                                 >  1:1924876/1‑64 (MQ=255)
                                           gTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGt                                     >  1:287725/1‑59 (MQ=255)
                                           gTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGt                                     >  1:2762000/1‑59 (MQ=255)
                                           gTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGtt                         >  1:4518488/1‑71 (MQ=255)
                                           gTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGtt                         >  1:299825/1‑71 (MQ=255)
                                           gTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGt                          >  1:2693535/1‑70 (MQ=255)
                                            tCTGCCTTATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTg                        <  1:781476/71‑1 (MQ=255)
                                            tCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGtttt                                  <  1:4154612/61‑1 (MQ=255)
                                            tCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTg                        <  1:3805633/71‑1 (MQ=255)
                                            tCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTg                        <  1:632972/71‑1 (MQ=255)
                                             cTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGaaa                            >  1:3487748/1‑66 (MQ=255)
                                             cTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGaaa                            >  1:1179273/1‑66 (MQ=255)
                                                ccTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAg                                                   >  1:2619142/1‑40 (MQ=255)
                                                ccTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTcgc                    <  1:1160607/71‑1 (MQ=255)
                                                ccTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTcgc                    <  1:3896008/71‑1 (MQ=255)
                                                ccTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTcgc                    <  1:3822706/71‑1 (MQ=255)
                                                 cTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTcgcg                   <  1:1052513/71‑1 (MQ=255)
                                                   gATGTTTGTGGAATCGAGGGGCATAGTATTCCTGAAGCAGGGTAGGGATGTTTTACGa                              <  1:1887723/58‑1 (MQ=255)
                                                   gATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGAt                 >  1:3613690/1‑71 (MQ=255)
                                                    aTGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATg                >  1:519332/1‑71 (MQ=255)
                                                    aTGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATg                >  1:4080572/1‑71 (MQ=255)
                                                      gTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGt               <  1:834206/70‑1 (MQ=255)
                                                       tttGTTGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGtt                         <  1:936093/59‑1 (MQ=255)
                                                         tgtgGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGa            >  1:3296178/1‑70 (MQ=255)
                                                         tgtgGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGa            >  1:2796822/1‑70 (MQ=255)
                                                         tgtgGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGa            >  1:2398010/1‑70 (MQ=255)
                                                         tgtgGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGa            >  1:2098810/1‑70 (MQ=255)
                                                             gAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACaaga       <  1:2155618/71‑1 (MQ=255)
                                                               aTCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACaag        >  1:715872/1‑68 (MQ=255)
                                                                tCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTc                      <  1:1210722/53‑1 (MQ=255)
                                                                tCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACaagaa      >  1:403665/1‑69 (MQ=255)
                                                                tCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAagaaga    <  1:421203/71‑1 (MQ=255)
                                                                tCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAagaaga    <  1:856505/71‑1 (MQ=255)
                                                                  gAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAgaga  <  1:1007609/71‑1 (MQ=255)
                                                                    ggggCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTc                      <  1:4368246/49‑1 (MQ=255)
                                                                      |                                                                  
AGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGA  >  minE/232720‑232856

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